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PREFACE 

The aim of this hook is to give a sound, systematic, unmistakably 
clear and reasonably complete treatment of the mathematical and 
mechanical aspects of the gyroscope and its more important 
applications. Since the book is concerned with fundamental 
principles, constructional details are not entered into except inso- 
far as they may aid in clarifying the principles. Friction has been 
left out of consideration in the motion of the gyroscope itself, 
because the construction, mounting and casing of modern gyro- 
scopes are such as to keep friction of all kinds as low as possible. 
Unavoidable friction is taken care of by the driving motor. 

Because I think the mathematical theory of the gyroscope and 
its applications are best treated with the aid of vectors, I have used 
vector methods throughout the book. And since a knowledge of 
vector analysis is not assumed on the part of the reader, the first 
chapter is devoted to the exposition of the amount of vector 
analysis needed in the subsequent chapters. 

The right-handed system of coordinate axes is used throughout 
the book, because of its advantages for this work. 

A complete and exact mathematical treatment of gyroscopic 
motion becomes intractable almost at the beginning. Approxima- 
tions of minor importance must be made in order to obtain 
tractable and solvable equations. When making simplifying 
approximations, I have pointed out their nature and in some 
cases have shown by numerical examples that the errors thus 
introduced were of no consequence. 

In the preparation of this book I have consulted the works of 
many previous writers, the most important of which are listed in 
the Bibliography at the end of the book. 

It is a pleasure to record my thanks and obligations to the 
Sperry Gyroscope Company, Great Neck, New York, for their 
unstinted cooperation in furnishing information and photographs 
relating to various gyroscopic instruments and applications, I also 
wish to record my thanks to the following other manufacturers for 
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furnishing information and photographs: The Arma Corporation, 
Garden City, New York; The Minneapolis-Honeywell Regulator 
Company, Minneapolis, Minnesota; and The National Engineering 
Company, Chicago, Illinois. Finally, I wish to thank Mr. Donald 
Trumpy, of the yacht-building firm of John Trumpy and Sons, 
Annapolis, Maryland, for gyroscopic data concerning a 29,000- 
mile cruise which he took on a yacht equipped with a Sperry ship 
stabilizer. 

J. B. SCARBOROUGH 
August, 1967 
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PART I 

THEORY OF THE GYROSCOPE 



CHAPTER I 
Some Necessary Vector Analysis 

The theory, behavior and applications of the gyroscope can be 
explained best by means of vectors. We therefore devote the present 
chapter to the exposition of the needed amount of vector analysis. 

1. Scalar and vector quantities 

The quantities which occur in physics are of two kinds: those 
which have magnitude only and those which have both magnitude 
and direction. The former are called scalar quantities and the latter 
are called vector quantities . Examples of scalar quantities are density, 
temperature and electric potential. Familiar examples of vector 
quantities are force, velocity and acceleration. Vectors may also be 
used to denote position, in which case they are calledposition vectors. 

Vector quantities are represented geometrically by segments of 
straight lines, the line segments carrying arrow heads at one end to 
indicate the sense of the vector. The magnitude and direction of the 
directed quantity are indicated by the length and orientation of the 
vector. 

In Fig. I is shown a vector AB, which is also denoted by the 

single letter r. The projections ofAB on the coordinate axes are the 
vectors X, Y and Z as shown. 

A vector is designated in print by a single boldface letter, as r in 
Fig. 1 ; or by the letters designating its end-points, with an arrow 

written over the two letters, as AB in Fig. 1. 

Vectors may be multiplied or divided at will by any numbers or 

scalars, the results always being vectors. Thus 3r, ABjl, etc., are 
vectors having the same direction and sense as the original vectors. 
However, the multiplication or division of a vector by a negative 
number always reverses the sense of the vector. 

Two vectors are equal when they have the same magnitude, 
direction and sense. 
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Fig. 1 

2, Geometric addition and subtraction of vectors 

To find the geometric sum of two vectors, place the initial point 
of the second vector at the terminal point of the first and then draw 
a line from the initial point of the first vector to the terminal point 
of the second. The vector thus drawn is the geometric sum of the 
given vectors. In Fig. 2, for example, AC is the geometric or vector 
sum of AB and EG. Such addition is denoted by either of the 
equations jt, ^n jfr (e> n 

jfi JLJ "T J-J O == Xt \j 3 \& */ 

p 4. = R f2*21 

The geometric difference between two vectors is found by 
changing the sign of the subtrahend vector and then adding it to the 
other vector. Thus P O = P + ( 0) (2*3) 



Reference to Kg. 2 will show that vectors are added and sub- 
tracted by the parallelogram law, the sum being given by one 
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diagonal of the parallelogram and the difference by the other 
diagonal. 

The geometric sum of any number of vectors is found by the same 
procedure as in the case of two vectors; that is, the initial point of 
the second vector is placed at the terminal point of the first, the 




Fig. 2 




initial point of the third is placed at the terminal point of the 
second, etc. The geometric sum of all the vectors is the vector 
drawn from the initial point of the first to the terminal point of the 
last, as indicated in Kg. 3. 

3* Analytical addition of vectors 

The sum or resultant of several vectors is found analytically by 
first resolving the vectors into rectangular components along 
coordinate axes, finding the algebraic sums of these components 
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(projections on the coordinate axes), and then finding the resultant 
sum in magnitude and direction by the methods of analytic 
.geometry. Thus, if the vectors P l3 P 25 P 3 , ..., P n make the angles 
%? AJ TiJ a 2J A> 72 > * a %> Pn> 7n ^^ ^ ne coordinate axes, we have 

SX = P I cos % + P 2 cos & 2 + . . . 4- P n cos a tt , 

(3-1) 




Then the magnitude of the resultant is 

M = V{(^) 2 + ( s Y 
and its direction cosines are 

27 



, x 
(3-3) 



4. The unit coordinate vectors i, j, fc 

For many purposes it is convenient and advantageous to express 
vectors in terms of unit vectors along the coordinate axes. Let i 
denote a vector of unit length along the #-axis, let j denote a unit 
vector along the y-axis, and let k denote a unit vector along the 
2-axis, as shown in Fig. 4. Let r be a position vector drawn from the 
origin to the point P(a?,y,s). Then the projections of r on the 
coordinate axes are the vectors xi, yj and zk; and since a vector is 
equal to the geometric sum of its projections (Art. 2 5 Fig. 3), we 
may write 

(4-1) 



r i> r & r s are the magnitudes of the projections of r 
on the coordinate axes, then r may be written in either of the forms 



(4-2) 
(4-3) 
The vector r is sometimes written in the form [r, r , r 



, 



Let us now consider two vectors p and q. They may be expressed 
in the forms 
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where p^ p 2 , p 3 and g l3 g 2 , g 3 are the magnitudes of the rectangular 
components of p and q, respectively. Adding p and q s we get 

p + q = (p l + q-j) i + (#2 + #2)] + fe + ?s) ^- (^*4) 

The geometric sum of any number of vectors may likewise be 
expressed in terms of the unit coordinate vectors. 




Fig. 4 



5. The scalar or dot product of two vectors 

We have already pointed out that a vector may be multiplied or 
divided by abstract numbers and scalars at pleasure, the resulting 
product or quotient always being a vector along the line of the 
original vector. But the product of one vector by another is an 
entirely different matter and must be explicitly defined before 
proceeding further. 

There are two distinct kinds of vector multiplication: scalar 
multiplication and vector multiplication. The scalar or dot product 
of two vectors a and b is indicated and defined as follows: 

where 6 is the angle between a and b (Fig. 5). The dot product is 
thus a scalar. 
Since a cos 6 is the length of the projection of a on b, and b cos 6 
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Is the length of the projection of b on a, the scalar product of two 
vectors is equal to the length of either vector multiplied by the length of 
the projection of the other upon it. 

If the vectors a and b are perpendicular to each other, then 
6 = 90, cos = 0, and a b = 0. Hence the condition for perpendicu- 
larity of two vectors is that their scalar product be zero. 




Fig. 5 

If a and b are parallel, then 6 = 0, cos#=l and a*b = a&. In 

particular, 

r a-a=& 2 . (5-2) 

From (51) it is plain that 

a-b=b.a. (5*3) 

The order of the factors in the dot product is thus immaterial. 

The following relations among the dot products of the unit 
coordinate vectors are of great importance and should always be 
k.ptinmind: j.j.j.j _ k . k =I , 

M-l.k-j.k-J (5 ' 4) 

We shall now show that the distributive law of multiplication 
holds for the scalar product of two vectors. Consider a vector a 
and two other vectors b and c. We are to show that 

a-(b-hc) = a-b + a-c. (5-5) 

Proof. In Fig. 6 consider the vector O4 = a and the projections 
of the other vectors upon it. We have 



(a) 
(b) 
(c) a.c=(Qj) (GH). 
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Adding (6) and (c), we get 



(d) a . b + a . c = (OA) ( FG + GH) = (OA) (FH). 
Since the right-hand members of (a) and (d) are the same, we have 



a (b + c) = a - b + a - c. 
c 



Q.E.D. 



(5-5) 




H 



Fig. 6 

This result and the relations among the unit coordinate vectors 
enable us to write the product a b in terms of the rectangular 
components of a and b. We have 

a - b = (%! + a 2 j + a 3 k) foi -f 6 2 j + 6 3 k) 



by (5-4), 

all cross-products such as %& 2 i- j being zero. We thus have the 
important relation - _ . - , 

6. The vector or cross product of two vectors 

The vector product of two vectors a and b is indicated and defined 
as follows: 



a x b = (ab sin 6) n, 



(6-1) 



10 
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where n Is a unit vector perpendicular to the plane of a and b, and 6 
is the angle through which a must be rotated to make it coincide 
in direction with b (Fig. 7). The vector product of two vectors is 
thus a vector perpendicular to the plane of the given vectors and of 
magnitude equal to the area of the parallelogram having the two 
given vectors as adjacent sides. 



axb 




Fig. 7 

The sense (indicated by the arrowhead) of the vector a x b is 
determined by the following rule, called the right-handed-screw rule: 

The product vector axb must be drawn in the direction a right- 
handed screw (corkscrew, for example) would advance if it were 
rotated in the direction that a would have to be rotated to make it 
coincide in direction with b. Thus, a glance at Fig. 7 shows that a 
right-handed screw would advance upward when rotated in the 
direction required to bring a into coincidence with b. 

Let us now consider the product b x a. In this case a right-handed 
screw would advance downward if rotated in the direction required 
to bring b into coincidence with a. Furthermore, the angle 6 in 
this case would be negative; and since sin ( - 6) = - sin 0, we have 

bxa = [6asin(~^)]n=-(6asin(9)n=-(axb), (6-2) 
thus agreeing with the corkscrew rule. 

Equation (6-2) shows that the multiplication is not commutative 
in finding the vector product of two vectors. The order of the factors 
must be observed. If two factors are interchanged, the sign of the 
product must be changed. 
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Prom (6-1) it is plain that if a x b = 0, a and b are parallel. In 
particular, axa = () (6 . 3) 

Since the unit coordinate vectors are mutually perpendicular to 
one another (0 = 90), we get from (6-1), (6*2) and (6-3), the following 
relations: . . - . . 



(6-4) 



jxk=i= -kxj, 
kxi=j= -ixk, 
ixi=j xj = kxk 

We now proceed to find the rectangular components of the 
product vector a x b. Let us denote this vector by c, and let the 
magnitudes of its rectangular components be c v c 2 , c 3 . Then, since 
c is perpendicular to both a and b, we have: 

a - c = %<?! + 2 c a + a 3 c 3 = 0, 
Solving these equations for c and c 2 in terms of c 3 , we get 

C-j = r- r 5 C 2 == f 7 . (6.5) 

Since c 2 = cf + c| + cf , we have 

/*2_* 



6! 6 8 -2a 2 

Now adding and subtracting af&f + a|6| + a|6| within the brackets 
of the right-hand member, we obtain 



which can be written in the form 



12 
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or 



Hence 



or 



c= 



(6-6) 




Fig. 8 

Now since c = cn = axb = (a6sin^)n j we get c = absmd. Sub- 
stituting this in (6-6), we get 



= + 1. 



(6-7) 



To determine which of the two signs on the right should be used, 
we utilize the fact that all the above relations hold for any values 
and positions of a and b. Hence we may take a and b as unit vectors 
along the positive x- and y-axes, respectively. Then c becomes a 
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unit vector along the positive g-axis. Hence the rectangular 
components of a, b and c in this case are 



Cl = 0, c 2 = 0, c 3 =l. 

Substituting these values in the left-hand member of (6*7), we find 

c 



that 






^^&g a 2 6 l5 



Hence 

and then (6*5) gives 



The product vector a x b is therefore 



This can be written more compactly in determinant form 

i j k 
axb = 



(6-8) 



% a 2 % 
&i & 2 6 3 

We are now in a position to prove that the distributive law of 
multiplication applies to vector products; that is, that 

a x (b + c) = (a x b) + (a x c). (6-9) 

To prove this extremely important fact we utilize the addition 
theorem of determinants, namely, 



&! 



&! 



Since the rectangular components of b -f c are b^ + c lt & 2 + C 2> 
we have, from (6-8), 



a x (b + c) = 



J 



i j k 



&i 6 2 



i j k 



= (axb) + (axc) by (6-8). 
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It is worth while to point out here that terms may be transposed 
from one side of a vector equation to the other, but that cancellation of 
a vector out of the two sides of the equation is not permissible. For 

example, suppose we have the vector equation 

a-c=b-c. 

We might be tempted to divide throughout by c and get a = b. But 
we get another solution by transposing b . c to the left-hand side 
and factoring out c, as we are permitted to do because the scalar 




c 
Fig. 9 




Fig. 10 



product obeys the distributive law of multiplication. Then the 
equation becomes 

(a-b)-c = 0. 

This equation is satisfied if a = b, but it is also satisfied if a =)= b and 
c is perpendicular to the vector a b. The situation in the latter case 
is represented by Fig. 9. 
In like manner, suppose we are given the vector equation 

axc = bxc. 

The left-hand side of this equation is a vector which is perpendicular 
to both a and c, and the right-hand member is a vector which is 
perpendicular to both b and c. Since the vectors axe and b x c 
are equal, the planes to which they are perpendicular must coincide 
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and therefore the three vectors a, b and c all lie in the same plane. 
Let us now transpose both members to the same side of the equa- 
tion and factor out c (permissible because the vector product obeys 
the distributive law). Then we have either 

(a-b)xc = or (b a)xc = 0. 

These equations are satisfied if a = b, but they are satisfied equally 
well if a =(= b and c is parallel to the vector a b. Fig. 10 represents 
the situation in this latter case. 



7. The vector triple product a x (b x c) 

This vector product is a vector lying in the plane of b and c, and 
is of great importance in the 
study of the gyroscope. Fig. 11 
represents the relative positions of 
the vectors involved in this triple 
product. The fact that the triple- 
product vector lies in the plane of 
b and c can readily be seen when 
it is remembered that the vector 
b x c is perpendicular to each of the 
three vectors b, c and a x (b x c) 
at the point 0. Hence all three of 
these vectors lie in the plane to 
which b x c is perpendicular at 0. 

It is to be noted that 




Fig. II 



ax (bxc)=f=(axb)xc; 

for the vector on the left-hand 

side lies in the plane of b and c, 

whereas the vector on the right-hand side lies in the plane of a and 

b; and these planes cannot be the same.* This fact shows that the 

vector triple product does not obey the associative law of 

multiplication. 

The vector triple product can be expressed as the difference of 
two vectors, as we shall now show. In Art. 6 the magnitudes of the 
rectangular components of the vector b x c were found to be 



* Since a, b, and c represent any three vectors in space, they are in 
general noncoplanar. 
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&2 c s - ^s c 2> &s c i - &i c 3 an( i &i c 2 "" Mr Hence by (6'S), we may write 

the triple product in the form 

i j k 

ax(bxc) = % 2 a 3 



Expanding this determinant and rearranging the terms slightly, 
we get 

a x (b x c) = BjfajCa + %c 3 ) i + 6 2 (% c i + % c s) J 



a 2 6 2 ) k. 

The form of the parenthetical expressions suggests that a third 
product term be added to each of the positive components and 
subtracted from the corresponding negative component. We there- 
fore add b^^i, b 2 a 2 c 2 j y 6 3 a 3 c 3 k to the positive components and 
subtract these terms from the corresponding negative components. 
We then have 

a x (b x c) = &!(%(?! + a a c a + %c 3 ) i -h & 2 K c i +&&+ 



ISTow reducing the parenthetical expressions by means of formula 
(5*6), we get 



or, finally, ax (bxc) = (a-c)b (a-b)c. (7*1) 

This is one of the most important formulas of vector analysis. 
Since the dot products in the parentheses are scalars, the right-hand 
member of (7-1) is the difference of two vectors. 

8. Differentiation of vectors 

If a vector r is a function of a variable parameter t (time, for 
example), then an increment A in t will cause a change Ar in r. 
The new value of the vector is therefore r + Ar, as indicated in 
Kg. 12. Then the derivative of r with respect to t is defined to be 

dr .. Ar 

-=T 

at 



SOME NECESSARY VECTOB ANALYSIS 17 

If r =rri-h#3 +zk s then 

dr __ . dx .dy dz 

~di =l ~dt +j ~di + dt> (8-1) 

' where i, j and k have been treated as constants in the differentiation. 
More generally, , , ,, , 

The derivative of a vector is thus a vector. 




Fig. 12 

To find the derivative of the scalar product a * b we write 
Then 



and 



AP 

_ 



Ab 

_ 



Aa 

.__ 



Ab 

._. 



_ 
~dt~~ Q '"dt + 



da 



or 



d(a-b) 

v _ : 

dt 



, 

__ L . _ 

dt dt ' 



(8-3) 



The derivative of the vector product a x b is found by a similar 
procedure. We write V = a x b 



and 



Ab + Aa x b + Aa x Ab . 
. ' . A V = a x Ab H- Aa x b + Aa x Ab, 

AV Ab Aa . , Ab 
= ax + x 
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dV db da 
Hence _ = ax _ +wxb? 

d(axb) db da ( . 

JJ.ax^+^xb (8-4) 

db K da 

==: a x r ox - 7 T~ 

d * 

Note that when the order of the factors is changed in a vector 
product the sign before the product must also be changed. 

If a vector is constant in magnitude, its derivative is not neces- 
sarily zero; for the vector is still free to change in direction. Sup- 
pose, for example, that 

r = a constant. 

Then r 2 = a constant. 

That is, r - r = a constant. 

Differentiating this, we have 

dr dr A 

r -dt +r -Tt =0 > 

or r * dr = 0. 

This shows that dr is perpendicular to r. 



CHAPTER II 
Some Fundamental Principles of Mechanics 

This chapter deals with certain topics in. mechanics which, are 
fundamental for the study of the gyroscope. 

9. Velocity and momentum 

Assume that a particle of mass m is moving along a curve in three- 
dimensional space, and let r denote a position vector drawn from 
a point to the particle m. Then r is evidently a function of the 
time, or r = r(). Let v denote the velocity of the particle along its 

path. Then 
* 



The momentum of the particle is mv=mt and is thus a vector 
quantity having the same direction and sense as v. 

10. The fundamental equation of dynamics 

The most fundamental and important principle of dynamics is 
Newton's second law of motion. It states that 



where F is the resultant of all the forces acting on the particle. If 
m is constant, as is usually the case, the above equation becomes 

F=w^7==m-7-==mf==ma, (10-2) 

at at 

where a denotes the acceleration imparted to the particle by the 
force F. 
Since V= 



and 
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equation (10*2) is the vector equivalent of the three Cartesian 
equations 

-C3L - iUitfJ-) I 

(10-3) 




where X, J and Z are the scalar components of F along the 
coordinate axes, 

11. D'AJemBert's principle 

In the fundamental equation Fm&, the quantity m& is called 
the effective force on m. The quantity m( a) is called the reversed 
effective force or inertia force. D'Alembert's principle is based on 
Newton's second and third laws of motion and states that: 

The inertia force is in equilibrium with the external applied force, 

or 



This principle has the effect of reducing a dynamical problem to a 
problem in statics and may thus make it easier to solve. Centrifugal 
force is a familiar example of an inertia force. 

12. Moment of a force about a point 

Let denote a point in three-dimensional space, let F denote a 
force, and let r be a vector drawn from to the point of application 
of F, as indicated in Fig. 13. Then the magnitude of the moment 
of F about is defined to be the product of the force by the perpen- 
dicular distance from to the line of action of F, or 

magnitude of moment = (r sin 6) F. 

The moment can therefore be represented by the vector product 
r x F, and we thus have the vector equation 

M = rxF, (12-1) 

where M denotes the moment of F about 0. The moment vector M 
is therefore perpendicular to the plane of r and F and pointing in 
the direction a right-handed screw would advance when rotated in 
the direction required for rotating r into the direction of F. In 
Kg. 13 the vector M would be directed vertically upward if the 
plane of r and F is horizontal. 
Let us take as the origin of a system of rectangular coordinates. 
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Fig. 13 




Fig. 14 
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Then x, y, z are the scalar components of r, Z, 7 5 Z are the scalar 
components of F, and M x , My, M z are the scalar components of M. 
Then the moment of F about the origin is, by equation (6-8), 

' i j k 



X 

X 



y 

j 



z = 



Hence 




(12-2) 



(12-3) 



The quantities M X9 M y) M z are the magnitudes of the moments of F 

about the x-, y- and 2-axes, re- 

spectively, as indicated in Fig. 14. 

It is clear that the component X 

can have no moment about the 

ru-axis. Viewed from the point 0, 

the moment of Z about the #-axis 

is positive and the moment of T 

is negative, because the moment 

yZ would tend to produce rotation 

to the right and the moment zT 

would tend to produce rotation 

to the left. Hence the resultant 

moment about the #-axis is 

yZ zY. It is to be noted that the 

moment of F about each axis is 

a vector along that axis. It is true 

in general that the moment of a force 

about any axis is a vector along 

that axis. 

13. Angular velocity 

Angular velocity is a vector 
quantity and is represented by a 

vector along the axis of rotation and pointed in the direction 
required by the right-handed-screw rule. 

Consider a rigid body rotating with angular velocity to about an 
instantaneous axis passing through a fixed point 0, as indicated in 
Fig. 15. Let P denote any point in the body, and let r be a position 
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vector drawn from to P. The linear speed of P is evidently 
(r sin 6) o). The linear velocity v of P is always perpendicular to the 
plane OMP* and is therefore given by the vector equation 

j k 
v = o)xr = i w, w, w, (13-1) 




where a) l9 <w 2 , &> 3 are the scalar components of to, Hence the scalar 
components of v along rectangular axes with origin at are 



(13-2) 



14. Moment of momentum of a rigid body about a fixed point 

Returning to Fig. 15, let us consider a particle of mass m at the 
point P. The momentum of m is wv, directed along the tangent to 
the path at P. The moment of this momentum about is, by (12*1), 

h = rxmv=mr x (to xr). 

This vector, it will be observed, is perpendicular to both r and v. 
It therefore lies along M P and is directed toward M. The moment 
of momentum of the entire body about is therefore 

H = Srxwv=Swrx(toxr) 

= Sm[co(r r)-r(r to)] by equation (7-1), 
or H = (Smr 2 )to~Swr(r.to). (14-1) 



* The proof that v is perpendicular to the plane OMP and therefore 
perpendicular to OP is as follows: 

For brevity, let G denote the plane of the circle described by P. Then 
OA, the axis of rotation, is perpendicular to C. The plane OMP is also per- 
pendicular to C; for if a line is perpendicular to a plane, every plane passing 
through it is perpendicular to that plane. 

Now v is directed along the tangent to the circle at P and therefore lies in 
the plane C. Moreover, it is perpendicular to the radius MP, which is also 
the line of intersection of the planes C and OMP. Hence v is perpendicular 
to OMP and therefore perpendicular to OP; because if two planes are 
perpendicular to each other, a straight line drawn in one of them perpendicular 
to their intersection is perpendicular to the other plane. 
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This equation shows that H Is the geometric difference of two 
vectors, one along <*> and the other along r. 

If we take a set of rectangular axes with origin at 0, then 



Then (14-1) becomes 

H = Sm(# 2 + y 2 + z 2 ) K i -f w 2 j 4- a> B k) 



(14-2) 



Hence the scalar components of H are 



H x = 



tf^ie 2 + 2/ 2 4- ^ 2 ) 



(14-3) 



The quantities Sm(2/ 2 4-^ 2 ), Sm(o; 2 +2; 2 ) and Sm(^ 2 + y 2 ) are the 
moments of inertia of the body about the x-, y- and 2-axes, respec- 
tively; and the quantities Swyz, %mzx and lumxy are the products of 
inertia about the same axes. 

Following the usual custom, we put 



(14-4) 



Then (14-3) becomes 




(14-5) 
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If the rigid body possesses three mutually perpendicular planes 
of symmetry which intersect at the origin, the products of inertia 
are zero. To see the truth of this statement, assume that the 
?/2-plane is a plane of symmetry. Then for each product such as 
mxy there is a symmetrically placed product m( x)y on the other 
side of the plane, and all these pairs cancel each other in the sum- 
mation. A similar result holds for each of the other two planes of 
symmetry. 

The lines of intersection of the three planes of symmetry are 
called the principal axes of inertia of the body. Hence, when the 
coordinate axes are the principal axes of inertia of the body, there 
are no product-of4nertia terms in equations (14-5) and these 
equations then become 

H x =Ao) v H y =Boj^ H 3 = Ca) 3 , (14-6) 

and therefore 

fi- = V(51 + ^+J?l) = V(^M+"B^I + ^I). (U-7) 

15. The equations of motion of a rigid body about a fixed point 

Differentiating the equation H =S(r x mv) with respect to time, 
we have 



=S[(r x ma) -f- (v x mv)]. 
But ma = F and v x v = 0. Hence the above relation becomes 

= S(rxF) = M by (12-1) 



where M denotes the sum of the moments about of all the 
external forces acting on the body. 
From (15'1) and (12-3), we get the following equations: 



^=I>(zX-xZ)=M v , } (15-2) 

dt 
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The vector form of these equations is 



and it refers to fixed space. 

Equations (15*1) and (15-2) are the mathematical statements of 
a most important dynamical law: 

The time rate of change of the angular* momentum of a rigid body 
rotating about any axis is equal to the moment of the applied external 
forces about the same axis. 

Since dM/dt is the time rate of change of H 5 it may be regarded 
as the velocity of the end of the vector H. Equations (15*1) and 
(15-2) therefore state: 

The velocity of the end of the angular-momentum vector is equal to 
the moment of the external forces. 

This statement of the law is sometimes convenient to use in 
dealing with the gyroscope. 

If no external forces are acting on the body, or if there are such 
forces but they have no moment about the axis of rotation, then 
M = and (16-1) becomes 



Hence H = a constant. (15-3) 

Note that in this case the angular-momentum vector remains 
constant in magnitude and fixed in direction. 

Equation (15-3) expresses the law of the Conservation of Angular 
Momentum. 

Furthermore, when H is constant equation (14-7) gives 

A*u\ + B*(i)\ -f O 2 ^| = a constant. (15-4) 

16. The general differential equations of motion. Baler's equations 

In the previous article the equations of motion were referred to 
a set of axes fixed in space. In order to study the general motion of 
a rigid body about a fixed point, however, it is necessary to employ 
two sets of rectangular axes, a fixed set and a moving set, but with 
both sets having the fixed point as origin. Let (Pig. 1 6) denote the 

* Moment of momentum is also called angular momentum and will often 
be called that in this book, because the name ' angular ' is less clumsy than 
* moment of. 
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fixed point, let X I 7 1 Z 1 denote a set of fixed axes (fixed in 
space), and let XYZ denote another set of rectangular axes 
which, moves as a rigid trihedral framework about in any manner. 
Let i, j, k denote unit vectors along OX, OT and OZ 3 respec- 
tively; let co denote the angular velocity of the moving trihedron 




Fig, 16 

with respect to (as viewed from) the fixed set of axes ; and let H 
denote the angular momentum of the body with respect to the 
fixed axes. Then if o) x , o) y , o) z and H x , H y , H z denote the magnitudes 
of the projections of to and H on the axes of the moving trihedron, 

wehave (16-1) 

. (16-2) 



Although i, j, k are unit vectors along OX, OY, OZ and are 
constants with respect to those axes, they are not constants with 
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respect to the fixed axes - X^Z^ Their directions change with 
respect to (as viewed from) the fixed axes and hence vary with the 
time. The time derivative of H is therefore from (16-2): 



To evaluate the derivatives dijdt, djjdt, dkfdt we regard i, j, k as 
position vectors drawn from 0. If r is a position vector drawn from 
to any point in a rotating body, then by (13-1) 



dr 



3 



Now applying this formula to the vectors i, j, k in turn, we have 

i J k 



di 



100 
i J k 



0) 



dk 



010 
i J k 



001 
On substituting these derivatives in (16-3) we get 

dH 

dt'' 



But by (15-1) 



dtt 

dt '' 
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On equating these two values ofdtl/dt, we get 
-^+H s o> v -H y a> z =M x , 

(16-4) 



These are the differential equations of motion referred to the moving 
axes, that is, referred to the axes of the moving trihedron. They are the 
most important differential equations connected with gyroscopic 
theory. 

If the axes OX, Y, OZ coincide with the principal axes of inertia 
of the body for the point 0, then 

H x =Aa) x , H y ^Boj y) H,=Ca) 8) 
and therefore 

dH x= A do>x ^ = ft^H dHg r do) g 
dt dt' dt dt' dt dt ' 

Substituting these into (16-4), we get 



(16-5) 



If the rotating body is dynamically symmetrical about the axis 
OZ, then B= A and equations (16-5) become 



(16-6) 



Equations (16-5) and their special cases (16-6) are known as Euler's 
dynamical equations for the motion of a rigid body about a fixed point. 
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Here we make a short digression to point out an interesting 
phenomenon. Assume that we have a solid of revolution which can 
be made to spin in a vertical position about its geometric axis. Then 
its differential equation of motion would be the third equation of 
(10-6). Suppose the solid is given a vigorous spin and let go. Then 
its equation of motion after that instant would be 



from which Ca) z =a constant =&, say, 

k 



or 



This equation shows that if the solid is of such construction that its 
moment of inertia can be increased or decreased at will, its angular 
velocity will immediately slow down or speed up so as to keep the 
product Co) constant. This is an example of the principle of con- 
servation of angular momentum; for, by (14' 6), Co) is the angular 
momentum in this case. Acrobats and other actors frequently utilize 
this principle to increase or decrease the spin of their bodies, which 
they can do by drawing in or spreading their arms or legs. 
The vector form of equations (16-4) is 

^4-toxH = M, (16-7) 

as may be verified by replacing H and to by their equivalent forms 
(16-2) and (164), respectively. In (16-7) M is the moment about 
of the external forces acting on the body. 

It is to be noted that (16-7) is the differential equation of motion 
referred to the axes of the moving trihedron. 

If H is constant, then dH/dt^O and (167) becomes 

toxH = M. (16-8) 

The vector product to x H is sometimes called the gyroscopic 
couple. It is also called the centrifugal couple. This vector is 
perpendicular to both to and H and is therefore perpendicular to the 
plane of those vectors. Hence in this case M is perpendicular to 
the plane of to and H. 
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If the moment of the external applied forces is zero, then (16-7) 
becomes jra 

-z-- = Hxc*>. (16-9) 

Uv 

Finally, if co x H = 0, which would be the case if <o and H should 
coincide, then (16-7) becomes 

f=M. (16-10) 

Attention is called to the dual utility of equations (16-4), (16-5) 
and (16-6) in the solution of gyroscopic problems: 

(1) If the external forces acting on the gyroscope are known, the 
moments M X9 M V9 M z can be found. Then the right-hand members of 
the equations are known, and the motion of the gyroscope can be 
determined by integrating the resulting differential equations. 

(2) If the motion of the gyroscope is given, the left-hand members 
of the equations will be known and then the moments M X9 M y > M z 
which could produce the given motion can be determined by 
substituting in the left-hand members the assumed or given data 
concerning the motion. 

These two uses will be shown in the next chapter. 

17. Kinetic energy of the body 

If a rigid body is rotating with angular velocity w about an 
instantaneous axis passing through a fixed point 0, and r is a 
position vector drawn from to a particle of mass m whose co- 
ordinates are (x, y, z), the linear velocity of the particle is given by 
the formula v = co xr = ^i + ^j+^k by (13-1), 

and the scalar components v x9 v y , v z are given by (13*2). The kinetic 
energy of the body is given by the formula 

T = pmv 2 = Pm(v| + v*+ v*) . 
Replacing v 2 , v* 9 v* by their values as given in (13*2) we obtain 

T = pm[(w 2 z-ce>3y) 2 4-(6> 3 ~^^ 
or 



- 2<y x 6> 3 IZmzx 2^ 2 ( 

Now if the principal axes of inertia are taken as coordinate axes, 
the products of inertia disappear and we then have, by (14-4), 

. (17-1) 
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From (16-1), (16-2) and (14-6) we have 

co-H = ^<4 + JRtfJ + <7(yf = 22* by (17-1). (17-2) 

Hence the kinetic energy is greatest when H is parallel to <o. This is 
the same case as o> x H = 0. 

Since ^'H = a)Hcos0, we have, from (17*2), 

Then since o), H and T are all positive quantities, 6 must be an acute 
angle, either positive or negative. 
From (17-1) and (14-6) we obtain the further relations: 






(17-4) 



If the moments of the external forces are zero, then Euler's 
equations (16-5) become 



(17-5) 



+ (B- A) 0)^^ = 0. 



On multiplying these equations by o) l9 &> 2 and o> 3 , respectively, and 
adding the results, we get 

d(jt) 9 . da)* 



Integrating this, we have 



= a 



2r by (17-1). 

Hence the kinetic energy is constant when the moments of the 
external forces are zero. A rigid body rotating about an instanta- 
neous axis through its center of gravity and acted on by no other 
forces than its own weight, would come under this case. 
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18. Eider's angles 

The orientation of a rigid body moving about a fixed point is best 
fixed by three angles 0, $ and ^, called Euler's angles. These con- 
nect the moving axes to the fixed axes and are shown in Fig. 17. 




Fig. 17 



The moving trihedron XYZ can be brought from the fixed 
position X^Z^to the position shown in Fig. 1 7 by the following 
procedure, carried out in the order stated: 

(1) Botate through angle t/r about OZ l} thereby bringing OY l 
into coincidence with ON. 

(2) Eotate through angle 6 about ON, thus bringing OZ to the 
position shown in Fig. 17. 

(3) Eotate through angle <f> about OZ, thus bringing OX to the 
position shown in the figure. 
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19. Components f the angular velocity in terms of Eider's angles 

To find the values of <y l9 fc> 2 , &> 3 in terms of 0, 0, ^ and their 
time derivatives, we consider the projections of all the angular- 
velocity vectors on each of the moving axes. In Fig. 18 the plane 




Fig. 18 



OZ^Z is perpendicular to ON, and BQ is the intersection of the 
plane OZ^Z and the XT-plane. This line BQ is to be used as 
a bridge for transferring vectors from the O-Z^Z-plane to the 
-XT-plane. 

Inspection of Fig. 18 shows that the projection of ^ on BQ is 
i]f sin 6, and this projection therefore has components ^ sin 6 cos <$> 
along OX and ^sintfsin^ along 07. The projections of 6 on 0-3T 
and OF are $sin^ and 0cos0, respectively. If co denotes the 
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angular velocity of the moving trihedron O XYZ about an 
instantaneous axis through 0, the components &> X9 a) y , &> s along the 
moving axes are therefore 

co^. = 6 sin <p ijf sin Q cos $ 9 \ 

*& y = 6 cos <j) + i/r sin 6 sin S , i (19*1) 



These equations are sometimes called Euler's Mnematical 
equations. 

The kinetic energy of the body in terms of 0, <f>, i/r and their time 
derivatives is, by (17-1), 



T = %[A (6 sin - $ sin 6 cos 0) 2 

+ J3 (<9 cos ^ + ^ sin sin 0) 2 + C(<f> + fi cos 5) 2 ]. (19-2) 



For a solid of revolution about the 2-axis, we have =A and 
then (19-2) reduces to 

(19-3) 



20. Rotation of a rigid body about a fixed axis 

The fundamental vector equation (15-1) and its Cartesian 
equivalents (15-2) hold for the motion of a rigid body about any 
axis passing through a fixed point. They therefore hold equally 
well for a, fixed axis passing through the fixed point. 

Let us take the #-axis as the fixed axis. Then from (14-3) we have 

H x = a) x *Iim(y 2 + z 2 ) o) y ^mxy o)^mxz. 

But since the rotation is entirely about the #-axis, o) y and o) z are 
zero. Hence in this case 



or, for any line, H=Io). 

Then, by (15-1) or (15-2), 



3-2 
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If 6 is the angle through which the body has rotated from a given 
position, then ft> = d0/<K and (20-1) becomes 

ig-Jf, (20-2) 

or Id = M. (20-3) 

This equation, in any of the three forms just given, is the funda- 
mental equation for the rotation of a rigid body about a fixed axis. 
It corresponds to the fundamental equation of Art. 10 for rectilinear 
motion. 



CHAPTER III 
Theory of the Gyroscope 

21. Definition and behavior of a gyroscope 

A gyroscope may be broadly defined as a solid body capable 
of rotating at high angular velocity about an instantaneous 
axis which always passes through a fixed point. The fixed point 
may be the center of gravity of the solid or it may be any other 
point. 

More specifically, however, the usual form of gyroscope is a 
mechanical device the essential part of which is a flywheel having 
a heavy rim and so mounted that, while spinning at high speed, its 
axis of rotation can turn in any direction about a fixed point on that 
axis. The flywheel can be mounted so as to turn about its center of 
gravity by means of double gimbals called a Cardan suspension 
(Fig. 19). ^ ^ 

The familiar toy top and 'spool dancer' are forms of gyroscopes 
in which the fixed point is the point of contact of the top with the 
floor or table on which it spins. 

A rapidly spinning flywheel, mounted as described above, 
resists any effort to change the direction of its axis of spin. In other 
words, the flywheel persists in maintaining its plane of rotation. 
For example, if the axis of spin is horizontal and we apply a steady 
downward force to the inner ring at one end of the spin axle, the 
end of the axle is not carried downward as we should expect but 
begins to skew around in a horizontal plane. This skewing around 
of the spin axis is called precession.} It arises in obedience to the 
fundamental relation dH/(B = M. Since the applied external 
moment cannot change the rate of spin of the flywheel, it increases 
the angular momentum of the mounted wheel and ring by making 
it rotate about a vertical axis. 

4- "The direction of this precession about a vertical axis depends on 
the direction of spin of the flywheel and also on which end of the 
axle is pressed downward. When the wheel is spinning in a given 
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direction, a downward pressure on one end of the axle will cause the 
axle to process in one direction, whereas a downward pressure on 
the other end of the axle will cause precession in the opposite 
direction. 




Fig. 19 

' f * 

f\ N i 

22. Steady or regular precession <* ^ 

Let us now assume that a heavy flywheel ^s mounted so that its 
axis of spin always passes through the center of gravity of the gyro 
as indicated in Fig. f|. Let us further assume that the 
gyro is spinning with constant angular velocity <j> about the axis 
OZ, that OZ makes, a constant angle 6 with the vertical axis OZ V 
and that the axis of spin precesses about OZ 1 with constant angular 
velocity X Constant precession under the conditions just described 
is called steady precession. It is required to find what external 
moment is necessary to maintain such steady precession after once 
started. 

As the spin axis precesses around OZ with constant angular 
velocity A, the plane ZOZ 1 (mlle^^ 



rotates about OZ 1 with angular velocity 1 
Hence if we take the rectangular trihedron - J YZ as the system 

* The mounted flywheel in a gyroscope is frequently called a gyro and 
will often be called that in this book. 
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of moving axes, the angular velocity of the trihedron about is 
also Q.* Hence the components of this angular velocity along the 
moving axes are 

iX = Osin0, 




Fig. 20 

If denotes the moment of inertia of the gyro about its axis of 
spin and A denotes the moment of inertia about the transverse 
axis OX, then, since OZ and OX are principal axes of inertia for 
the center of gravity and since the total angular velocity vector 
along OZ is D cos 6 + <fi, we have 



* Note that the axes of the trihedron are not fixed in the gyro. 
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Now since 9, <j> and Q are all constants, we get 



w J 7 JL v 3 JT , v 

a (X5 f 

Substituting the above values into (16-4), we get 



(22-1) 



The required moment is thus entirely about the y-axis, or in the 
azimuthal plane ZOZ^ 

The expression for M y can be rearranged and written in the 
following equivalent form: 



Jf = 

<p J 



~. (22-2) 

As G > A for a flywheel of the usual type,* the quantity in 
brackets in (22*2) is positive. Hence M y is a positive quantity and 
as a vector it would be drawn from along the positive y-axis. 
Then, by the right-handed-screw rule, this moment tends to 
increase the angle 6. Hence for steady precession with a constant 6 
there must be a turning moment in the azimuthal plane, and this 
moment must tend to increase 6. 

When a gyroscope is spinning at high speed, the precessional 
velocity D is negligible in comparison with the spin velocity ^5. 
Hence in that case we may neglect the fraction O/$5. Then (22-2) 

becomes ' > 

(22-3) 



* For a flat disk, C = |mr 2 and A = Jwr 2 ; and for a wheel in which most of 
the mass is in the rim, O = mr 2 , approximately, and J. = Jmr 2 , approxi- 
mately. Hence, for the usual rotor of a gyroscope, C>A. 
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If (9 = 90, this becomes 



#). (22*4) 

This special equation (22-4) is of great practical importance. 
When written in the form 

Q = ^ 3 (22-4) 

it states that the velocity of precession varies directly with the 
external moment and inversely as the velocity of spin of the gyro. 
Thus, for a given external moment the greater the spin velocity, 
the slower the precession. 
Another aspect of (22-4) will be discussed in Art. 24. 




Fig. 21 

23. Direction of precession 

In practical applications of the gyroscope, it is of the utmost 
importance to know the direction in which the spin axis will 
precess under given conditions. In the highly important case 
where 6 = 90 and M = CCltf*, reference to Fig. %ill show that the 
spin axis OZ will lie along OX, Then, according to the right-handed- 
screw rule, a precession to the right about OZ^ will swing OX 
around toward OF; that is, it will swing the spin vector around 
toward the moment vector, as indicated more clearly in Fig. ML 
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Hence we can state the following rule for determining the direction 
of precession: 

An applied external moment causes the positive end of the spin axis 
toprecess toward the positive end of the moment axis, the positive end 
of each axis being that determined by the right-handed-screw 
rule. 

An extension of this rule is given at the end of Art. 33. 

24. Free precession 

The formulas of Art. 22 give the external moment required 
to maintain a steady precession after it 
is once started. It is conceivable that 
for some velocity of spin the preces- 
sion might continue after the external 
moment is removed. To find such a 
velocity we put M y = in (22-1). We 
then get 



or 



Equation (24-1) expresses the condition 
for free precession. When this condition is lg " 

satisfied, a steady precession would continue indefinitely without 
the action of an external moment if no frictional forces interfered. 

The negative sign before the right-hand member of (24-1) should 
be noted. This sign shows that when C>A and #<90, O and <p 
have opposite signs. This means that the precession would be 
backward or retrograde; whereas if C< A, Q, and <j> have the same 
sign and precession would be direct. 

Equation (24*1) shows another interesting result. If\C A\<O, 
the fraction Cj(C-A)Gosd is greater than L Hence in that case 
the precession .Q is greater than the spin <j>. Furthermore, if 6 is 
near 90, O is very large. 

At this point we wish to call attention to the fact that when a 
gyroscope is processing, the axis of instantaneous rotation does 
not coincide with the axis of $pin. The total velocity about the 
instantaneous axis is the geometric or vector sum of the spin 
velocity and the precessional velocity, as indicated in Fig. 22. 
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25. Unsteady precession* General motion of a gyroscope 

In Art. 22 we discussed the motion of a gyroscope under the 
conditions that the spin velocity, the precessional velocity and the 
angle 8 were all constants. We now drop these restrictions and 




Fig. 23 

assume that the angle 6, the spin velocity (j) and the precessional 
velocity O may all vary in any manner. We are to find the necessary 
external moments in this case. 

Referring to Fig. 23, let the axis of spin lie along OZ, and let this 
axis precess about OZ with variable angular velocity ^. The time 
rate of change of 6 will be denoted by 6 and is seen to be a vector 
along the positive direction OY. The spin velocity with respect to 
the azimuthal plane ZOZ is denoted by $ and is laid off as a vector 
along OZ. 

The position of the moving trihedron XYZ is given at any 
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instant by the angles i/r and 6 (Fig. 23). The projections of the 
angular velocity of the trihedron on the moving axes are i/rsind 
on the x-axis, 9 on the y-axis, and i/r cos 6 on the si-axis, as indicated 
in the figure. Hence we have 

(25-1) 



Because of the spin velocity along OZ, the total angular velocity 
along that axis is (j> + ft cos 6. 

Now if the axes OX, 07, OZ are principal axes of inertia of the 
gyro for its center of gravity 0, the components of the angular 
momentum of the gyro along the axes are 



Consequently we have 



On substituting into (16*4) the components of the angular 
velocity and angular momentum, and the values of the derivatives 
given above,, we get 



Ad-Aijr* sin 6 cos + C($ + ft 



These equations can be simplified somewhat and written in the 

forms 

Aty sin v ZAVi/r cos V + C'($> + ty cos V ) V = M x , 

(25-2) 



These are the fundamental differential equations of motion of a 
symmetrical gyroscope. It will be seen that they reduce to (22-1) 
when 6, <p and t/r are constants. 
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It will further be seen on comparing the equation for My with equa- 
tion (22-1) that My in the general case equals the moment required 
for maintaining steady precession plus the accelerating moment A 8. 

When a gyro is spinning at high speed, the quantity <j> is the 
dominant term in the left-hand members of equations (25-2), and 
all other terms are negligible in comparison with it. Hence if we 
neglect all terms not containing ^, those equations become 



(25-3) 



And if the spin velocity is constant, or nearly so, we may neglect 
. Then equations (25-3) become 

=M x> } 

\ (25-4) 

My.) 

In some applications of the gyroscope there is no external 
moment about OZ\ in those cases M z = 0. Then when the axes of the 
moving trihedron coincide with the principal axes of the gyro, we 

have H x =Aa> x , 

and 

On substituting these quantities into (16-4), we get 
A(b 






dH x .. dEL . . dH, _ d 
--^ -=**> -- 



A(b y + Aa) x a) z - 



From the third of these equations we get 

(j) + (j) z = a constant = O, say. 

Now substituting this value of <j) + o) z into the first two equations, 
we have Afi + CQ-A-M^ 



Here it is to be noted that O represents the total angular velocity 
about OZ, the spin velocity plus the ^-component of the angular 
velocity of the moving trihedron. 



46 THE GYROSCOPE 

26. Steady precession of an nnsymmetrical gyroscope 

In the preceding articles it has been tacitly assumed that the 
gyroscope was symmetrical, that is, that B=A. This, however, is 
not always the case. In some gyroscopes of practical importance 
B^A. We now proceed to find the equations for steady or regular 
precession in this case. 

Since for steady precession the quantities 6, <j> and i/r ( = O) must 
all be constant, equations (19-1) become 

a) x = iff sin 6 cos <f>, 



d) z = (j)~{-^f cos 6. 
Hence 

do)~ ;;./) j 

= (pyf sm u sin 0, 

^r, , 

= <^sin#cos0, 



On substituting into (16*5) the above values of o) x) o) y , o) z and their 
derivatives, we get 



(A - B) ijr* sin 2 d sin ^ cos $5 = 

(26-1) 

The above equations (26-1) are referred to the axes of a trihedron 
which is fixed in the body and therefore spins with the same 
angular velocity as the body. Since it is rather absurd to use 
coordinate axes which may rotate thousands of times per minute, 
we choose a set of axes which rotate at a speed comparable to the 
speed of precession of the gyroscope. We therefore take the axis of 
spin, the axis of nodes, and a transverse axis perpendicular to each 
of these as the axes of the moving trihedron. The moving axes will 
therefore be OZ, ON and OB (Pigs. 17 and 18). Then since angle 
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XON = 90 - 0, the moments, about ON and OB in terms of M 9 
J f . are , _ ,,. 



M OB = M x cos <j> My sin 0.] 
Substituting into (26-2) the values ofM x and M y from (26-1), we get 



H- (B J.) (cos 2 sin 2 0) 00" sin 6 

-f (7(sin 2 + cos 2 ) 0^ sin Q = 
(^ B) ^ 2 sin ^ cos 6 sin i cos <p -f 2(-4 5) $ijr sin ^ sin 5 cos <j> 

=J1 OJ?> . 

(J[ - S) ^ 2 sin 2 6> sin <p cos ^ = lf g . 

Now changing to double angles in ^ by means of the relations 
2 sin 5 cos <p = sin 20, cos 2 (j) sin 2 ^ = cos 20, 
sin 2 = |(1 - cos 20), cos 2 J = |(1 + cos 
and then rearranging the results slightly, we get 



(26-3) 



J 

A ) ($ cos 6 + 20) T/r sin 6 cos 20 = M ON , 

- B) (i/r cos 6 + 20) ^ sin 6 sin 20 = Jf o ^ , 

\(A - B) fi* sin 2 6 sin 20 = M z . 



These equations (26-3) show that the moments about all the axes 
depend on the angle 0, that the moments are pulsating with period 
TT, and that the pulsation frequency is therefore twice the rotational 
spin of the gyro. The moments about OB and OZ are zero twice 
during each rotation of the gyro, but the moment about O^V" is 
never zero unless 



or 

which would rarely or never occur. 
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On rewriting the first two equations of (26*3) In the forms 






-f i(.B - A ) | T cos 6 + 2 1 (j>\jf sin 6 cos 2$ = M ON , 
%(A - B) |T cos 6 + 2 J <f>\jr sin 6> sin 2^ = M OB , 



we easily see that when the spin velocity <f> is very large and the 
precessional velocity ^r is consequently very small, the fraction 
ijrl<j> is negligible. Furthermore, in this case the quantity 



is also very small and may be neglected. Hence, on dropping the 
negligible quantities just mentioned, we get 

G(j)i}f sin 6 + (B A ) <f) $ sin 6 cos 2<f> = M ON , } 



for the moments in a fast-spinning unsymmetrical gyroscope. 

If we put EA in equations (26-3) and (26-4), we get equations 
(22-1) and (22-3), as we should. 

27. Gyroscopic resistance 

A rapidly spinning gyro is a stubborn affair. It resists any effort 
to turn it aside from its plane of rotation. To state it otherwise, it 
resists any force that tends to change the direction of its axis of spin. 

In Art. 22 we found that an external moment which tended to 
increase the angle 6 had to be applied in order to make the gyro 
precess at a constant rate. By Newton's third law of motion the 
gyro resisted the precession by an equal and opposite moment. 
This opposing moment is known as the gyroscopic reaction moment, 
the force of which is a reaction against the bearings and supports of 
the gyroscope. It is an inertial reaction ('gyroscopic inertia') 
analogous to centrifugal force and is the property that gives gyro- 
scopes their practical usefulness. Henceforth in this book the term 
'Gyroscopic Reaction Moment' will at times be abbreviated G-.R.M. 

As an illustration of this gyroscopic reaction moment, let us 
suppose that a gyroscope is rigidly installed on some moving body, 
such as a vehicle or a ship, and that its spin axis is horizontal. Then 
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if the moving body changes its direction of motion while the gyro 
is spinning at high speed, such a change in direction will change the 
direction of the spin axis by the same amount. This enforced 
change in direction of the spin axis is a forced precession about a 
vertical axis, and this precession will in turn automatically call 
forth a reaction moment against the bearings and supports of the 
gyroscope. This reaction moment is equal and opposite to the 
external moment that would be required to produce the forced 
precession. This instant gyroscopic action on the supports of the 
gyroscope tends to prevent the moving body from changing its 
direction of motion and thereby helps to steer it straight ahead. 

Since the gyroscopic reaction moment is equal and opposite to 
the external applied moment, the equations giving the reaction 
moments can be obtained from the corresponding moment equa- 
tions of Arts. 22, 25 and 26. Denoting gyroscopic reaction moment 
by jfir, we have, by definition, 

K=-M, K X =-M X , K y =-M y , etc. 
Hence from equations (22-l)-(22-4) we get 

(27-1) 



= |" 



. (27-2) 



c r + (C-A)cosd (^ xfi),. 



&), (27-3) 

(27-4) 
Likewise, from equations (25-2)-(25-4) we get 



(27-6) 
} 
K X =-C6<J>, i 
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If the quantities in the right-hand members of the above 
equations (27*I)-(27*7) are known, the reaction moments can be 
found. 

The gyroscopic reaction moments for an unsymmetrical gyro- 
scope are found from equations (26-3) and (26-4) to be 



(27-8) 
K OB = - [\(A - JS) (ft cos 6 -h 2<j>) ft sin 6 sin 2$ , 



^* --[ 
and 

jg: oiv = - [0^ sin 6 + (B - A) ^ sin d cos 2U 

V (27-9) 

n 2rf] . J 



We see at once from these equations that the G.R.M/S in an 
unsymmetrical gyroscope are pulsating with a frequency twice the 
rotational speed of the gyro. If these pulsations happen to be nearly 
the same as the natural vibration frequency of some part or parts 
of the machine to which the gyroscope is attached, resonance will 
take place and may cause serious damage. 

Because of the very great importance of the principle of gyro- 
scopic reaction, we state again by way of emphasis that if a gyro 
is forced by any means to precess with an angular velocity O about 
any axis, it will immediately exert on its bearings or supports a 
reactionary moment equal and opposite to the external moment 
that would be required to produce the precession. The magnitude 
of this gyroscopic reaction moment is given by the formula 

K^C&t (27-10) 

when the axis of precession is at right angles to the axis of spin. 

The fundamental relations between external moment and pre- 
cession, and between precession and gyroscopic reaction moment, 
have been succinctly stated by Gramme!.* 

To a moment M the gyro responds by a precession SI; a precession 
Si induces in the gyro a reactionary moment K. 

From the standpoint of practical usefulness equation (27*2) and 
its special case (27-4) are the most important formulas in gyroscopic 

* Der Kreisel, p. 72, 1920 edition. 
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theory. Particular attention is called to (27-4). Reference to Fig. 20 
will show that a moment about Y tends to change the direction 
of the axis of spin and therefore to deflect the gyro from its plane 
of rotation. When a heavy gyro is spinning at high speed, the 
right-hand member of (27-4) is large and therefore the reaction 
moment K is also large. This means that the gyro offers great 
resistance to being deflected from its plane of rotation. The higher 
the spin velocity and the greater the moment of inertia of the gyro, 
the greater is the deviation resistance. 

A glance at equation (27*4), or (27-10), shows that there can be no 
gyroscopic resistance unless there is precession, for when O is zero, 
K is also zero. 

28. Direction of the gyroscopic reaction moment 

We come now to a consideration of the all-important question of 
the direction of the gyroscopic reaction moment, for in the applica- 
tions of the gyroscope it is even more important to know the 
direction of the reaction than to know its magnitude. The vector 
equation K = C(<f> x SI) shows, according to the definition of the 
vector product, that the reaction K tends to turn the spin vector 
into the direction of the precession vector. This fact is also indicated 
in Fig. 21. Hence we can state the following rule for determining 
the direction of the gyroscopic reaction moment: 

The gyroscopic reaction moment (G.R.M.) tends to turn the positive 
end of the spin axis toward the positive end of the precession axis, the 
positive end being determined in each case by the right-handed- 
screw rule. 

Note that the dotted curved arrows in Fig. 21 are drawn from one 
known vector to another known vector in each case. 

29. Motion of a free gyroscope 

A free gyroscope is one on which no external forces act, or one 
suspended at its center of gravity and acted on by no other force 
than its own weight. In either case the external moment about 
its point of suspension is zero. Then, since M = 0, equation (151) 
becomes 



Consequently, H = a constant. 

4-2 
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Since equation (1 5-1) refers the motion to a set of 'fixed axes, and 
since H is constant in both magnitude and direction, it is evident 
that H has a fixed direction in space. The fixed line on which H 
lies is called the invariable line, and any plane perpendicular to 
it is called an invariable plane.* 




Fig. 24 



To find the motion of a free gyroscope we therefore take the 
invariable line as the fixed vertical axis. The #-axis is taken in the 
, and the t/-axis is taken perpendicular to this plane, as 



* The solar system is not acted on appreciably by any external forces. 
Hence its moment of momentum is constant. It therefore has an invariable 
line. The plane through the center of the sun and perpendicular to the 
invariable line is the invariable plane of the system. It is inclined to the 
plane of the ecliptic by about 2. 
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indicated in Pig. 24. Then since the external moments are all zero, 
equations (25-2) become 






Our problem now is to integrate these equations. 

From Fig. 24 we have 



Now integrating the third of equations (29-1), we get 

C(<j> -f ijr cos 6) = a constant. 
But H 3 = H cos = G((j> + i/r cos 6) = a constant. 

Therefore, since H and H cos $ are both constant, cos 6 must also 

be constant. Hence n , , 

6 = a, constant, 

and = 0, = 0. 

Substituting these values for 6 and 8 into the first two of equations 
(29-1), we get , 

(29-2) 



From the first of these equations we get 

^ = a constant = O, say. 
Substituting this for ^ into the second of equations (29-2), we get 



C A 

Hence <j> = -- -^j O cos 5, a constant. (29-3) 



We have now shown that 6, i/r and are all constant for a free 
gyroscope. The motion is therefore that of steady precession. It 
will be noted that equation (29*3) is the same as (24-1), and the 
precession is therefore a free precession. 
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We now show that in this motion the instantaneous axis of 
rotation makes a constant angle with the invariable line. Since the 
rotation about the instantaneous axis is the geometric sum of the 
spin velocity and the precession velocity, the instantaneous axis is 

H 



H 




Fig. 25 Fig. 26 

represented by the vector to in Fig. 25. If a denotes the angle 
which GJ makes with the invariable line we have, by the law of sines, 



sna 



from which 



(29-4) 



Since $5, 6 and o) are all constants, sina is constant and therefore 

a = a constant. 

The instantaneous axis thus describes a cone of semivertical angle a 
about the invariable line as an axis, as shown in Eig. 26. Equation 
(29-4) shows that a depends on the velocity of spin, the velocity of 
precession, and the angle 6. 
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30. Stability of the morion of a free gyroscope 

We now investigate the stability of the motion of a gyroscope 
on which no external couples are acting. Consider a gyroscope that 
is free to move in any manner about its center of gravity as a fixed 
point. To simplify the problem we assume that the axis of spin is 
horizontal, that the gyro is spinning with angular velocity 5, and 
that the spin axis may or may not be processing with angular 
velocity i/r. Assume now that the axle of the gyro is struck a sharp 
impulsive blow I, directed vertically downward as indicated in 
Fig. 27. Our problem is to determine the effect of this blow on the 
subsequent motion of the gyroscope. 

The immediate effect of the blow is to give the spin axis an angular 
velocity 6 Q in the azimuthal plane. Hence the initial conditions of 
the problem are 



On putting = J?r and M x = M y = M z = in equations (25-2), these 

equations become 
H 



(30-1) 
C0-O. 

From the third of equations (30-1) we get 

C<j) = constant = Co), say. 
Then ^ = 6> = ^ and $4 = &tf-f0 . 

Substituting this value of <j> into the first two of equations (30-1), 



Integration of the first of these equations gives 

-A^ + Co)8 = G 1 . 
To find C v put ijr=i/r and 6 = \ir. Then 
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The above equation now becomes 



or 



(30-3) 




Fig. 27 



Substituting this value of ^ into the second of equations (30-2), we 
obtain 



or 



Solving (30*4) by the method of undetermined coefficients, we 
have 

r* + 



whence 



r = 



THBOBY 57 

Then the complementary function is 



Since the right-hand member of (30*4) is a constant, we find a 

particular integral by putting 

6 $ = P, a constant. 
Then = 0, 0=0. 

Substituting these values of 6 and 8 into (30-4), we get 
fC<oV IC^\^ _G^, 

\A) (A) 2 A r ' 

whence P= |_^ . 

The complete solution of (30-4) is therefore 



mi A n Q (@\ n Q 

Then u = (7 j- cos \ -r- 1 1 Cg r- sin 

jfji. \ jOL / -cl 

To find Cg and (7 3 we utilize the initial values 6 = JTF, <9= ^ 0? when 
t = 0. Then we find C 2 = Ad /Ca) ) C 3 =Ai/r Q /Ca). The value of 6 is 
therefore 



a" OiTo 

(30-5) 



2 



where /? = tan" 1 
Now substituting this value of into (3G-3), we get 



(30-6) 
v-o- / 

Integration of this gives 
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To find (7 4 we utilize the initial condition i/f = when t = 0. This gives 



Then becomes 



The value of 6 in (30-5) can be written in the equivalent form 

(30-8, 



which shows that the periodic parts of 6 and i/r have the same 
period but are out of phase by a quarter period. 
The final solution of equations (30-1) is therefore: 



(30-9) 



These equations show that after the impulsive blow the gyro 
continued to spin with the same angular velocity a& before, but 
that the axis of spin was displaced from the position 6 = f TT, i/r = to 
the mean position d = ^7TAt/rQ/Ca), T/r=A0 Q ICo), and that it there- 
after oscillated about this mean position with simple harmonic 
motion both vertically and horizontally, the vertical and horizontal 
oscillations being out of phase by a quarter of a period. It is to be 
noted here that if the gyro axis was precessing at the instant of the 
blow, the progressive precession ceased after the axis reached the 
position for which T/r=A6Q/C(0', for equation (30-6) shows that the 
only precession is periodic. We note further that the displaced 
position of the axis depends on the intensity of the blow (# ) and on 
the angular velocity of spin the greater the velocity of spin, the 
smaller the displacement. 

The amplitude of the harmonic oscillations is (AIC(t))*J(6 2 + i/r*) 
and their frequency is Ca)/2nA, thus showing that the higher 
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the spin velocity, the smaller the amplitude and the higher the 
frequency. 

The preceding investigation shows that the motion of a free 
symmetrical gyroscope is stable. Moreover, if A and C be inter- 
changed, the motion is still stable. Hence the motion is stable 
whether the axis of spin is about the A -principal axis of inertia 
or about the C f -principal axis. 

To investigate the stability of a free unsymmetrical gyroscope, 
we utilize equations (16-5) and (19-1). We first consider the case in 
which the spin axis coincides with OZ and assume that the axle is 
struck a sharp blow in the ^0^-plane and at right angles to the 
axle. The immediate effect of the blow is to impart an angular 
velocity $ to the axle (see Fig. 20). 

Turning now to equations (19-1), we notice that a) x and c^ are 
very small in comparison with o) sy because 6 and i/r are very small 
in comparison with the high spin velocity (j>. Hence we shall 
neglect the product (*) x (*) y in the third equation of (16*5). Then, 
since M x M y = M z in a free gyroscope, equations (16-5) become 



(3040) 



From the third of these equations we get 
o) z = a constant = fl, say. 

Substituting this value of o) g into the first and second of equations 
(30-10), we have , 



Now differentiating the first equation in (30-11) and then substi- 
tuting the resulting value of da) y ldt into the second equation, we 
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Here the auxiliary equation is 

ABr* + (A - O) (B - C) O 2 = 0, 

from which r _ + < 

The period of oscillation resulting from the blow is therefore 



AB 



2n 

J- "~7T 



This period is real either for A>B>G or for C>B>A, and the 
gyroscope is therefore stable whether the axis of spin is along the 
A -principal axis or the O-principal axis of the gyro. 

We next consider the case in which the axis of spin coincides with 
the y-axis and therefore the ^-principal axis of the gyro. In this 
case o) x and o) a are very small and their product negligible, while <o y 
is large. Hence equations (16-5) become 






(30-12) 



From the second equation of (30-12), we get 
(i) y =a constant Q, say. 

On substituting this into the first and third of equations (30-12), 
we get 

A + (C-)^ a = 0,\ 

(30-13) 



Differentiating the first equation and substituting the resulting 
value of do) z /dt into the second, we find 
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Here the auxiliary equation is 

-B)(B-C) = 0, 
I(A-B)(B-O 



from which r = + 1 



AC 



The values off are real, since the factors under the radical are both 
positive for A>B>C or G>B> A. Hence w x is not periodic and 
the spin axis does not oscillate after the blow. The gyroscope is 
therefore unstable when spinning about the intermediate or 
J5-principal axis of inertia. 

We may therefore state that a free unsymmetrical gyroscope is 
stable when spinning about either its greatest or least principal 
axis of inertia, but is unstable when spinning about its intermediate 
principal axis. 

31. Motion of the spin axis in unsteady precession 

We next consider the motion of a gyroscope which is spinning at 
high speed and also precessing under the action of a constant 
external moment. Let N denote the external moment and let it be 
applied about the axis of nodes (y-axis). Then 

My = N and M X =M S = Q., 
With 6 = 90, as in Art. 30, equations (25-2) now become 

(31-1) 




These equations are to be solved subject to the initial conditions: 

. = J*r, = 0, fi = fi Q , when = 0. 
From the third equation of (31-1) we get 

(7$5 = a constant = <7o> 5 say. 
Then ^ = w and <fi = (i)t + <f> Q . 

Substituting this value of <j) into the first two of equations (3 1 1 ) , we 

have - > = 

(31-2) 
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Integration of the first of these equations gives 

^Ai/f + Ca)0 = C l . (31-3) 

To find Oj we put 6 = JTT and i/r = ^ . Then 



and (31-3) becomes 



from which ^ = ~ (0 - 1 j + ^ - (31-4) 

When this value of ^ is substituted into the second of equations 
(31*2), we have 



which becomes 



Solving this equation by the method of undetermined coefficients, 

as in Art. 30, we have 

a /I - Oft), OOK 
__ fj g|j^ _i_ (j QOS t 

To find a particular integral we put 9 p = P and substitute it and its 
second derivative into (31-5). We thus get 



/^\* # /<?6A*7T 

(A) ^-A + \A) 2 



A 



The value of # is therefore 

(7w, JVJ. TT 



m-i A f^ \J\U \JVJ ~ \JW . \SW 

A A ^ A A 

jtJL ud jOL jd 
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To find C 2 and C s we put $ = f 71, $ = when t = 0. Thus 



TT ^ NA n A 



whence 



NA 



And 



Then 6 and 6 become 



or = 0. 



i i Ctojt 
--^nlsm-jtf. 



(31-6) 



Now substituting this value of 6 into equation (31-4), we get 



Integration of this gives 



Putting ^=0, = 0, we get (7 4 = 0. Hence 

N A / , N\ . Co) 

The final solutions of equations (31-1) are therefore 
(ft = fa = a constant or <f> = o)t, 
n TT A [ , N\ A f , N\ Cut 



(31-7) 



N At, N 



(31-8) 



These equations show that the velocity of spin remains constant, 
that the precession continues at a uniform rate except for a periodic 
disturbance which alternately speeds it up and slows it down, and 
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that the axis of spin oscillates in the azimuthal plane. The mean 
position of the spin axis is not in the plane 6 = \TT but lies in the 
conical surface A , , T 

-IT A (N 



If I/TQ = N/C(JI), or N CtoifrQ, all vibration of the spin axis dis- 
appears and it precesses at a constant rate in the plane B^\TT. 
Since the equation N = Ca>ifr is the same as (22-4), the motion is 
steady precession. There is therefore no vibration of the axis when 
the precession is steady. 

On writing the equation for 8 in the form 

. TT A , N A , N . Cut 



and comparing this result with the value for T/T, we note that the 
oscillations have the same periods and amplitudes in both cases but 
that they differ in phase by a quarter of a period. 

Note. This book deals with the ideal gyroscope, free from friction 
and all imperfections of construction. Although the spin velocity of 
the rotor is kept constant by an electric motor, the gimbal bearings 
are not absolutely perfect, and it is not possible to keep the 
center of gravity of the rotor and gimbal assembly at the exact 
point of intersection of the gimbal axes. These departures from 
the ideal state cause the spin axis to depart slightly and slowly from 
its initial orientation. This slow change in orientation of the spin 
axis is called drift. Designers and manufacturers of gyroscopes 
make every effort to reduce drift to a minimum. 



CHAPTER IV 

Theory of the Gyroscope (Continued) 

Motion of a Gyroscope under the 
Action of Gravity. The Top 

32. The differential equations of motion 

Many gyroscopes are suspended in such a manner that the point 
of suspension is not at the center of gravity. In such cases the weight 
of the gyroscope has a moment about the point of suspension. The 
ordinary top is a gyroscope of this type. Fig. 28 represents a 
gyroscope suspended in the manner just mentioned. 

To study the motion of a gyroscope under the action of gravity, 
we take a set of moving axes X TZ with the 2-axis along the axis 
of spin, the #-axis in the ZOZ^ plane, and the g/-axis perpendicular 
to this plane, as indicated in Fig. 28. Let W denote the weight of 
the gyroscope, or top, and let I denote the distance OG. Then since 
the 2/-axis is perpendicular to the plane ZOZ ly the moment of the 
weight is about the 2/-axis. Hence we have M x = 0, M M Wl sin 6, 
M 3 = 0. Substituting these values into equations (25-2), we get 

cos 0)0 = 0, 



as the differential equations of motion. Before attempting to find 
the general solutions of these equations, we shall consider particular 
cases and aspects of them. 

33. Forced steady precession 

For steady precession we must have (Art. 22) = constant, 
96 = constant and ^ = constant. Then $ = $ = 0, ^r=0 and 

(j) + ^f cos 6 = constant = &>, say. 
5 65 SG 
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When these values are substituted Into (32-1), these equations 
become = WL (33-1) 




Fig. 28 

This is the condition for the maintenance of steady precession under 
the action of gravity. Writing (33*1) in the form 






we see that for a high velocity of spin the fraction ijfjo) is negligible 
and that the equation thus reduces to 



or 



Wl 



(33-2) 
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regardless of the magnitude of 0. This is the same equation as 
(22-4). It shows that the velocity of precession varies directly as the 
distance I ( = OG in Fig. 28). 

If 

On putting G = M Jc 2 = & 2 , we may write (332) in the equivalent 

' 



ff 

where k denotes the radius of gyration of the gyro about its 
geometric axis. This equation shows that the rate of precession, for 
a given rate of spin, is determined by the shape of the gyro. When 
the gyro is in the form of a flat circular disk, or flywheel, the rate 
of precession is very slow; whereas if the gyro has the form of a 
slender body such as a rod, the rate of precession is very fast. 

For 6 = 90, equation (33-1) reduces to (33-2) regardless of the 
ratio ft/a). 

Since equation (33*1) is a quadratic in ^ 3 there are two rates of 
steady precession. These are found by solving (33*1) to be 

^ ' 

If < < 90, the radicand is less than Co) and both values of ijr 
are positive; whereas if 90 < 6 < 180, the radicand is greater than 
Co), thereby making one value of ijr positive and the other negative. 
It is to be noted that when i/r is positive the precession has the same 
direction as the spin of the gyro, whereas if ^r is negative the 
direction of precession is opposite to that of the spin. 

The roots of (33-4) are equal when the radicand is zero; that is, 



^ ( <n. W . (33 .5) 

C 

Any smaller value of o) will make the radicand of (33-4) negative 
and thereby give imaginary values for ijr. Hence (33-5) gives the 
minimum spin for precession. In this case the value of ^ is, by 
(33-4) and (33-5), 



5-2 
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Solving (33*1) for cos 6, we get 

Cwjr-Wl 



cos 6 = - 



Aft 2 



(33-7) 



TMs equation shows that in steady precession an increase in G causes 

a decrease in 6, and an increase in 

A causes an increase in 6. This 

means that the geometric axis of a 

spinning disk will make a smaller 

angle with the vertical than will 

the axis of a slender rod spinning 

about its geometric axis. TMs fact 

can be verified experimentally by 

using different shapes of c spool 

dancers' and noting their behavior 

when spun. 

The direction of precession in this 
case is determined by the following 
rule: 

The external moment causes the 
positive end of the line I sin 6 to 
precess (rotate about 0) toward the 
positive end of the moment vector, as indicated in Fig. 29. 




Fig. 29 



34. General motion of the top 

To investigate the general motion of a top under the action of 
gravity, we must use equations (32-1). From the third of those 
equations we get 

C(<}> + ft cos 6) = constant = Co), say. (34- 1 ) 

Then, since H S . = C(($> -\-ftcos 6), we have 

H z = Go) = constant. 
We now write equations (32-1) in the form 



(34-2) 



Q, | 

A0+ (Cd) - Aft cos 6) ft sin 6 = W I sin 6. ) 

These are the differential equations of motion of the top. 
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Instead of attempting to Integrate these equations directly, we 
find first integrals of them by other methods.* As the applied force 
W acts vertically downward, there is no external moment about the 
vertical axis OZ-^. Hence the angular momentum about OZ^ is 
constant (Art, 15), The projections on OZ l of the angular momenta 
H X9 H y and H z are (Fig. 23) AtysuPd, and 0&>cosd, respectively. 
Hence we have 

Co) CQS0 + A$ sin 2 B = constant = <7 1? say, (34-3) 

We find another first integral by utilizing the fact that the total 
energy of the spinning top is constant. By (17-1) the kinetic energy 
of a symmetrical gyroscope (BA) is 



In our present problem, (o x = ^r sin d,a) y = () and o) g = &>. Hence we 
have 



The potential energy at any time with respect to a horizontal 
plane through is IFZcos 6. Hence the total energy is 

sin 6) -i + Wl cos 6 = constant = C 2 , say. 

i 2, 

(34-4) 

* First integrals of these equations can be found by direct integration as 
follows: 

Multiply the first equation throughout by sin 6, thus obtaining 



(a) 
This is an exact differential, and integration gives 

^^sin 2 + Ceo cos 6 = C l} 
which is (34-3). 

To find a second first integral, multiply equation (a) by 2^, multiply the 
second of equations (34-2) by W, and add the results thus obtained. Then 
we ave 



Each member of this equation is now an exact differential, and integration 
glVe8 Afc + A$* sin 6 = - 2WI cos 6 + C a , 

which is the same as (34-4). 
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The constants C and G 2 in (34-3) and (34-4) are to be determined 
from the initial conditions of the problem. If the gyro is released 
with a spin velocity ^ , the initial conditions are 

6-6 Q , $ = and fi = 0, when = 0. 
When these values are substituted into (34.-S) and (34-4), we get 



d = Co) cos , C% = + Wl cos 8 g . 

2i 
Hence (34-3) and (34-4) now become 

Ai/r sin 2 6 = Ca)(co8 6 Q - cos 0), (34-5) 

and A(d 2 + ^8m 2 e) = 2Wl(G08d -Gosd). (34-6) 

These two equations determine the subsequent motion of the top. 

Equations (34*5) and (34-6) can be reduced to quadratures by 
the following procedure: 

From (34*5) we get 

/ _ Go)(co& cos 8) __ <7o>(cos 6 cos 6) . 

~ ' ( ] 



Substituting this value of ^ into (34*6), we get 



or 



<9 = ^Vl 2 -^- ^ ( cos 0o ~" cos 5 ) sin2 ^ - G2 ^ 2 ( cos ^o - cos 9 Y}- 
jA. sin u 

(34-8) 
For brevity let us put cosd = u. Then 



--=7 = -=7 or v= -j -- . 
a w smo' ^ 



du A I du 

Hence (34-8) now becomes 

dt A 
whence 



This integral on the left is an elliptic integral, and its inversion 
gives u as a function of t. 
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Then on writing (34-7) in the form 



_ ___, 

we get ^ = 



U 



And since u is assumed to be a function of t from (34-9), we can find 
T]T as a function of t from (34-10). 
Finally, from (34-1) we have 



Hence = ft ,- M 

dt dt 



or 



We can therefore find <f> and ^ as functions of t. 

Since elliptic functions are periodic, the above investigation 
shows that cos 6, 6, i/r and $5 are all periodic functions of the time 
and therefore repeat their values at equal intervals of time. As a 
knowledge of elliptic integrals and elliptic functions is not assumed 
here, we shall investigate the motion of the spinning axle by other 
methods. 

35. The general precession 

To study the precession of the top, we solve (34-5) for $, obtaining 

o cos 0) _ Ck>(cos GQ cos 6) . 

= .4(1- cos 2 0) ' ( } 

This equation shows that tjr and a) have the same sign when 6 > $ 
and opposite signs when 6 < # , which means that the precession is 
in the same direction as the spin in the former case and in the 
opposite direction in the latter case. The equation also shows that 
^=0 when # = # , as it should. 

Equation (35*1) shows that when #># 0; ^ has a finite positive 
value. This means that when the spin axis starts falling, the 
precession starts. The equation shows further that the precession 
has a periodic motion because of the periodicity of cos 6\ 
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On differentiating (35*1) with respect to time, we get 

- 2 cos cos 6 + cos 2 6) 

-~ __ - 



This equation (35-2) shows that $ = when = 0, but that as soon 
as $ becomes different from zero ijr has a positive value, because 
1 2 cos cos (9 + cos 2 6 > ( I cos 0) 2 . Hence $ begins to increase 
at once. The train of events is this: The weight W starts to fall,thus 
increasing 6 and giving rise to 6\ at the same time, ijf starts up at an 
increasing rate. 

36. Nutation of the spin axis in the azimuthal plane 

To find the motion of the spin axis in the azimuthal plane, 
we eliminate ^ between equations (34*5) and (34-6), thereby 
obtaining 

<7 2 w 2 (cos d - cos 0) 2 + J. 2 2 sin 2 6 = 2A Wl(cos G - cos 6) sin 2 6. 

(36-1) 

The maximum and minimum values of 6 are found by putting 
<9 = in (36-1). We thus get 

(cos # - cos 0) [C 2 w 2 (cos - cos 0) - 2A Wl sin 2 6] = 0. 
It follows that 6 = 6 Q and 



Since sin 2 5 is always positive, cos 6 Q cos must be positive and 
6 > Q -i This means that 9 Q is the minimum value of and that the 
spin axis is highest in the azimuthal plane for this value. 
To find the other values of 0, we put 



2AWI ""' 
and replace sin 2 9 by 1 cos 2 6. Then (36-2) becomes 



from which cos 6 = A ^/(l - 2A cos 6 + A 2 ). 
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Now since cos0 <l, the radicand Is greater than 1 2A4-A 2 , or 

( 1 A) 2 . Then for the positive sign before the radical we have 



or cos$>l, 

which is impossible. Hence this value of 6 must be discarded. 

We must therefore take the negative sign before the radical. 
Then we have cos ^ = A _ V(1 _ 2Acos5o+A2)j (36 . 3) 

where x= Swr 

Equation (36-3) gives the maximum value of 6 and therefore the 
lowest position of the spin axis in the azimuthal plane. This 
equation shows that 6$ and 6^ are not Independent of each other. 

From (36-3) we have 

cos X - A = - V( 1 - 2A cos + A 2 ). 

Squaring this and rearranging the result slightly, we get 

cos 2 0-1 2.4TFZ(cos 2 0-l) 



The right-hand member of this equation decreases rapidly as the 
spin velocity increases. Hence cos 6 1 cos # likewise decreases at 
the same rate and consequently 6^6$. The altitude of the 
cycloidal loops (Fig. 31) thus decreases rapidly as the spin velocity 
increases. 

If the maximum value of is to be JTT, thus making the spin axis 
horizontal, we find from (36*3) that 

(3M) 

The angular velocity of the spin axis in the azimuthal plane is 
found from (36-1) to be 

6 = * VI 2 A ^(cos G - cos 6) sin 2 d - (7 2 o> 2 (cos <9 - cos ^) 2 } . 



(36-5) 

Since the radicand is a cubic function in cos 6, cos 6 is an elliptic 
function of t. 
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37. Path of the spin axis on a unit sphere 

Assume that a unit sphere is described about as center, and 
that a tangent plane 5 is drawn through the point P where the axis 
of spin pierces the sphere (Fig. 30). The velocity of P along the 




Fig. 30 



meridian is therefore 6, and its velocity along the dotted latitude 
circle is i/r sin 6. 

Let x denote the angle which the resultant velocity vector makes 
with the meridian-velocity vector 6. Then tan % = ^ sin #/$. Now 
replacing i/r and 6 by their values from (35-1) and (36*5), we obtain 

COS0A COS0 , _. 



This equation shows that % = when 6 = . The point P is therefore 
moving along the meridian when 6 = # . 
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If % = 90 ? so that the motion of P is entirely along the dotted 
parallel circle, then tan i/r co and thus the denominator of the 
fraction in the right-hand member of (37-1) must be zero. Hence 

TfTF* M A Ve 



. <7 2 6> 2 (COS0 -COS0) 

m *- - 2AWI - ' 

which is (36-2). Hence P is moving horizontally on the lower 
circle 6 = X . 

By (35* I ), ^ = when 6 = $ ;. and for 6 = & the value of i/r is, by 
(35-1) and (36-2), 

cos X ) __ Co)(cos 6 6 cos ^j) 

~~~~~ ~ 



"" Co) * 

This angular velocity of precession on the lower circle 6 = 6^ is twice 
the value found in (33-2) for regular precession. 

The preceding discussion has brought out the following facts 
concerning the motion of the spin axis on the unit sphere: 

(1) It reaches its maximum height on the circle = # and is 
moving along the meridian when reaching that circle. 

(2) It reaches its lowest point on the circle = 6^ and is moving 
horizontally at that point. 

(3) Equations (35*1) and (36-5) show that the motion is periodic 
in the azimuthal plane and also along the latitude circles. 

These facts enable us to trace the path of the spin axes on the unit 
sphere as shown in Fig. 31. The path consists of a periodically 
recurring series of loops which are sometimes called spherical 
cycloids. A numerical discussion of these loops will be given at the 
end of the next article. 

38. Periods of the cycloidal loops. Pseudo-regular precession 

In order to study the motion of the spin axis in more detail, we 
go back to equation (36-5) and write it in the form 



/rt A TTr , . on ^o <>/ * ^ 

dt ~~*~~ ( ~ CW(cos 6 - cos 6)). 



76 THE GYROSCOPE 

We simplify this by putting cos 8 = u, from which 



_ __ __ 
dt "" sin dt * 

Then the above equation becomes 




Fig. 31 



We simplify it still further by putting 

2WI , C 
a-j-, 6 = ^. 

Then the equation becomes 



(38-1) 



Since the right-hand member is a polynomial of the third degree 
in u, an attempt to integrate (38-1) will lead to an elliptic integral 
and therefore the inversion of it will give u as an elliptic function 
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of I; and since elliptic functions are periodic, the rise and fall of the 
spin axis is periodic. We can avoid the elliptic functions., however, 
when the spin velocity of the top is high. 

Since (dujdtf is a positive quantity, the right-hand member of 
(38-1) must also be positive. The first factor % u is positive 
because # < 6, and consequently cos 6$ > cos 6, or n > u. The factor 
in brackets must also be positive, but this fact is not evident at 
sight. In order to get a better picture of the situation, we construct 
the graph of equation (38*1). 

m. 



Fig. 32 

Let us put (duldt)*=*f(u). Then 

f(u) = a(u Q -u)(l~<u?)- 6 VK - uf. (38-2) 

Now since au* is the dominant term for large values of u, we have 
/(-co) =-oo, /(oo) = oo, 



Furthermore, since cos# must lie between 1 and 1, and since 
f(u) is positive, there must be a region between 1 and 1 where/(it) 
is positive. These facts enable us to determine the general form of 
the graph off(u), which is shown in Fig. 32. 

The three roots are thus all real, and the locations of u^, u 1} u^ 
are as shown on the graph. Here % is the value 



78 THE GYROSCOPE 

and % is the value A + */( 1 2A cos Q + A 2 ) found in Art. 36. These 
two roots can be found by equating to zero the bracketed expression 
in the right-hand member of (38*1). On putting that expression 
equal to zero, we have 

< 38 ' 3 > 



This equation shows that when a) is large, u u is very small and is 
of the order 1/oA This fact enables us to find 6 and ^ as functions of 
the time t. 

As 6 must be very near to 6 when a) is large, we may put 



2 



where y is an arbitrary variable. Then cos# = cos(# + ?//to 2 ). We 
now expand cos(0 + ^/&> 2 ) by the Taylor formula: 

f m (fl \ 

s ' (38 ' 5) 



and retain only the first two terms of the expansion. We have 
f(0) = cos (9, hence /($ ) = cos 6 Q ; 

/'(#)=- sin 6>, hence /'(0 )= -sin0 ; 

f(0) = - cos 0, hence /"(# ) = - cos . 
Now since A = 9//w 2 , we have by (38-5) 







Hence 

v 0)* ' " | 

(38-6) 



Substituting into (38-1) these values of u and dujdt, replacing 
1 ul by sin 2 , and simplifying slightly, we get 
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As the term 2aw ^ 2 is negligible in comparison with the other two 
terms (because it does not contain the large factor 2 ) ? we neglect it 
and then have 



from which -- - 



Integration of this gives 

. , / 
sm-i 

\ 

Now 5 = 9 Q when $ = 0. Then by (38-4), 97 = when * = 0. Hence from 
the above equation we get 



Then 
from which 



j: 
= sm 6* ~ i^ = - 



Solving for ?/, we get 

(38-7) 



Then by (38-4) we get 

(38 ' 8) 



The period of this motion is 

*-*?-*, (38-9) 

cob o)C \ > 

which evidently decreases as the spin velocity increases. 
To find ^ as a function of t, we have, from (35-1), 

_ Ca)(co8 6 Q cos 6) __ CV(UQ u) 

*~ 4(1-008*5) "" 4(l-tt 2 ) * 

On replacing u by its value given in (38-6), we have 



" A (o> 2 sin 2 + 2u<fl sin ) * 
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Neglecting the term 2w ^sin0 because of its smallness in com- 
parison with, (jt) 2 sin 2 , we get 

, Cj = by 

Ao) sin OQ a) sin 6 Q ' 

since Of A = b. Now replacing ^ by its value given in (38-7), we have 



(38-10) 
Integration now gives 



= -tt -- -sin cobt], 
Ca)\ o)b ) 



the constant of integration being zero because ^ = when = 0. 
Equation (38-10) shows that the angular velocity of precession is 
the algebraic sum of the steady precession WljCa) (the same as 
given by (33-2)) and the periodic velocity (Wl/Ca))coso}bt, This 
periodic term causes the velocity of precession to be slightly un- 
steady, and for this reason the precession is called pseudo-regular 
precession. 

Our results for the general motion of the axis of spin are thus 
expressed in the two equations 

l-ooB06f), (38-8) 

Wl /I \ 

i/r^^-lt -- rsin^Ztf , (38-11) 

r Co) \ o)b ) 

where b = C/A. These are the parametric equations of the spherical 
cycloids in Fig. 31. 

Numerical Example. In order to get an idea of the magnitude 
and frequency of the oscillating motion of the end of the axle of a 
spinning gyroscope, let us consider a steel disk 6 inches in diameter 
and 1 inch thick, mounted at the center of a 12-inch spindle of 
negligible weight. The spindle is held so as to make an angle of 15 
with the vertical, its lower end resting in a smooth bearing. If the 
gyro is given a spin of 7500 r.p.m. and the spindle is free to move in 
any manner about its point of support, find the rate of precession 
and the period and amplitude of the upper end of the spindle. 

Solution. Before proceeding with the numerical calculation, we 
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shall derive formulas for -4, and A in terms of the given quantities. 
Denoting the mass of the disk by M, its radius by r, and its thickness 
by A, we have 

Jfr* 






Substituting these values into the formula, 



we get 



Since, by (36-3), 6^ depends only on and A, we have the interesting 
result that 6-^ is independent of the mass of the gyro. This means 
that the lower circle for a gyro weighing a ton is the same as for 
one weighing an ounce, if the upper circle is the same in both 
cases. 

We now proceed with the numerical work. The pertinent data in 
this problem are 

(j) = 7500 r.p.m. = 125 r.p.s. = 2507T radians per second, 



Z = |foot, r = ifoot, h =tV foot, g = 32-16ft./sec. 2 
Substituting these data into the expression for 2A above, we have 

_ 



/. A = 35.18194. 

Now substituting into (36*3) the values of 6 Q and A, we have 
008(9! = 35. 18194 -34.21699 = 0.96495. 
/. 6> 1 =1512 / 52 // . 
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The amplitude of nutation in the azimuthal plane is thus only 
X - = 12' 52" = 0.00374 radian. 

The arc distance between the parallel circles and 6 on the unit 
sphere is thus 

s = r(0! - ) = 1 x 0.00374 = 0.00374 foot 
= 0.0449 inch=^ inch. 

The top of the spindle thus oscillates in a zone % inch wide on the 
unit sphere and deviates only about ^ inch from its mean pre- 
cessional path. 
From the data of the problem we have 

M 

(12Z 2 -j-3r 2 + A 2 ) 

jui \.2t X / 

C = M^ = 6\ 



= 230 
6 1 (i) 2 ) 27" 

Hence by (38*9) the period of the oscillations is 



which means 14.7 vibrations per second, or about three times as 
fast as the ticking of a watch. 

By (38'iO) the steady or constant part of the angular velocity of 
precession is 



a 



2(32 16) (4) 

.\* n j^\' = 0.655 radian per second. 
(|) z (2507T) 

This means one revolution in 9.6 seconds, or 6.26 r.p.m. 
The number of cycloidal loops on the unit sphere is thus 



0.068 
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The circumference of the upper circle is 

277(12 sin 15) = 19.5 inches. 
Hence the width (base) of each loop is 

19.5 



141 



= 0.1 38 inch. 



Since the base of the ordinary plane cycloid is TT times its maximum 
height, the base of a cycloidal loop 0,0449 inch high would be 
0.04497T, or 0. 141 inch. The spherical cycloids of the present example 
are therefore slightly more slender than plane cycloids. 

The results of this example show that the axle of a spinning top 
precesses with a rapid tremor of small amplitude, the tremor being 
hardly noticeable to the naked eye. 

39. The sleeping top 

When a top is spun rapidly with its axis nearly vertical and then 
released, it soon settles down to a steady spin with its axis seemingly 
vertical. The top in this apparently motionless state is said to be 
'sleeping'. We now proceed to determine the necessary condition 
for the top to sleep. 

On substituting into the second of equations (34-2) the value of 
i/r given by (35-1) and then replacing sin 2 6 + cos 2 6 by 1, we get 

Now since the spin axis of the top is practically vertical when the 
top is sleeping, we put # = 0. Then (39-1) becomes 



- =m ^ 

3 v 



Here 6 is a very small angle, and this fact enables us to simplify the 
equation still further. Erom the series expansion 



6* 
we get (1 cos 0) 2 = .. 



6-2 
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Then after replacing sin 5 by 6, equation (39-2) becomes 



, A A /OAOX 

or $ + - 5 - = 0. (39-3) 

This is a harmonic motion equation, the solution of which is 
= G I cos kt + C 2 sin K, 

, 7 

where A 



Then we get 0^-C^Jc sin kt + O^k cos R (39-4) 

To find the value of ^ for a sleeping top we put = in equation 
(35-1). Then we have 



or ^ = T/, ^ /A- ( 39 * 5 ) 

r -4(1 + cos 0) x ' 

Equations (394) and (39-5) show that when a top is sleeping the 
free end of the spin axis moves with simple harmonic motion in the 
azimuthal plane and that it precesses with a periodic motion about 
the vertical. It is thus a stable motion about the vertical position. 
If we let 0->0 in (39-5), we see that i/r^Ccj/^A, a constant. The 
irregularity in the precession thus fades out as the spin axis becomes 
vertical. 

The period of the motion of the sleeping top is 

P- 27r - 

P ~~k- 

If this period is to be real, that is, if the spin axis is actually to 
oscillate and the motion is to be stable, we must have 



or > 2 >. (39-6) 

Equation (39-6) thus gives the velocity of spin necessary for a stable 
sleeping top. A smaller spin velocity will allow the top to wobble 
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and eventually fall. This equation shows that if is In 

comparison with A , the spin velocity must be very high in order for 
the top to sleep. This fact explains why it is almost impossible to 
spin a pencil on its point. 

For the case of a sleeping top the graph off(u) has the appearance 
shown in Fig. 33, the roots U Q and % each being equal to 1. 



Fig. 33 

Historical note 

The early history of the gyroscope is rather obscure. Probably 
the earliest gyroscope of the type now in use was constructed by 
Bohnenberger, of Germany, about 1810. An improved type was 
constructed by Walter Rogers Johnson, of Philadelphia, Pa., in 
1832, and used by him to illustrate the dynamics of rotating bodies. 
Johnson called his device a 'rotascope'. 

In 1852 Leon Foucault, of Paris, constructed a very refined 
gyroscope to show the rotation of the earth. In order to bring the 
rotor to very high speed and keep it spinning for several minutes, 
Foucault constructed a gear train machine on which the gyro was 
placed for speeding up. When the speed was sufficiently high, the 
spinning gyro was lifted from the machine and placed in its gimbals. 
Because of the fact that the refined device exhibited the rotation of 
the earth before his eyes, Foucault named it the 'gyroscope', from 
the Greek words gyros (yvpos), circle or ring, and sJcopein 
to view. 
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In 1878 George M. Hopkins, American inventor, applied the 
electric motor to the gyroscope and thereby kept it spinning at any 
desired speed as long as desired. Gyroscopes were not put to 
practical and industrial use until the first decade of the twentieth 
century, at which time the gyroscopic compass, the ship stabilizer, 
and the monorail car were all invented. 

The mathematical foundations of gyroscopic theory were laid by 
Leonhard Euler in 1765, in his Tkeoria motus corporum solidorum sue 
rigidorum. 



PART II 

APPLICATIONS OF THE GYROSCOPE 



CHAPTER V 

Gyroscopic Action in Vehicles and 
Rotating Bodies 

Because of the stubborn tendency of a rapidly spinning gyro to 
maintain its plane of rotation, gyroscopes are used as direction 
indicators and as stabilizers. The forced precession of the gyro axle 
also produces interesting reactionary effects in rotating machinery 
and in the astronomical world. 

40. Gyroscopic effects in car wheels rounding a curve 

When a vehicle such as a train, automobile or racing car is 
moving at high speed along a straight road, the gyroscopic effect of 
the wheels is to keep the vehicle moving straight ahead. But when 
an external moment is applied to force the vehicle to change its 
direction, gyroscopic reaction immediately comes into play. To 
determine this gyroscopic reaction in the case of vehicles, let us 
consider a pair of wheels connected by a straight axle. This wheel- 
axle ensemble may be considered a gyroscope. When the vehicle 
is forced to round a curve, such rounding is a forced precession of the 
wheel-axle ensemble around a vertical axis through the center of 
curvature of the curved track. For simplicity we shall assume 
that the roadbed is not banked and, in the case of a railroad, that 
the outer rail on the curve is not elevated higher than the inner 
rail. 

Let Fig. 34 represent an axle and its two wheels rounding 
a circular curve of radius R with speed v. Let O be the center 
of curvature of the curve, and let r denote the radius of the 
wheels. If co denotes the angular velocity of the wheels and Q 
denotes the angular velocity of precession about the vertical axis 

OF, then 

(40-1) 



and v = ra). (4O2) 

89 
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By equation (27-2), with 6 = 90, the gyroscopic reaction moment 

due to the forced precession is 

A r = (7o>a (40.3) 

where C denotes the moment of inertia of both wheels about their 
common axle and K is the reaction couple (moment) acting in the 



ft 




Fig. 34 

vertical plane through OF and the axle. Substituting into (40-3) 
the values of O and a> from (40-1) and (40-2), we get 



(40-4) 



Jfr- 



The gyroscopic reaction thus varies as the square of the velocity of 
the vehicle along the curve. 

Since gyroscopic phenomena depend only on rotation and not at 
all on translation, and since the angular velocity of precession of 
the ensemble about a vertical axis through the mid-point of the 
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axle Is exactly the same as its velocity of precession about Of 7 ", we 
may consider all the gyroscopic phenomena as taking place with 
reference to the mid-point of the axle as origin. Then by the right- 
handed-screw rule the precession vector must be drawn vertically 
downward; and since the gyroscopic reaction moment tends to 
turn the spin vector into coincidence with the precession vector 
(Art. 28), we see at once that the reaction vector K must be drawn 
backward, as indicated in Fig. 34. This shows that the effect of the 
gyroscopic reaction is to decrease the pressure on the inner raJI or 
wheel and increase it on the outer rail or wheel. The gyroscopic 
moment is thus in the same 'direction as the moment due to centri- 
fugal force and aids the latter in its tendency to overturn the 
vehicle. 

In order to compare the gyroscopic and centrifugal moments, let 
us consider a vehicle having n pairs of wheels. Let 

Jf=mass of entire vehicle, including wheels. 
m = mass of each wheel. 
&= height of e.g. of vehicle above roadbed. 

r = radius of each wheel ( = distance from center to rail tread in 
the case of locomotive and train wheels). 

c = a constant such that cr= radius of gyration. For most 
wheels, 1/^/2 < c < 1. 

Then for a single wheel the moment of inertia about the axle is 
mcr 2 , and for n pairs of wheels the moment of inertia is 

C=2mncr*. 

By (40-4) the gyroscopic reaction moment for n pairs of wheels is 
therefore 



Since the moment of the centrifugal force about the outer rail is 

Mv 



K 2mncr /2mn\ 
have _- __. ^__- . (40 . 5 ) 
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Now since the mass of the wheels Is less than the mass of the entire 
vehicle and h is greater than r, we have 2mnjM< 1, r/A< 1, c< 1. 
Hence ^ 

N <1 > 
or K<N. 

Note that the ratio KIN is independent of the speed of the 
vehicle and of the radius of the curve. In all ordinary vehicles the 
gyroscopic moment is so much less than the centrifugal moment 
that the effect of the former need not be considered. In auto- 
mobiles, for example, the gyroscopic moment constitutes less than 
5% of the total overturning moment when rounding a curve. 
Calculations from known data for one electric locomotive showed 
that the gyroscopic moment was only 6% of the centrifugal 
moment, or about 5f % of the total overturning moment. 

Although the gyroscopic and centrifugal moments act in the same 
direction, the forces of those moments do not have the same 
direction. The forces due to gyroscopic moments act vertically 
downward on the outer rail or wheel and vertically upward in the 
case of the inner rail or wheel, whereas the centrifugal force acts 
horizontally outward on both wheels. 

In order to make a quantitative comparison of the gyroscopic and 
centrifugal moments on a vehicle, let us consider the extreme case 
of a racing car. Assume that the wheels are 3 feet in diameter and 
weigh 100 pounds each, that the combined weight of car and driver 
is 1800 pounds, that the center of gravity of car and driver is 
20 inches above the ground, and that c = 0-9. Then 

r = f feet, & = ff=f feet, Jf=1800/^, m*=IQQ/g 3 n = 2. 
Substituting these data into (40-5), we get 

| 

' 



N~ lQ\lBOOIg 
or ^" = 



The gyroscopic moment is thus less than one-fifth of the centrifugal 
moment. 

Another gyroscopic effect in automobiles is that due to the 
flywheel. Since the axis of rotation of the flywheel is along the 
longitudinal axis of the chassis and since the flywheel rotates to the 
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left as viewed from behind, the spin vector must be drawn 
along the axis of spin. When the car makes a turn to the right or to 
the left, the axle of the flywheel is forced to precess in the 
direction about a vertical axis. This precession brings into play a 
gyroscopic reaction moment about a transverse horizontal axis, 
The effect of this reaction moment is to increase the downward 
pressure on one end of the flywheel axle and increase the upward 
pressure on the other end, thus tending to bend the axle in a vertical 
plane. The gyroscopic moment hi this case is therefore a bending 
moment. 

The magnitude of this gyroscopic moment is given by the usual 
formula 



where o) is the angular velocity of the flywheel, O is the angular 
velocity of precession, and C is the moment of inertia of the flywheel 
about its axis of spin. The gyroscopic reaction moment due to the 
precession of the spin axle of the flywheel in automobiles is usually 
of minor importance because of the relatively small mass of the 
flywheel. 

41. Derivation of the differential equations of motion of a gyroscope 
by application of the principle of gyroscopic reaction 

In the applications of the gyroscope the differential equations of 
motion can be derived in either of two ways: (1) by determining the 
angular velocities a) x , o) y , o) z about the axes of the moving trihedron 
and then using equations (16-4), or (2) by determining the G.R.M.'S 
and then substituting them in the fundamental rotational equation 
(20-2). One method is preferable in some problems and the other in 
others. The G.R.M. method is the shorter, but is less foolproof. The 
reader is advised to use the first method until he feels safe in using 
the second. 

As an illustration of the G.R.M. method, let Fig. 35 represent a 
gyro spinning about a horizontal axis and let the axle be acted on 
by a horizontal couple F, F. Let us consider the moments about 
the axes OZ and OX. 

(a) For OZ we have first the moment M due to the external 
couple F, F. This couple tends to increase the angle <j> and is 
therefore positive and must have the same sign as d 2 <f>]dt 2 . The 
couple also causes the axle of spin to precess about the horizontal 
axis OX as indicated in the figure. This forced precession about OX 
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induces a gyroscopic reaction moment G^dd/dt about OZ, this 
moment being negative according to the rule laid down in Art. 28. 
Hence the reaction vector K must be drawn downward. The total 
moment about OZ is therefore M Co) dd/dt. Then from the funda- 




Fig. 35 

mental moment equation (20-2) we have for the motion about OZ 
the differential equation 



or 



d6 



where A denotes the moment of inertia of the gyro about OZ. 
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(6) As. to the moments about OX, we see that there is BO external 
moment. The first effect of the external couple F, F is to an 

angular velocity d$jdt about OZ. This forced rotation (a precession) 
will immediately induce a gyroscopic reaction moment of magnitude 
Ca>d$/dt about OX, and by the rule of Art. 28 this moment is 
positive. Hence from the fundamental equation for rotation about 
an axis we have 



ft 

*---<>. 

for the motion about OX. 

The differential equations for the motion of the gyro under the 
action of external couple F, F are therefore 



T 

dt 
dd 



where M is the external moment tending to deflect the axis of spin. 
If we derive the differential equations of motion for this example 
by using equations (16'4), we have 



On substituting these into (16-4), we have 






as before. 



96 



THE GYROSCOPE 



42. Gyroscopic effects of dynamos, motors and turbines inslaEed on 
ships 

The rolling, pitching and yawing of a ship on a rough sea can 
cause gyroscopic action of motors and turbines in three mutually 
perpendicular directions. Boiling is rotation about the longitudinal 




Fig. 36 



axis of the ship, pitching is rotation about a transverse axis, and 
yawing is rotation about a vertical axis, all these axes intersecting 
at the center of gravity of the ship as origin. Let us consider the 
gyroscopic effects of these actions on a motor or a turbine when its 
axis of spin is (a) parallel to the longitudinal axis of the ship, 
(6) parallel to the deck and perpendicular to the longitudinal axis, 
and (c) perpendicular to the deck. Bearing in mind that each of the 
above motions of the ship will cause a forced precession of the axis 
of spin, applying the G.K.M. method of the preceding article or the 
formulas of Art. 27 and rule of Art. 28, and referring to Fig. 36, we 
obtain the following results: 
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(a) Axis of spin parallel to longitudinal of : 

Rolling: no gyroscopic action. 

Yawing: K t = Co^. 
Pitching: K S =-C(M-. 

Wb 

(b) Axis of spin transverse to axis of ship : 

Rolling: K z = Cw^f. 
tftt 

Yawing: K y =-cJ^. 

Pitching: no gyroscopic action. 

(c) Axis of spin perpendicular to deck : 

Rolling: K 9 =-aA. 
d/t 

Yawing: no gyroscopic reaction. 
Pitching: K v = Cv^. 

In the above formulas G denotes the moment of inertia of the gyro 
about its axis of spin. 

Since rolling is the most common form of motion of a ship, it is 
very important that motors and turbines be installed with their 
spin axes parallel to the longitudinal axis of the ship, for in that 
position rolling does not cause any gyroscopic action. 

To get a quantitative estimate of the gyroscopic moment due to 
rolling and pitching, let a denote the angle which the ship mast 
makes with the vertical (the plumb line) at any time t, and let <I> be 
the maximum numerical value of a. Then since rolling and pitching 
are periodic phenomena, we may write 

a = sin# (42-1) 

where y? is an auxiliary variable such that /?= JTT when a = $. Then 
for the angular velocity of rolling or pitching, we have 

. M) 
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In terms of $ the period, T 9 of rolling is given by the well-known 
formula 

^Jjjdt' 

_.- dB 2n 

Hence -^7 = "m- 

Substituting into (42-2) this value ofd/t/dt, we get 

g = ^ cos/? . (42-3) 

This angular velocity has its maximum value when /? = 0. Hence 

* < 42 ' 4 > 

Now since the rolling and pitching of the ship forces the axles of 
dynamos and turbines to process with the same angular velocity, 
we have from (42*4) - 

fl = ^p, (42-5) 

where O denotes the angular velocity of precession. 

If O denotes the moment of inertia of the turbine rotor about its 
axis of spin and o) is the spin velocity, we have 

(42-6) 



for the gyroscopic reaction moment due to the rolling or pitching of 
the ship. 

Numerical Example. Let us assume that a steam turbine weighs 
3 tons and is rotating at 3000 r.p.m. or lOOyr radians per second. 
Then 

(a) If the turbine were installed with its axis at right angles to the 
longitudinal axis of the ship, the gyroscopic reactions would be due 
to the rolling of the ship. If we assume an angle of roll of 20 to 
either side from the vertical and that the period of rolling is 
8 seconds, then 

= 20 = -radians 
9 

and T = 8 seconds. 



GYROSCOPIC ACTION IK ROTATING BODIES 99 

If the average radius of the turbine rotor is 1 foot and the moment 
of inertia is assumed to be fiM> 2 ? then by (42-6) 



K=. w v ' V - = 12,036 lb.ft. 

b 

This reaction moment reverses its direction every 4 seconds. It acts 
in a plane parallel to the deck of the ship and tends to skew the axle 
of the turbine back and forth ever 4 seconds, thus exerting a 
pulsating pressure against the sides of the bearings. Since the rotor 
is rigidly keyed to its axle, the axle is subjected to an alternating 
bending moment which changes its direction (with respect to a 
diametral plane of the axle) 3000 times per minute. 

(b) If the turbine be installed with its spin axes parallel to the 
longitudinal axis of the ship, as it should be, then the only gyroscopic 
action of any importance that will occur is that due to the pitching 
of the ship. The angle of pitching is less than the angle of rolling and 
will be taken as 10 in the present example. The period of pitching 
will be taken as 10 seconds. Then 

< = 1 = radians, 
18 

T=10 seconds. 
Hence by (42-6) we have 



K- . -- 48 i51 b . fl . 

The gyroscopic moment due to pitching is thus very much 
smaller than that due to rolling, because of the smaller angle in the 
case of pitching. 

(c) If the ship turns completely around (through 180) in 
1 minute, the angular velocity of precession is O=^7r radian per 
second. Hence the gyroscopic reaction moment in this case is 



This moment acts in a plane perpendicular to the deck. If the 
turbine is installed with its axis parallel to the longitudinal axis of 

7-2 
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the ship, the above reaction moment will act in a perpendicular 
plane which contains the longitudinal axis and will raise one end 
of the ship and depress the other. For example, if the turbine is 
rotating to the right as viewed from the rear end of the axle, the 
gyroscopic moment will raise the stern if the ship turns around to 
the right and raise the bow if it turns around to the left. Similar 
gyroscopic effects can be observed when making a quick turn with 
a motor boat. 

43. Rolling hoops and disks 

A rolling hoop or disk (coin, for example) may be looked upon as 
a gyroscope whose center of gravity has an irregular motion of 
translation and whose rotor rolls on a rough surface without 
slipping. Such a motion is more complex than that of the usual 
type of gyroscope, and its mathematical treatment is considerably 
more complicated. Inasmuch as the mathematical treatment of 
the most general case of hoop motion would take us beyond the 
scope of this book, we do not enter into it. Instead, we confine our 
treatment to some of the simpler and more familiar aspects of the 
motion.* 

Everyday observation shows that if a hoop is started rolling on 
the ground in a vertical plane, it will steer itself for a time as it 
encounters small obstructions in its path. Thus if it passes over a 
small obstruction which tilts it from a vertical plane, the moment 
of its weight immediately causes the hoop to precess about a 
vertical axis through the point of contact (about a vertical diameter, 
practically) and thus steers it in the direction in which it started to 
fall. This change in direction sets up a centrifugal force which tends 
to bring the hoop back to a vertical position and allow it to roll on 
in the new direction. 

To study the behavior of a rolling hoop, let us assume that it is 
started rolling with constant angular velocity o). Tig. 37 represents 
it in a slightly inclined position after encountering a small obstacle in 
its path. This slight tilt causes it to travel temporarily in a curved 
path. Let MN denote the tangent to the path at the point of 
contact of the hoop, and let d^> denote the horizontal angle through 

* The interested reader will find a good treatment of the general case in 
Appell's Traiti de Mfoanique Hationnelle, tome 2 (third edition, 1911), 
pp. 248-54. 
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Fig. 37 



which the plane of the hoop turns in time dt. Then the angular 
velocity of the precession of the hoop about a vertical axis is d<f>jdt. 

Let 

a) = angular velocity of hoop about its geometric axis, 

v = linear velocity of point of contact 

= linear velocity of e.g. of hoop, 
r = radius of hoop, 
R = radius of curvature of path, 

6 = angle of inclination of plane of hoop with the vertical, 
m = mass of hoop, 

G = moment of inertia of hoop about its geometric axis, 
A = moment of inertia of hoop about a diameter, 
/ Mjv = moment of inertia of hoop about MN. 
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Then v=*ra) (43-1) 

and v = JR^' (43*2) 

at 

Since ce> is assumed to be constant, v must also be constant. Hence 

JJ-^- = a constant. (43-3) 

at 

Now let us consider the moments acting about the line M N when 
the plane of the hoop is inclined at a small angle 6. The moment of 
the weight is mgr sin 0, which is practically mgr 6, since 6 is assumed 
to be small. This moment tends to increase 0. 

The moment of the centrifugal force is 



- 5 
M at at at 

and this moment tends to decrease 0. 

The precession of the hoop about a vertical axis induces a 
gyroscopic reaction moment of magnitude Ca)(d(j>ldt). This moment 
also tends to decrease 6. Hence from the fundamental moment 
equation (20-2) we have 



We next consider the moments about a vertical diameter through 
the point of contact. The only moment about this line is the 
gyroscopic reaction moment tending to increase $ to keep the hoop 
from falling. This moment is Co)(d0ldt), and it has the same sign as 
Hence by the fundamental moment equation we have 



The differential equations for the dynamical equilibrium of the 
hoop are therefore 



(43-4) 
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These equations will be solved with the initial conditions 6=0 
and 6 = 6 Q when t = 0. In this case the hoop starts out in a vertical 
plane. Since (D is assumed to be constant, we can integrate the 
second equation at once, thereby obtaining 

A~ C(j)6 = constant == O lt say. 

Cut 

If d$/dt = when & 0, we get C- L 0. Hence we have 



. (43-5) 

at 

This equation states that the rate of change of direction in its 
path is proportional to the inclination of the hoop ? and this, in 
view of equation (43-3), means that the more the hoop or coin is 
inclined, the greater is the curvature of its path. The behavior of 
rolling coins agrees with this statement. 

On substituting in the first of equations (43-4) the value of 
d$/dt from (43-5), we get 



For a hoop we have 

mr 2 



Substituting these values into (43-6), we get 



The auxiliary equation of this differential equation is 



or 



(43-8) 



If the motion of the hoop is to be stable in a vertical plane, so that it 
will oscillate about a vertical position as it rolls along, the values of 
a must be imaginary. Hence we must have 
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(43-9) 



If the spin velocity becomes less than J ^(g/r ), the hoop will begin 
to wobble and will soon fall over. 
From (43-8) we have 



say. 
Consequently, the general solution of (43-7) is 



Then 6 = C 2 Jc cos M C z Jc sin H. 

Now substituting the initial conditions into these equations, we 
find 



Hence the values of 6 and 6 are 

(43-10) 

(43-11) 
The frequency of oscillation about a vertical position is 



and thus varies practically directly as the spin velocity. 

To find $5 we substitute into (43-5) the value of 6 from (43-10), 
giving 



dt L# j 

Integrating this, we get 
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To find C t , put $ = when t = 0. Then 



Hence we get 



For a rolling eKsfc or com we have 

C = |wr 2 , J. = |wr 2 , /^ = mr 2 -f %mr 2 = fmr 2 . 
Substituting these into (43-6) and simplifying, we get 



-O. (43-14) 

Here the auxiliary equation is 



If the values of oc are to be imaginary, we must have 



or fc> 2 >. (43-15) 

3r 

On comparing this value of co 2 with that' given by (43*9), we see that 
the spin velocity of a disk or coin must be greater than for a hoop 
if they are both to roll on without toppling over. This result agrees 
with everyday observation. 

44. Gyroscopic grinding nulls 

The principle of gyroscopic reaction has been applied in two types 
of grinding mills where intense pressure is desired. The two types 
are known as the edge mill and the Griffin or pendulum mill. Both 
types grind by crushing. 

(a) Edge Mills. The edge mill is a gyroscopic grinder used for 
crushing ore, seeds, etc. It consists of a large steel pan in which 
heavy rollers, called mutters, roll on the bottom of the pan and at 
the same time revolve about a vertical shaft passing through the 
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center of the pan. The mullers (usually two placed symmetrically 
about the shaft, but sometimes only one) are in the form of heavy 
flywheels with thick rims. They rotate about axles (either horizontal 
or inclined) which are attached to the vertical shaft by means of 
short cranks that permit the mullers to have enough vertical motion 




Fig. 38 
(Photograph courtesy National Engineering Co.] 

to roll over uneven lumps in the pan. The material to be ground by 
crushing is deposited on the bottom of the pan, through a hopper, 
and the mullers are forced to roU over it without slipping. A kind 
of plow is rigidly attached to the vertical shaft and keeps the 
material raked into the path of the mullers (see Fig. 38). 

In the following investigation of the gyroscopic action in an edge 
mill, * we consider the general case in which the axle of the muller 

* The gyroscopic theory of edge mills was first given by B. Grammel in 
1917. 
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Is inclined at any angle to the vertical shaft. Since the results for a 
single muller apply equally to the case of two millers, we confine 
our discussion to the case of one muller. 

Fig. 39 represents the general case of one muller. Since the muller 
axle is not horizontal, the face of the muller must be a frustum of a 
right circular cone in order to rest flatly on the bottom of the pan. 




'////^///////////^ 



Fig. 39 




If the vertical shaft is rotated at constant speed and the muller 
rolls without slipping, the gyroscopic reaction of the muller on the 
bottom of the pan is given by equation (27-1), in which A now 
represents the moment of inertia of the muller about an axis 
perpendicular to OG and passing through the fixed point 0. Since 
by hypothesis the muller does not slip on the plate (bottom of the 
pan), the instantaneous axis of rotation is the line passing through 
the fixed point on the vertical axis and the point IS, which may be 
taken as the point of contact of the muller on the plate. Then if 6 is 
the angle between the axle and the downward vertical, counted 
positive from the downward vertical to the axle, and a is the angle 
subtended at by the mean radius GE of the muller, we have the 
vector parallelogram as shown in Fig. 40. 

Applying the law of sines to this figure, we get 



Q. 



sin(# a) sin a' 



(44-1) 
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Since rotation about the vertical axis is to be considered the 
independent motion in this case, we eliminate $5 between (44-1) and 
(27-1). Thus, from (44*1) we have 



sin (8 a) 

sina 



(44-2) 




Fig. 40 

Substituting this into (27-1) and disregarding the negative sign, we 
have 



sma 



J 



Expanding the factor sin(#~a) and simplifying the result, we get 
K = O 2 [<7 cot a sin 2 - \A sin 26]. (44-3) 

If the muller is rotating to the right as viewed from 0, it will roll 
to the right (clockwise) around the vertical axis as viewed from 0. 
Hence the vectors representing <j> and O will be drawn as shown 
in Fig, 39. Then since the gyroscopic reaction moment tends to 
turn the spin vector into the direction of the precession vector, the 
gyroscopic reaction exerts a downward pressure on the grinding 
plate. This reaction moment is thus in the same direction as the 
moment due to the weight of the muller. Since the moment due to 
the weight of the muller is W I sin 6, the total crushing moment is 

M = Wl sin 6 + a 2 (C cot a sin 2 6 - \A sin 26). (44-4) 
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A i 2/ji, I~eos2i9 , 

On replacing sin* & by , we nave 
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Ccota 



The terms in pa.renth.eses can be telescoped into a single term by the 
well-known procedure of multiplying and dividing by the square 
root of the sum of the squares of the coefficients of sin W and 
cos 20. This procedure gives 



02 

~ 
2 



sin 26 



Ccota 



<7O 2 cota 



or 



(44-5) 

where tan /? = (G/A) cot a (Kg. 41). Since tan/? is necessarily positive, 
ft must be less than 90. 

In order to show just how the total moment M varies with the 
angle 6, we construct the graph of equation (44-5). That equation 
shows that the value of M for any given value of 6 is the algebraic 
sum of the three quantities: 

(a) M^ 

(b) Jf 2 = - 



(c) 
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Hence If we construct the graphs of (a), (b) and (c), we can obtain 
the ordinate M of (44*5) for any value of 6 by adding the ordinates 
of (a), (b) and (c) for that value of 6. The period of (a) is 2n and that 
of (6) is n. Since Jf 3 is a constant, its graph is simply a straight line 
parallel to the $-axis. 

When W, I, Q, A, (7, a and ft are literal quantities, it is not possible 
to construct the composite graph of (44-5) without additional 
information. We need the values of If and dMjdd for certain values 
of 0. Consequently, we next find dM/dd and then make a table of 
values of M and dM[dt at 45 intervals from # = to = 1 80. 
Also dMjdd for 6 - 135 - J/?. 

From (44-4) and (44-5) we get 



dM_ 
dd 



= Wl cos 6 + Q 2 (0 cot a sin 20 -A cos 20), from (44-4), j 

, from (44-5) J 
(44-6) 



Then from (44-4) and (44-6) we get Table I. A glance at this table 
shows that the slope of the graph is positive for 6 = 45 and = 90 
and that it is negative for 0=135-i# 0=135 and 0=180, 
whatever be the values of the constants involved. Hence the 
maximum value of M occurs between = 90 and = 135 J/?. 
Having this information and the graphs of (a), (6) and (c), we 
construct the graph for M as indicated by (5) in Fig. 42. 

TABLE 1 



8 


M 


dflf /cM 








TF^ulQ 




1 


j. 


45 


W 1 + JCO 2 cot a - p. (I 2 

V 2 


V 2 


90 


TFZ+CWcota 


J.02 


135~i/? 





TFZcos(135-^/?), negative 




1 


-, 


135 


y-M + iC^cota + I^O 2 


Wl CO, 2 cot a 




V 


V 2 


180 





-m-^n 2 
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Curve (3) is the graph of (44*3) and shows that in the interval 
6 = to Q = 75 (approximately) the gyroscopic moment is negative 
and actually tends to lift the muller from the pan. The total 
moment is also negative in the interval $=0 to 0=55 
(approximately). 

Note. There are crushing machines similar to that of Fig. 38 in 
which the pan rotates under the mullers, causing them to rotate on 



Afi 



(51 
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60 



90 120 

Fig. 42 



150 180 



their axles but not revolve about the vertical axis. In such machines 
there is no gyroscopic action whatever, the crushing being due 
entirely to the weight of the mullers. The crushing pressure 
is therefore much less than that in the gyroscopic type of 
machine. 

Numerical Example. Let us consider an edge mill having a single 
muller weighing 2200 pounds and with a mean radius of 1 .5 feet, and 
let us assume that the muller rotates about the vertical shaft at 
40 r.p.m. Then Q = (40 x 27r)/60=f7T radians per second. If 



then 



=l^ 2.5 feet, 



1 5 



and 



a = 30 58' = 31, say. 



112 THE GYROSCOPE 

Taking the radius of gyration of the muller about its axle to be 
0.8r = 0.8 x 1.5 = 1.2 feet, we have 



2200 
A = -55-5- (2.5) 2 + 1(98.4) approximately 

32.2! 

= 476 approximately. 



mu 4. a 98.4cot31 A Q _ 

Then tan/? = -5 = - - =0.345, 

jQ. 4:7 O 

^=19 02 7 = 19, say. 

We next find the value of 6 for maximum pressure moment M. 
Putting dM/dd^Q in the second of equations (44-6), we have 



On substituting the above numerical values for the constants in this 
equation, we get 



5500 cos - Y V{^ 

which reduces to 

cos (20 + 19) - 0.622 cos 6 = 0. 

This equation must be solved by trial, and it will be found that 
0= 117 14' satisfies it. 

Substituting into (44-5) this value of 9 and the numerical values 
of the constants as given above, we have 

M = 5500 sin 117 -8.78x503 sin 253 + 8.78x164 
= 4900 + 5660 = 10,560 lb.ft. 

The gyroscopic moment is thus greater than the moment due to the 
weight of the muller. 

For 6 = 90 we would have 

M = 5500 + 2880 = 8380 lb.ft. 

In this case the gyroscopic moment is only slightly more than half 
as great as before. But in the case of a grinder such as that pictured 



GYROSCOPIC ACTION IN ROTATING BODIES 113 

In Fig. 38, for example, the G.R.M. would still increase the pressure 
on the pan by more than 50 % above that due to the weight of the 
rollers. 

Equation (44-3) shows that the gyroscopic moment varies as the 
square of the speed of revolution about the vertical shaft. 

The rotational speed of the muller about its axle is found from 
equation (44-2) to be 



(b) The Griffin Grinding Mill. Although this mill grinds by intense 
pressure due entirely to gyroscopic action, it is quite different in 
construction from the edge mill. In the Griffin mill the grinding is 
done by a pendulous cylinder attached to a vertical shaft by means 
of a Hooke or universal joint. The rotating cylinder presses against 
the inner wall of a heavy steel pan as it rolls around inside the pan, 
as indicated in Fig. 43. This grinder is also called a pendulum mill 
because of the pendulous roller it employs. 

To derive the formula for the crushing moment in this mill, we 
again apply equation (27-1). In this grinder the rotational speed of 
the vertical shaft above the Hooke joint is the same as the speed (*> 
about the instantaneous axis, as can be seen by imagining the joint 
replaced by a flexible shaft. Furthermore, a moment's reflection 
or a simple experiment will show that when the cylinder rotates to 
the right (clockwise) and rolls around the inner wall of the pan 
without slipping, it will move around the central vertical axis in a 
counterclockwise direction as viewed from the fixed point (Hooke 's 
joint) above. In other words, the axle of the spinning roller pro- 
cesses to the left. This means that 9$ and O have opposite signs and 
that the spin oj of the vertical shaft is the algebraic sum of these 
speeds of opposite sign. Hence we have 

(*) = < + 1. (44-7) 

Since the speed of the vertical shaft is the controlling motion and 
known speed in this grinder, we must express <p and O in terms of 
o) before using (27-1). 

From equation (44-1) we have 
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Fig. 43 
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Adding 1 to each side of this equation, we have 

Q i 1 = sinos i i 
' 



from which 



<p sin (6 a) 

But O + ^ = ^, by (44-7). Hence 

a) __ sin a -f sin (6 a ) 
^ sin($ a) 



Now replacing <j> in (44-1) by its value as given by (44-9), we get 



Q-. *?* ,. (44-10) 

-- a ^ 



On replacing $ and <Q in (27-1) by their values given in (44-9) and 
(44*10), disregarding the minus sign, and then expanding the factor 
sin (6 a) in the numerator of the result and simplifying slightly, 

we get w 2 sin a sin # 

K = =-. - : -rfi - r=- 9 ( G co s a sin 6 A sin a c o s 6 ) ( 44 1 1 ) 
- 2V ; v J 



for the gyroscopic reaction moment. 

If the cylinder is rotating to the right (clockwise) as viewed from 
above, the spin vector must be drawn downward along the axis of 
the cylinder, and the precession vector must be drawn vertically 
upward. Then since the rule for the direction of the gyroscopic 
reaction moment requires that it tend to turn the positive end of 
the spin vector into the direction of the positive end of the preces- 
sion vector, we see that the moment K tends to press the cylinder 
against the wall of the pan. Since the gravity moment due to the 
weight of the cylinder tends to pull the cylinder away from the wall, 
it is plain that the gyroscopic moment must exceed the gravity 
moment in order to keep the cylinder in contact with the wall. The 
crushing moment of the cylinder is the algebraic sum of these two 
moments, so that we have 

- A sina cos 0), 



:r-: - r~7B - ^To 
[sin a -f sin (6 a)] 2 

(44-12) 

8-2 
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or 



(44-12a) 

This moment, in the case of a given cylinder, depends on o) and 6. 
In determining the sign and magnitude of 6 we must keep in 
mind the physical conditions of the problem: 

(1) <j> and O must be of opposite sign. 

(2) <j> > O, because the circumference of the pan is necessarily 
greater than that of the cylinder which rolls inside it. 

(3) A^>G in this problem. 

(4) The gravity moment Wl&ind and the gyroscopic moment K 
must be of opposite sign. 

Equations (44-8) and (44* 12 a) show that these conditions will be 
fulfilled only when 6 is a negative acute angle, or 90 < 6 < 0. 
There is no best value of 6 in this case as was true with the edge mill. 
The gyroscopic moment continues to increase as 6 decreases in 
magnitude. 

Numerical Example. The muller of a Griffin grinding mill is in the 
form of a frustum of a right circular cone and weighs 980 pounds. 
The altitude of the frustum is 1 foot, the radius of the upper base is 
1.077 feet, and the radius of the lower base is 0.5 foot (Fig. 43). 
(The generating line of the conical surface makes a 30 angle with 
the geometric axis of the frustum.) 

The center of gravity of the muller is found to be 0.383 foot from 
the upper base. Taking Z = 3.383 feet and the radius through the 
center of gravity as 0.856 foot, we have 

0.856 rt ^ 
tan* = = 0-253, 

Further computation gives 
0=11.52, 

980 
^ = 322 (3 - 383)2 + 7 ' 98 = 356 ' 

Assuming that the vertical shaft (the part above the Hooke 
joint) rotates at 240 r.p.m., we have 

240 
o) = -^r- x 27T = STT radians per second. 
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Taking Q =30 and substituting these data into (44*12) or 
(44- 12 a), we find 

Jf =-1658 + 30810 ==29, 150 lb.ft. 

Formulas (44-9) and (44- 10) give 

<j> = 38.75 radians per second = 370 r.p.m., 
Q= 13.65 radians per second = 130.3 r.p.m. 

To find the angular velocity w necessary to keep the muller in 
contact with the wall of the pan, we must have 



Hence o) > 5.83 radians per second = 55.7 r.p.m. 

45. Gyroscopic action on oblong projectiles fired from rifled guns 

A rifled gun barrel is one having spiral grooves on its interior 
surface or bore. Such grooves impart a rotary motion to projectiles 
fired from the gun. The pitch or twist of the grooves is such as to 
give the projectile a high rotational spin about its geometric axis. 

When an oblong projectile is fired from a gun having a right-hand 
twist, it swerves slightly to the right of the plane of fire (the vertical 
plane containing the axis of the gun barrel), as viewed from the rear 
of the gun; and when such a projectile is fired from a gun having 
left-hand twist, it swerves slightly to the left of the plane of fire. 
This swerving to the right or to the left of the plane of fire is called 
drift and is due mostly to gyroscopic action. 

Consider an oblong projectile fired from a rifled gun having 
right-hand twist. At the instant when the projectile leaves the gun, 
its axis coincides with the axis of the barrel and is therefore tangent 
to the trajectory (path of the center of gravity of the projectile) at 
that instant. Immediately after the projectile leaves the gun, 
however, its center of gravity is pulled downward continuously by 
the force of gravity, thus causing the center of gravity to trace a 
curved path in a vertical plane. Since the high rotational spin tends 
to keep the axis of spin oriented in its original direction in space, 
and gravity is all the while changing the direction of the tangent to 
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the path, of the e.g., it is plain that the axis of the projectile will 
soon make an angle a with the tangent, as indicated in Fig. 44. 

The resultant air resistance opposing the motion of the projectile 
is practically parallel to the tangent and always intersects the axis 
of spin at some point A in front of the e.g. of the projectile. In Kg. 44, 
designates the e.g. of the projectile, OT the tangent to the 
trajectory 5 and R the air resistance. 




Fig. 44 

The force R can be replaced by a force R acting at and a couple 
R, R as shown. The moment of the couple is Rasina, where 
a = OA, and tends to rotate the projectile about a transverse axis 
passing through 0, this axis being normal to the vertical plane 
containing OT. This moment is represented by the vector M and 
must be drawn to the right (as viewed from the rear) of the vertical 
plane through OT in order to conform to the right-handed-screw 
convention. 

Let & denote the angular velocity of the projectile about its 
geometric axis. Then the moment M will cause the point of the 
projectile to precess with angular velocity O about a line through 
and drawn in the direction of the resultant air resistance, which 
line is practically OT. By formula (22-3) the rate of precession is 

/ _ 



Co)smoc Cwsma Co)' 

where G denotes the moment of inertia of the projectile about its 
axis of spin. 

It might seem from the above considerations that the point of 
the projectile would describe a circular helix about the trajectory, 
but this is not the case. It can be shown both graphically and 
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analytically that the point of the projectile describes a eyeloidal 
space curve which always lies on the right-hand side of the vertical 
plane through O T. The proof of this fact lies properly in advanced 
works on exterior ballistics and will therefore not be given here. * 
The fact that the axis of the projectile is always pointing to the 
right of the vertical plane through the tangent shows that the air 
force (which is parallel to the tangent) has a small rectangular 
component which is perpendicular to the axis of the projectile and 
acts toward the right. It is this small force which causes the pro- 
jectile to swerve to the right. 

46. Gyroscopic effects on propeller-type airplanes 

An airplane propeller spinning at high speed has the dynamical 
properties of a spinning gyro. In the case of an airplane equipped 
with a single propeller a turn in azimuth, for example, is a forced 
precession of the spin axis about a vertical axis. If the propeller is 
rotating clockwise as viewed from the rear, this forced precession 
will induce a G.B.M. which will cause the nose of the plane to move 
upward or downward according as the turn in azimuth is to the 
left or to the right. Turning the nose upward rather suddenly will 
cause the plane to turn to the right, thus tending to follow an 
upward spiral. On the other hand, turning the nose downward 
rather suddenly will cause the plane to turn to the left, thus tending 
to follow a downward spiral. 

When an airplane is equipped with two propellers rotating in 
opposite directions and at the same speed, the gyroscopic effect of 
either propeller annuls that of the other, and there is therefore no 
appreciable gyroscopic effect on the plane as a whole. 

* The interested reader should consult the following: New Methods in 
Exterior Ballistics, ch. 6, by F. R. Moulton, 1926; Ballistics of the Future, 
ch. 8, by Kooy and Uytenbogaart, 1946; Aussere Ballistik, sections 57 and 
58, by Cranz, Von Eberhard, and Becker, 1925. 



CHAPTER VI 

The Gyroscope as a Direction Indicator 
and Steering Device 

A. THE GYROSCOPIC COMPASS 

47. The principle of action 

The stubborn tendency of a fast-spinning flywheel to maintain its 
plane of rotation early suggested the use of a gyroscope as a compass 
to replace the magnetic compass. If the axle of a fast-spinning free 
gyroscope is pointed toward a fixed star, for example, it will 
continue to point toward that star (because of the principle of the 
conservation of angular momentum) as long as the gyro is kept 
spinning at high speed. But if the spinning gyro is set with its axle 
horizontal and in the meridian plane at any point on the earth's 
surface, the axle will remain horizontal and in the meridian plane 
if it is compelled to precess about the earth's axis at the same 
angular velocity as the earth rotates, or precess about the vertical 
with an angular velocity equal to the vertical component of the 
earth's angular velocity at that place, thus making the change in 
direction of the gyro axle keep pace with the change in direction 
(with respect to fixed space) of the meridian at that point (see 
Fig. 45). This is the principle of action of the gyrocompass. Hence, 
if the axle of a spinning gyro is pointed toward the north point 
of the horizon, it will continue to point to the true north if it 
is compelled to precess in the right direction and at the proper 
rate. 

The gyro axle can be made to precess at any desired rate by 
attaching a pendulous weight to the frame which supports the gyro 
axle, as we learned in Chapter IV. In order to determine on which 
side of the frame the weight must be attached, it is first necessary to 
know in what direction the axle must precess. Gyrocompasses are 
mounted on board ship with the gyro axle practically horizontal and 
the compass card (dial) facing upward. Hence the compass must 
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be read by looking down on It from above. A glance at Fig. 45 
shows that the spin axis must precess to the left (counterclockwise) 
as viewed from above, in order to keep pace with the changing 
direction of the meridian. Consequently, the pendulous weight 
must be attached as indicated in Fig. 46, for in that position the 
moment vector of the weight is directed as indicated in the figure 
and the spin vector must be rotated counterclockwise (as viewed 
from above) in order to bring it into coincidence with the moment 
vector (Art. 23). 

The precession of the gyro axle may be looked upon either as a 
precession of angular velocity Q about the polar axis of the earth 
or of an angular velocity Q, sin A about a vertical axis at latitude A 
(see Fig. 45). In either case we have, from Fig. 46 and (22-3) in the 
first case or (22-4) in the second case, the relation 



(47-1) 

where co = angular velocity of spin of gyro in radians per second, 
O = angular velocity of earth's rotation in radians per second, 

C = moment of inertia of gyro about its axis of spin in pound- 
feet-second units, 

A = latitude of position of gyrocompass, 
w = weight of pendulous weight in pounds, 
a = moment arm of w in feet, 
ft = angle which arm a makes with vertical. 

If /? is the value of ft which is necessary to give the gyro axle the 
proper precession to keep it in the meridian, we have from (47-1) 

. n O^QsinA 

sm/? = - . 

u wa 

But y? , as will be seen later, is a very small angle of & to 10' in 
magnitude. Hence we may replace sin/? by /? . Then we have 

v ' 



wa 



The gyro axle is therefore not horizontal in its equilibrium position 
but is tilted slightly upward. 
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The angular velocity of the earth's rotation, in radians per second, 
is equal to 2/r divided by the number of seconds in a sidereal day, 

or 27T 

Q=r7^ ,. A - = 0.000072921158 radian per second. 
86164.091 

This number, rounded to four significant figures, will be used 
frequently in the work that follows. Note that 
O 2 = 0.000,000, 005,317,3. 

This quantity is so small that it will henceforth be neglected wherever 
it occurs, except when it is to be multiplied by a very large number. 

48. Differential equations of motion of a disturbed gyrocompass axle 

To determine the motion of the gyro axle when disturbed from its 
position of equilibrium, we take a right-handed set of rectangular 
axes rigidly attached to the frame of the gyroscope and with origin 
at the center of gravity of the rotor (which is practically the center 
of gravity of the movable system), as indicated in Kg. 46. The 
#-axis coincides with the axis of spin, with positive end toward the 
north. The y-axis coincides with the axis of nodes (moment axis of 
pendulous weight), the positive end being directed toward the 
reader. The 2-axis is perpendicular to the plane of the other two 
and is positive downward, because the compass card is graduated 
with at the north and with the numbers increasing toward the 
east (clockwise as viewed from above) to 360 and is read by looking 
down on it. Horizontal or azimuth angles will be taken as positive 
from the north toward the east. Since the earth rotates from west 
to east (as one looks toward the north), O will always be taken as 
positive. All the above conventions agree with the right-handed- 
screw rule used throughout this book. 

The mounting of the gyro is such that its axis of spin is free to 
move in any manner, but the axis of nodes (OF) is constrained to 
remain approximately horizontal. 

Figure 46 shows that the components of the angular velocity of 
the earth and of the frame along the coordinate axes are as follows: 

c = D cos (A /?) cos a, 

W - 

(48-1) 
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Hence the components of the angular momentum on the are 
H x = CO cos (A /?) cos a + Cw, . 

H y =A--AQGO&(\--/})Bwa, 

(M> " \^!to m ij 

j-j A doc 

z dt * 

where A denotes the moment of inertia of the gyro and frame about 
an axis through the center of gravity of the gyro and perpendicular 
to the axis of spin. Hence we have 



~ = 
cLt 



cosasin(A /?)-~--smacos(A-/?)-^ L 

dt dt_j 



dH 



It will be seen later that dajdt and dfijdt are very small. Hence the 
products l(da/dt) and l(dj3ldt) are negligible. Neglecting such 
products in the above derivatives, we get 



dt ' dt dt*' dt dt*' 

On substituting into (16-4) these derivatives and the quantities in 
(48- 1) and (48-2), then neglecting all terms containing the negligible 
quantities Q 2 , Q(da/<ft), &(dpldt), and (dajdt) (dfi/dt), we get 



dt* dt 



or, after expanding sin (A/?) and cos (A ft): 

A ^ + Ow ^ - <7wO sin A cos ft + (7wQ cos A sin ft = - wa sin ft, 



A -T - Co) -~ + Ca)l sin a cos A cos B + CuQ. sin a sin A sin /? = 0. 
dt 
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Now since /? is a very small angle, we may neglect the terms 
containing the product Qsin/?. We may also put cosy?=l and 
sin /?=/?. Then the above equations become 



(48-4) 

10 j n i ^ ' 

^rz C f w-~ 
dt* dt 

These are the fundamental differential equations which determine 
the motion of the gyro axle when disturbed from its position of 
equilibrium. Since the quantities d 2 oc/dt 2 , dccjdt, d 2 fijdt* and dffldt 
are all zero in the equilibrium position, we can find that position 
by putting these quantities equal to zero in equations (48-4). We 
thus get 



Hence a=0 

, ~ 

and p n 



wa 



This value of /? is the same as already found in (47-2). Hence in the 
equilibrium position the gyro axle is in the meridian plane and its 
north end elevated at the small angle /? above the horizon. 

Equations (48-4) can be reduced to a symmetrical form by 
making the substitution 

(48-5) 



wa 

a n Co) 

or p 



wa 



Making these substitutions into (48-4) and putting sin a = a, since a 
is a small angle, we get 



(48 ' 6) 



doc 
(a) A- + C< t >- 



(6) ^-<7w 
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49. Derivation of the differential equations by the GJtM. 

Equations (48-4) can also be derived by considering all the 
moments, external and G.B.M., about the axes OF and OZ. 

Let us first consider the moments about OF. The moment of 
the pendulous weight is wasm/}, negative because it tends 
to decrease ft. The precession Qsin(A~ ft) about OZ induces 
ao.B.M. -[C f w(-asin(A-^))] J orC f ^asin(A~-/?),aboutOr. Also, 
the angular velocity da/dt induces a G.B.M. -Ca)da/dt about OF. 
Hence the total moment about OF is 

wa sin ft + Ca)l sin ( A ft) Co) -j- . 

(Jut 

Then by the fundamental equation of rotation (20-2) we have 
A -r~ = CcoQ, sin (A /?) Cw-j- wa sin/?, 

Cvt Cut 

or A -- + Go)-y-~ Co)l sin (A-/?) + wasin/? = 0, 



which is the first of equations (48- 3). 

Now considering the moments about OZ, we note first of all that 
there is no external moment. The forced rotation (precession) 
dfijdt induces a G.E.M. G^d^jdt) about OZ. Also, the precession 
Q, cos (A ft) sin a about Y induces a G.B.M. 

Co)[ O cos (A /?) sin a] 
about OZ. Hence by the rotational equation (20-2) we have 



A -j = Co)---C(*)l cos (A ft) sin a, 
ctt ctt 



or -7~ ? ~--~ 

at* at 

which is the second of equations (48-3). 

Although the G.E.M. method is much shorter than the use of 
equations (16-4), it is less foolproof. The student is advised to write 
down the angular velocities o) x , a) y) o) s on the three axes and use 
equations (16-4) until he becomes expert in the use of the G.R.M. 
method. 
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56. Solution of the equations 

The simultaneous equations (48-6) can be solved by eliminating 
one of the variables and thereby obtaining an equation of the fourth 
order in the other variable. The elimination can be effected in 
several ways. A simple and direct method of eliminating x> fo*" 
example, is the following: 

Differentiate (a) with respect to t and then substitute into the 
result the value of dx/dt obtained from (6), thus obtaining 

/ v A d*x ~ d 2 a Awa d*a ^ , 

(c) 4-j + O f <>-= s + - 7r -- -j~z + lwaQ08A.a = Q. 

v ' dt* dt* Ooj dt* 

Now differentiate (6) twice with respect to t, obtaining 

*-**.-* 

Eliminating d 3 xldfi between (c) and (d) and then simplifying the 
result slightly, we get 



(50-1) 

By following a similar procedure we can eliminate a from 
equations (48-6) and get 



* TOT + (V^* + -AC&Q cos A + Awa) - + Cu&wa cos A . x = 0. 



(50-2) 

This equation is exactly the same as (50-1) except that a has been 
replaced by x- Hence the motion of the gyro axle in the vertical 
direction is of exactly the same character as its motion in the 
horizontal direction (in or parallel to the plane /?=/?<)) 
The auxiliary equation for either (50-1) or (50-2) is 

A V 4 + ((7V 4- ACa)l cos A + Awa) r 2 + Ca)lwa cos A = 0. 
Hence by the quadratic formula we have 
- (<7V +AC(*)l cos A +Awa) 

2 = V{(<7 2 ^ 2 + ACwCl cos A + Awa) 2 - 4=A*Ca)lwa cos A} 
T ~~ 2A* 

Now since (C^+AGcoQ.Go&X+Awa) 2 is very much greater than 
, the radicand is positive. Consequently, the two 
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values off 2 are real. Moreover, they are both negative, because the 
radical is less than C*a) 2 + ACa)lco8& + Awa. 

Let these negative values of r 2 be k\ and jfcf , so that 

r 2 i* 2 r 2 2 

*l ^i? '2 K> 2 t " 

Then r i= *&i r* = ^2- 

Hence the general solution for a is 

a = c x sin k^ -f c 2 cos ^2 + c 3 sin & 2 -f- c 4 cos Jc^t. 

A similar solution holds for %. 

The above solution represents two simultaneous and superposed 
simple harmonic motions of periods STT/^ and 2njk 2) respectively. 

51. A numerical example 

At this point it is instructive to work out a numerical example in 
order to obtain an idea of the magnitudes of the several terms in 
(50-1). The following data are approximately correct for some of 
the later Sperry gyrocompasses: 

a) = 8500 r.p.m. = 890 radians per second, 
(7=0.25 in foot-pound-second units, 
A = 0.15 in foot-pound-second units, 
wa = 6.5 pound feet. 
Hence C f 6> 



Also 1 = 0*00007292 radian per second. 

Take A = 41 N. 

On substituting these quantities into (501) we have 



0.0225 ~ + (49506.25 + 0.0018370 + 0.975) + 0.079592 a = 0, 
(MI 



. , 

or 0-0225 7 + 49507.226837 -^ -f 0. 079592 a = 0. 

at* at* 

Hence 0.0225 r 4 -f- 49507.226837 r 2 + 0.079592 = 0, 

and therefore 

2 __ -49507.226837 +V{(49507.226837) 2 - 0.007 16328} 
r _ 0.045 ' ^ J 

2 ^ -49507.226837 ( _ F 0.00716328 1 

ri "~ 0.045 { L (49507.226837) 2 J 
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Expanding the quantity in square brackets into a binomial series 
by means of the formula 



... 

and retaining only the first two terms of the expansion, we get 
- 49507. 



y.2^. 

Hence 



507.226837 f F 0.00358164 "H 

OM5 II (49507.226837) 2 ] j * 



2 _ -49507.226837 I 0. 

fl ~" "" + 



00358 164 



O045 1 "" (49507.226837J2j 

-0.00358164 



0.045 x 49507-226837 
.*. r lS= 0.00126795*. 
The period of oscillation is therefore 

27T 



= 82-6 minutes. 
For the other two roots we have 



49507-226837 f 0.00358164 } 

-0.045 | + (49507.226837) 2 J 

99014.453674 0.00358164 



0.045 0-045 x 49507.226837 

= -2200321.192754. 

Hence r 2 = 1483.348i, 

and the period in this case is 



This exceedingly short period indicates that this type of motion 
consists of an extremely rapid tremor of the gyro axle. It is of no 
importance practically and is soon wiped out by friction. 
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52. Simpler forms of the differential equations solution 

The magnitudes of some of the terms In the above numerical 
example are seen to be negligible in comparison with others. The 
coefficient of d^a/di*, for example, is very small in comparison with 
that of d 2 a/dt*. Furthermore, the terms ACo)Q,QOB A and Awa are 
exceedingly small in comparison with (7 2 &> 2 . These facts suggest that 
the small terms just mentioned may be neglected without appreci- 
able error. On neglecting these terms in (50-1) and (50-2), and then 
dividing the resulting equations throughout by Co>, we get 



==Q, (52-1) 

+ Qwa cos A x = 0. (52-2) 



Here the auxiliary equation in each case is 
Cwr 2 4- lwa cos A 0. 
!lwa cos A 



The period is therefore 

(52-3) 



. ,. 

Llwa cos A 

Now using the numerical data of Art. 51, we have 

222 ' 5 



0.00007292 x 6.5 cos 41 
= 4955.4 seconds 
= 82.6 minutes, 

which is the same value as found for the long period in the previous 
case. The dropping of the negligible quantities in (50-1) thus caused 
only the loss of the rapid tremor oscillations, which is of no import- 
ance. We shall therefore use equations (52-1) and (52-2) instead of 
(50-1) and (50-2). The period of oscillation both vertically and 
horizontally is given by (52-3). 

We shall now solve for a and i/r with assumed initial conditions. 

From (52-1) we have n 9 t r\ in 

v J CW 2 + OwacosA = 0. 

9-2 
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Hence . /OwacosA . T 

r -. -j- t / . = + ^ say. 

>y G<y 
Then a = C^ sin kt + C<> cos kt, 

and = C^fe cos kt 2 A sin M. 

Now assuming that a = a and da/d = when t = 0, we have 



Hence C^ = 0, and 

(52-4) 



which represents simple harmotion with amplitude a . 

We wish to find x without further integration, in order to avoid 
bringing in additional constants of integration. Consequently we 
go back to equations (48-6). From (52-4) we have: 

-=-= a Q k&mkt y - = a & 2 cos&, -5-^ 

Clt (Li (M> 

Now differentiating (48-6) (6) with respect to time, we have 



. 
at 

Eliminating d 2 xldt* between this equation and (48-6) (a), we get 



Replacing d^ajdt 5 and dajdt by their values as found above and 
solving the result for %, we have 

y= |C^ + ^QcosA-- 7r -~)sinR 
A wa \ Co) ] 

When k is replaced by its value / - -^ , the value of x becomes 

a //i^acosA\ ( Af . . -4 2 O^acosA\ . 7 
% = - - / - TC - |(7w + -4QcosA -- 7W ~^ - Ism &. 
A wa V V #6) / \ 
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Neglecting the second and third terms in the parentheses because 
of their smallness in comparison with Co> and then reducing the 
result slightly, we get 

(62 ' 6) 

This equation (52-5) shows that ^=0, or /?=/? , when * = 0. Hence 
when the gyro axle reaches its extreme position to the right, it is then 
in the plane /} =/? . The above equation (52-5) also shows that the 
spin axis executes simple harmonic motion in the vertical direction. 

Eliminating t by solving (52-4) and (52-5) for cosJU and sinB, 
respectively, and then squaring and adding the resulting equations, 

we get 2 

<* 2 X 2 



which is the equation of an ellipse in a plane perpendicular to the 
spin axis in its settled position. The north end of the spin axle 
therefore traces an ellipse as it executes simple harmonic motion in 
the horizontal and vertical directions. 

53. Comparison of the amplitudes of the horizontal and vertical 
oscillations. Locus of the end of the gyro axle 



o- *u J.M. /CtoQcosA . 

Since the quantity / - is very small in comparison 

^ J *j wa J r 

with unity, the amplitude of % is very much less than a , the 
amplitude of a. For the data of the numerical example worked in 
Art. 51 we have 

(7o>O cos A /222.5 x 0.00007292 cos 41 



wa A/ 6.5 

= 0. 0434 = ^g-, approximately. 

Hence the amplitude of x is only ^ that of a. 

Equations (52*4) and (52*5) show that a and x are out of phase by 
a quarter of a period; for when Jet = 0, a = a and # = 0, whereas when 
& = JTT, a = and ^ is a maximum. This means that the north end 
of the gyro axle describes a very elongated ellipse in a plane 
perpendicular to the axle, the major axis of the ellipse lying in the 
plane /?=/? . 
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To find ft for the data of the example In Art. 51, we have, from 

equation (47-2), 822.6x0.00007292 sin 41 

A)= _ 

= 0.001 6376 radian 
= 5'.6. 

54. A third method of deriving the differential equations of motion 
If a = 2, for example, the maximum value of x will be 

120' 



23 



= 5'.2. 



Since the period of oscillation of gyrocompasses is about 84 minutes, 
the time required for the axle to move through 5'.2 will be 
21 minutes. The average angular velocity of the upward motion for 
a quarter period is thus 

5'.2 0.00151 radian AAAAAA10 ,. , 

^ = r~ = "~n^ - T- = 0.0000012 radian per second. 
21 minutes 1260 seconds 

Since this angular velocity is so very small, the upward angular 
acceleration d^x/dt* is also exceedingly small and is therefore 
negligible in equation (48*6) (a) in comparison with the other terms 
in that equation. This fact gives us a short method of finding the 
final differential equation of a disturbed motion and the period of 
oscillation. 

On putting d 2 #/W = in (48-6) (a), we get 

Co)^+wax = . (54-1) 

at 

Hence by differentiation, 

d 2 a 

Co) d*ot 
. 

wa dt* " 

Substituting this value of dx/dt into (48-6) (6), we get 



or 

dt wa 



- 1 - 5 - 
wa ) dt 2 
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Now A is totally negligible in comparison with 2 o> 2 /i0a. Hence we 

neglect the A in the parentheses and then divide the resulting 
equation throughout by 7<y, thereby obtaining 

dP& 

Co} ~df* + ' Wa COS ^ " a = J (54-2) 

which is (52-1). After solving this equation for a, we could find 
doc/dt and then substitute it into (54-1) to find #. 

55. Motion of the gyroscopic compass with damping 

In the preceding pages of this chapter we have considered the 
motion of an undamped gyrocompass. Such a compass, when once 
disturbed, would continue to oscillate with undiminished amplitude 
if not slowed down by frictional forces. For practical reasons it is 
important that the oscillations be damped out as quickly as 
possible. Several methods of damping have been devised and used. 
The simplest and one of the best is the method devised by Elmer A. 
Sperry, which we shall now explain. 

In Arts. 47 and 48 we assumed that the pendulous weight w 
acted in the XZ-plane. In order to produce the precession necessary 
to keep the spin axis in the meridian and at the same time damp out 
the oscillations that might occur, Sperry devised an eccentrically 
weighted bail (semicircular ring) and fastened it rigidly to the inner 
ring of the suspension as indicated in Fig. 47. The damping was 
accomplished by placing the weight w on the bail a short distance to 
the east of the XZ -plane. In that position the weight will exert a 
small moment about OZ , the effect of which will be a slight down- 
ward precession of the spin axis. This precession will tend to 
decrease the angle of inclination ft and at the same time auto- 
matically decrease the azimuth angle a, because of the fact that the 
ratio of the semiaxes of the ellipse described by the north end of the 
spin axis is constant, this ratio being <J(Co)l cos A/wa) (Art. 52). The 
path of the north end of the spin axis with damping is therefore no 
longer an ellipse but is a spiral whose radius vector continually 
decreases to zero. 

To find the equations of motion of the spin axis under damping, 
let S denote the small angle (about 1) which the off-center arc of 
the bail subtends at (Fig. 47). Then the distance of the weight w 
from the z-axis is practically aS. If w be resolved into rectangular 
components wcosy? and wsin/?, the latter is directed backward 
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toward the south and thus has a clockwise moment about OZ as 
viewed from above and is therefore positive. The moment of 
weosft about either OF or OZ is zero, but the moments of w sin ft 
about these axes are 

M y = - w sin ft x a cos 8 = - waft, practically, 
= waftS, practically. 



On substituting these values into the right-hand members of 
equations (48-3), which came from (164) ? we get 




~*s -& 

Now expanding sin (A/?) and cos (A ft), neglecting the terms 
containing the product O sin ft, putting cos ft = 1, and sin ft=ft, we 

get da 



A -3-5- - G(J) -j: + (7o)Q sin a cos A - waftS- 0. 1 
c& 2 ^ ^ J 

To find the equilibrium position of the spin axis we put 

d a da dp dp 

dt dt dt dt 

in (55-2), thus obtaining the equations 

Co)l sin A = waft Q , 
Co)l sin a cos A = waft Q S. 
From the first of these we get 



_ 

ro wa v ) 

Substituting this value into the second equation we get 



But since S is very small, sin a is also small and may be replaced by 

a . Then we have ** i /~~ ^ 

(55-4) 
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The spin axis therefore points slightly to the east of the meridian 
when it comes to rest, but this small deviation is easily corrected. 
If we put A = in (55*3) and (55-4), we get a = and /? =0> 
which shows that at the equator the spin axis is horizontal and 
points true north. 




Fig. 47 

Since ft changes very slowly, eZ 2 /?/c& 2 is negligible in comparison 
with the other terms in the first of equations (55-2). Also, since a 
is a small angle, we may replace sin a by a in the second equation. 
Then equations (55-2) become 



__ Ca>i sin A + waft = 0, 
dt 



(55-6) 



= 0. 
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From the first of these we get 

doc . \ 
"Si sin j . 

Tt^ d(t Co) d 2 a 

Inen -=-= 375-. 

dt wa at* 

Substituting these into the second equation, we get 



_ + QMS + C(j)Q,a cos A = 



wa / dt* ' dt 

Now since A is negligible in comparison with C^jwa, we neglect 
it in the above equation and then divide throughout by Co), 
thereby obtaining 

_^ + _? ^L^L ? -a = . (55*6) 

Solving this equation, we have 

n waS 



-n~ 7? 
Co) Co) 



. wad /l/waS\ 2 waQcosA] 

whence r=~^ 7r ~ / hrFT" " 7^ } 

2Co) V I V 2 ^/ ^ w / 

If the spin axis is to oscillate, the radicand must be negative. 
Hence we have 

waS . /(waiicosA IwaS}^ 



and therefore the complementary function is 

oc c = exp ^ t (Ci sin kt + C 2 cos Jet), 
I 2t\^0) I 

where Jc = 



Since the right-hand member of (55*6) is constant, we can find 
a particular solution by writing 

a = N, a constant. 

Substituting this and its derivatives into (55*6), we get 

N = S tan A. 
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The complete solution of (55-6) is therefore 

a = exp - -^r-tf (Cj sin Id + C 2 cos kt) -f tan A. (55-7) 
The period of a complete oscillation is 



k llwaQ, cos A /waS \ 2 ] 
and increases with S. When 



flwaQ, cos A /wad \ 2 | 
VI 05 \2Ca)l 



Cu 
or ^=2 



the spin axis becomes aperiodic or dead beat. A value of S as small 
as 7 is sufficient to make the Sperry gyrocompass dead beat. 

The period of a damped gyrocompass is about 90 minutes. 

To find the damping factor of the oscillations, let T denote the 
period of a complete oscillation .Then if we count time from the instant 
t = , the time at the end of the first period is t = -f T. Hence the 
amplitude of the oscillation at the end of the first period is, by (55*7), 



= exp I - 1^ T\ exp I - 1^ t \ (C^ sin M Q + 2 cos kt Q ) + S tan A. 
The damping factor for a complete oscillation is therefore 



and for a half period (the time of swing from extreme right to 
extreme left) it is exp -~^- T . This latter is called the damping 
factor and is denoted by /. We thus have 



The value of /is usually about |. About 3 hours is required for the 
disturbed spin axis to settle to its equilibrium position. 
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56. The speed and course error 

In tlie preceding discussion of the gyroscopic compass we have 
assumed that it was installed ashore on a stationary base. A com- 
pass installed on a ship is subject i 
to certain errors due to the motion I N 
of the ship. For example, when j 
a ship is sailing due north at a I 
constant speed of v feet per second, I 
it has an angular velocity of vjR 
radians per second about an axis I 

through the center of the earth W V I R I 

and perpendicular to the meridian 
plane, R being the radius of the 
earth in feet. This angular velocity 
is represented by a vector directed 
to the west and therefore along 
the negative direction of the ?/-axis 
(Fig. 46). It could therefore be in- 
serted as a negative term in o> y in 
equations (48-1) and its deviating 
effect found from the solution of 
the amended equations.The devia- 
tion of the gyro axis can be found 
more simply, however, by the 
following method: 

In latitude A the horizontal pro- 
jection of the vector representing 
the earth's angular velocity is 
Q cos A, directed true north. The 
gyro axis points in this direction 
when the ship is at rest. When the 
ship is sailing north with speed v, Fig. 43 

however, the gyro axis points in 

the direction of the resultant of the combined vectors v/B and 
XcosA, as seen from the vector triangle Fig. 48. The westward 
deviation of the gyro axis from the meridian is therefore given by 
the equation 

v/E v 

T-Qyi f\i ~ ' ,. ____________^__, 

OcosA -BQcosA* 
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If the course of the ship is not true north but an angle 6 

with the meridian, the northward component of the speed is 
vGOsd. Hence in this case the deviation angle is given by the 
equation 



___. 

Since y is a small angle not exceeding 2 or 3, we may replace 
tanv by y and thus get 

* 6 



7 -BOcosA* 

If the course of the ship is south, or in a southerly direction, the 
compass deviation will be to the east of the meridian, as can be 
seen from a vectorial triangle. 

The speed and course error is easily corrected by mechanical 
means. The Sperry gyrocompass is equipped with a corrector which, 
when set by hand for the ship's speed and latitude, makes the 
required correction automatically, thereby causing the compass to 
indicate the true course of the ship. 

57. The Sperry mercury ballistic gyroscopic compass 

The bail type Sperry gyrocompass is subject to error on a rolling 
ship when the ship is proceeding in an intercardinal direction. The 
error is caused by accelerations produced by the rolling. To prevent 
this type of error the weighted bail has been replaced in the more 
recent Sperry compasses by a pair of U -tubes with cylindrical cups 
at the ends and partially filled with mercury, as indicated in 
exaggerated size in Fig. 49. The planes of the IPs are parallel to the 
spin axis of the gyro, and the U -tubes are rigidly fastened to the 
gyro frame so as to tilt with the gyro axle. 

It is plain from Fig. 45 that when a free gyro is spinning at high 
speed the north end of its axle will rise as the earth turns eastward. 
This upward tilting of the axle and the U -tubes with it causes 
mercury to flow southward, increasing the weight of the south cups 
at the expense of the north cups. This unequal weight of the cups 
produces a moment about the axis of nodes and thus causes preces- 
sion of the gyro spin axis. If the rotation of the gyro is counter- 
clockwise or westward, the precession of the spin axis will also be 
westward. When the size of the cups and the quantity of mercury 
are so adjusted that the westward precession is just sufficient to 
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keep the gyro axis in the meridian, the compass will indicate true 
north. The reader will note that whereas the gyro rotates to the 
eastward in the bail type of compass it must rotate to the westward 
in the mercury type. 

The inertia and viscosity of the mercury prevent the flow from 
being instantaneous from the north cups to the south cups, and this 




Fig. 49 

means that there is a time lag between the tilting and the moment 
it produces. This time lag and the decreased amplitude of the 
mercury surging are the factors that prevent intercardinal rolling 
error. The mercury ballistic gyrocompass is comparatively in- 
sensitive to quick alternating accelerations such as are produced 
by a rolling ship. 

The damping of this type of gyrocompass is accomplished by an 
eccentric arrangement that produces a small moment about the 
2-axis, just as in the bail type. 

The mathematical theory of the mercury ballistic Sperry compass 
is very similar to that of the bail type andleads to thesame results.The 
interested reader will find this compass fully described and treated 
in pp. 88-100 of the The Theory of the Gyroscopic Compass (1944), by 
A. L. Rawlings, one of the inventors of the mercury ballistic type. 
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B. GYROSCOPIC STEERING 
58* Rate gyroscopes 

A rate gyroscope Is a gyroscope constrained to one degree of 
freedom. Although such a gyroscope has only one degree of 
complete freedom (the freedom to spin about Its axis of rotation). 




Fig. 50 

a restrained and limited rotation is permitted about its torque axis, 
or axis of nodes. Such rotation is limited by springs which can exert 
a moment about the torque axis. 

Fig. 50 is a schematic representation of the mounting of a rate 
gyroscope. The outer gimbal frame is rigidly fastened to a platform, 
which may be the deck of a ship or the floor of an airplane. If an 
airplane, for example, is flying horizontally and makes a turn to the 
right or left, the platform turns with it and thus rotates about a 
vertical axis through the center of gravity of the gyro. Such a 
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rotation is a forced precession of the spin axis and therefore induces 
a gyroscopic reaction moment about the torque axis. The rotation 
caused by this induced Q.R.M. is resisted by the springs shown in the 
figure. 
Let: &= spring constant, 

a = distance from torque axis to point of attachment of 
spring, 

F = force exerted by spring, 

^=tilt angle. 

Then a<j) = stretch of spring (approximately), 

and F=ka<j>, by Hooke's law. 

The moment of F about the torque axis is therefore 



<> J say. 

But this moment is equal to the G.R.M. = OOJI/T Co) (di/r/dt), where \Jr 
denotes the precession angle. We therefore have 

&# = Cfc^. (58-1) 

Civ 

This equation shows that the tilt angle is directly proportional to the 
rate of turn of the platform. When the tilt angle is shown on a dial, 
the aviator can therefore see at a glance the magnitude and 
direction of the turn his plane is making. 

Since the rotation of the platform will usually be at a variable 
rate, we find from (58-1) ^ ^^ 

tti-TW (58 ' 2) 

The rate of change of the tilt angle therefore indicates the accelera- 
tion of the platform rotation. 

Rate gyroscopes have many aeronautical and military applica- 
tions. They are the controlling elements in the following devices: 
the Turn Indicator used in blind flying, the Rate Gyro Gunsight 
that automatically computes the 'lead' required for hitting a 
moving target, the Gun-Bomb-Rocket sights used on jet fighter 
planes, the control gyro for the Sperry Ship Stabilizer, the Yaw 
Damper for stabilizing yawing oscillations in jet aircraft with 
swept-back wings, and the Hermetic Integrating Gyro (HIG) used 
in platform stabilization and in rate measurements. 
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59. Motion of a damped rate gyroscope 

Imperfections in a rate gyroscope may cause tremulous vibration 
or quiver of the torque or output axis. To reduce such disturbance 
as much as possible, rate gyroscopes are equipped with a damping 
mechanism (dashpot or otherwise). To find the motion of a damped 
rate gyroscope, let 

c = damping coefficient, 

/ = moment of inertia of gimbal system (gimbal, rotor and axle, 
and crossbar as shown in Fig. 47) about torque axis. 

Then the G.R.M. Cwijr induced by forced precession must maintain 
the angle ^, produce an angular acceleration of the gimbal assembly 
about the torque axis, and overcome the damping moment about 
that axis. Hence by equation (20.2) the differential equation of the 
oscillations about the torque axis is 






Solving this equation by the usual procedure, we have 
^2 _|_ *. _j_ _... Q^ 

whence c /// c \ 2 ^1 

r = . = -2/y {j- 12J 

Hence the complementary part of the solution is 
(j) Q = e-^/ 21 ((?! cos m^ -h C 2 sin m^) ? 
where 



The precessional velocity ijr may here be treated as a constant, 
since a small disturbance about the torque axis could have little 
effect on the precession of the base to which the gyroscope is 
attached. Hence to find a particular solution of (59-1) we put 

< = A, a constant, 



Then 

at at* 
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Substituting these quantities into (59*1), we find 

bA 



I " I ' 

whence A = -=- ^. 

The complete solution of (59*1) is therefore 

7 X cos mt + Cj sin m). (59-2) 



The disturbance thus soon dies out, leaving the necessary relation 



The period of the oscillation is 

*-- 

m 



2/ y 

practically, since (c/2I) 2 is negligible in comparison with bjl. 

To find the damping factor we note that the damping coefficient 
after a complete period is 



Hence the damping factor for a complete period is e~ p/21 and for a 
half-period (from one extreme of the swing to the other) is e~ CPIU . 
When P is replaced by its value from (59-3), the exponent cP/4J 
becomes 



cn 



and we then have for the damping factor 



60* The rate gyroscope as a differentiator and integrator 

The gyroscope, as is the case with all machines, is an instrument 
which receives an input and delivers a resulting output. If the input 
is a moment or torque, the output is a precession; whereas if the 
input is a forced precession, the output is a reaction moment. 
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The equation ^ 

W = cA 

Uv 

of Art. 58 shows that If the input of a rate gyroscope is an angle <j> 
(due to an applied moment), the resulting output is an angular 
velocity di/r/dt which is proportional to the input angle; and that if 
the input is an angular velocity d^/r/dt (due to a forced precession), 
the output is an angle (produced by the gyroscopic reaction 
moment) which is proportional to the input angular velocity. These 
operations and results correspond to the differentiation of an angle 
in the first case and to the integration of an angular velocity in the 
second case. The rate gyroscope can thus be used to compute 
angular velocities from given angles and to compute angles from 
given angular velocities. 

These properties of the rate gyroscope are utilized in automatic 
gun sights and other computing devices where angles and angular 
velocities are involved. The constants in the equation can be 
adjusted to obtain the desired angle or angular velocity. 

The equation 



shows in like manner that the rate gyroscope can compute an 
angular acceleration from a given angular velocity and compute an 
angular velocity from a given angular acceleration. 

61. Steering of torpedoes 

One of the earlier applications of the gyroscope was to the 
steering of torpedoes. A torpedo is equipped with two pairs of 
rudders, a horizontal pair for controlling the depth of the torpedo 
as it moves through the water to its target and a vertical pair for 
keeping the torpedo on its course to the target. Both pairs are 
controlled by a single gyroscope. 

The gyroscope carried by a torpedo is mounted in gimbals 
(Cardan suspension) with its spin axis parallel to the longitudinal 
axis of the torpedo, the outer ring of the gimbals being free to 
rotate about a vertical axis in a frame rigidly attached to the 
torpedo. The gyro is small, only 3 or 4 inches in diameter, but it 
spins at 1 4, 000-2 0,000 r .p .m. At the instant when the torpedo is fired, 
or a few seconds before, the gyro is started spinning and reaches 
its required speed within a few seconds. If the torpedo encounters 
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any force wMch deflects it horizontally, the spin axis of the gyro 
maintains its direction unchanged, but the rotation of the torpedo 
about the vertical axis of the gyro frame actuates the vertical 
rudders- and causes them to deflect the torpedo in the opposite 
direction.. As soon as the torpedo passes the direction of the spin 
axis of the gyro in this opposite swing, the rudder action is reversed 
and deflects the torpedo back in its original direction of swing. This 
back-and-forth deflection continues throughout the run of the 
torpedo, with the result that the torpedo travels in a serpentine 
path to its target. 

The depth of the torpedo is controlled in a similar manner by 
horizontal rudders which are actuated when the torpedo noses 
upward or downward in its run. 



62. Other steering devices 

The Gyro-Pilot is a gyroscopic device for steering ships auto- 
matically. It is used in connection with the gyroscopic compass and 
is controlled by the compass. 

The Directional Gyro is used in connection with the magnetic 
compass in blind flying. It is an azimuth instrument and gives a 
stable indication of direction. 

The Artificial Horizon, though not a steering device, provides an 
artificial horizon to enable the aviator to keep Ms airplane in level 
flight when the natural horizon cannot be seen. It is controlled by 
a gyro having its spin axis vertical, the axis remaining vertical 
however the plane may be tilted. 

The Automatic Pilot is a combination of the Directional Gyro and 
the Artificial Horizon. It employs two gyros, one installed with 
spin axis horizontal and the other with spin axis vertical. The 
former maintains the direction of the airplane in azimuth by 
controlling vertical rudders and the latter keeps the plain hori- 
zontal by controlling horizontal rudders. The automatic pilot can 
steer an airplane for a considerable time without attention from the 
human pilot. 

Another type of automatic pilot (under development or recently 
completed) for airplanes employs three rate gyroscopes, one for 
each of the principal axes of the plane. The outer gimbal of eachis 
rigidly fixed to the plane, so that angular motion (forced precession) 
about any principal axis will be clearly indicated. 
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Gyroscopic gun sights enable antiaircraft to hit 

targets with ease and accuracy. They calculate automatically 
instantaneously the proper aiming of the gun in order to hit the 
target. 

The controlling element in the gun sight is a rate gyroscope 
which computes the 'lead' of the gun required for hitting the 
moving target.* In aiming the gun the spin axis of the gyro is aligned 
parallel to the line of sight from gun to target. Then as the line of 
sight is rotated at an angular velocity dtjrjdt to keep up with the 
moving taget, the rate gyro computes the lead angle <j) in accordance 
with equation (58-1). This appKcation of the gyroscope to gun 
aiming is due to Dr C. S. Draper of the Massachusetts Institute of 
Technology. 

In the first application of a gun sight of this type (the Sperry 
Mark 14) in World War II, 32 attacking airplanes were shot down 
in a single engagement by the antiaircraft guns on one battleship. 

An application of the gyroscope to the automatic firing of naval 
guns on board ship is briefly discussed at the end of Art. 72. 



C. SOME RECENT TYPES OF GYROSCOPES 

In recent years gyroscopes have undergone many improvements and 
refinements, and their applications have been widely extended. 
Refined methods of suspending the rotors and the wide use of 
electric, magnetic and electronic accessories have increased the 
accuracy and applicability of modern gyroscopes. Much research 
has been done on inertial systems of navigation, systems based 
entirely on the inertial property of a spinning gyro the property 
that keeps the spin axis pointing in the same direction with respect 
to fixed space. The inertial system is self-contained and does not 
depend on any external object or source of energy. Some of the 
more recent types of gyroscopes are described briefly below. 
Doubtless many other advances in the gyroscopic field have not 
yet been made public. 

* In order to hit a moving target with a projectile, the gun must be aimed 
at a point ahead of the target, because the target will move a certain distance 
ahead while the projectile is in flight. The distance from the target to the 
point aimed at is called the lead. The amount of required lead depends on 
the distance from gun to target and on the speed and direction of flight of 
the target. 
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63. The Sperry Mark 22 gyrocompass 

This is the smallest gyroscopic compass yet developed. It is 
simple and rugged in construction, occupies less than half a cubic 
foot of space s and weighs 9 pounds. It has been designed for use on 
small vessels: of all types speed boats, landing craft, amphibious 
vehicles, etc. and will indicate true north when a speeding boat 
is making fast and sharp turns. In performance, the Mark 22 
compares favorably with the large standard gyroscopic compasses 
used on large ships. 

64. The Arma miniatiire gyrocompass' 

This is a small gyrocompass designed especially for military 
vehicles, such as tanks and P.T. boats. It indicates true north with 
an accuracy of about half a degree under all vehicular speeds and 
conditions, and therefore provides accurate bearings and courses 
for land and water vehicles in uncharted areas. 

The rotor spins at 12,000 r.p.m. inside a helium-filled and 
hermetically sealed ball made of stainless steel. The ball floats in 
neutral equilibrium in a fluid inside a fluid tank, and is held in 
position by very fine wires. 

The Arma Miniature Gyrocompass is designed to operate at 
temperatures from 65 to -f 130 F. and to give reliable indications 
in all latitudes up to 75. Important features of this compass are its 
ruggedness, small size, light weight (67 pounds) and reasonably 
high accuracy. 

The Arma Corporation also has under development a very small 
and precise gyroscope of a similar type for use in an inertial system 
of navigation for airplanes and guided missiles. 

65. The Draper hermetic integrating gyro (RIG) 

This is a very small one-degree-of-freedom gyroscope of high 
sensitivity, high accuracy and high resistance to shock and vibra- 
tion. It is used mainly for platform stabilization and rate measure- 
ment applications. 

The original research leading to the design of the HIG unit was done 
under the direction of Dr. C. S . Draper at the Massachusetts Institute 
of Technology. The development was continued by the Minneapolis- 
Honeywell Regulator Company, by whom it is manufactured. 

The chief components of the HIG are the gyro wheel, a torque 
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generator, a signal generator, and a viscous damper. The gyro 
wheel is mounted in a sealed container which floats in a viscous 
of the same average density as the gimbal. The gimbal rotates on 
jeweled bearings. A magnetic pickup is mounted at one end of the 
container and a magnetic torque generator is mounted at the other 
end. The fluid is kept at constant viscosity by temperature control. 
TMs fluid opposes the turning of the gimbal about its axis, serves 
to float the gimbal and thus reduce the pressure on the gimbal 
bearings, and provides a cushion against shock. The instrument is so 
designed that the damping torque is equal to the processional torque. 

The smallest of the HIG units weighs only IJ pounds and is 
compactly housed in a small cylinder slightly more than 2 inches in 
diameter and 4J inches long. The rotor spins at 8000 r.p.m. 

Gimbal displacement is proportional to the time integral of the 
input turning rate, and the secondary voltage output of the signal 
generator is proportional to the angular displacement of the gimbal 
from a null position. The signal generator serves as the signal 
source to indicate the angular position of the gimbal with respect 
to the case. The sensitivity of the HIG is such that a null signal of 
1 second of arc can be measured. 

66. The Sparry gyrotron vibratory gyroscope 

TMs is an entirely new type of gyroscope and does not even come 
under the definition of gyroscope given early in this book. Its 
fundamental element is not a spinning rotor but a kind of tuning 
fork. It is still in the development stage, but in principle it will 
probably remain as at present. Its chief application will probably 
be for measuring the magnitude and direction of the rate of turning, 
or angular velocity, of the moving body on which it is mounted. 
For example, when it is mounted vertically on an airplane, it will 
indicate whether the airplane is turning horizontally about a 
vertical axis how fast and in what direction. 

Fig. 51 represents schematically the principal parts of the 
Gyrotron. The tuning fork is attached rigidly to the torsional rod or 
member, as shown, and the latter is fastened rigidly to the platform. 
(The platform represents the airplane or other vehicle on which the 
Gyrotron is mounted.) The torsion rod is designed so as to be in 
mechanical resonance with the tuning fork; that is, it is designed 
so that its natural torsional vibration frequency about its central 
axis is the same as the natural vibration frequency of the fork tines. 
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The fork vibrations will then produce torsional of the 

same frequency. 

The device is provided with three sets of magnetic coils, as indi- 
cated in the figure. The tine pickoff coils pick up the natural 
vibrations of the tines and send them through an amplifier. The 
amplified signals go to the magnetic drive coils, wMeh keep the tines 
vibrating at their natural frequency. The third set of coils is the 
torsion pickoff coils from which the output signal voltage is obtained, 
The tines must vibrate in the same plane, with the same frequency 
and amplitude, and in opposite directions; they must move toward 
and away from each other simultaneously. The tuning fork and its 
electrical pickoffs are all hermetically sealed in a suitable container. 

The operation of the Gyrotron is as follows: When the platform is 
not turning, the fork tines vibrate in one plane. But when the 
platform turns, it forces the vibrating fork to rotate at the same 
rate about its central axis. Then when the tines are approaching 
each other, the moment of inertia of the fork about the central axis 
is decreasing. This decrease in the moment of inertia causes the 
angular velocity of the fork to speed up in obedience to the law of 
conservation of angular momentum (Art. 15, equation (15-3)); and, 
conversely, when the tines are moving away from each other, the 
moment of inertia is increasing and thereby slowing down the 
angular velocity of the fork about its axis. * This alternate speeding 
up and slowing down of the angular velocity of the fork about its 
axis produces relative motion between the plane of the tines and 
the supporting platform, and this means that the particles com- 
posing the tines are periodically accelerated with respect to the 



* From the equation H = Io) (Art. 14, equation (14*6)), we have at = 
and since H is constant in this case, it is obvious that a) varies inversely as I. 

The moment of inertia of a particle m of one of the tines with respect to 
the central axis is evidently given by the relation [see Figs. 51 and 52 and 
equation (66.1)] 



f AJ2 I 2 

1 + - sinfc^ 

1 1 + 2 sin w x t \ , practically, 



since (A.R/JR) 2 is negligible hi comparison with 1. The moment of inertia of the 
tines about the central axis is given by a similar expression. 
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platform and also with respect to fixed space. The periodic accelera- 
tion of the mass particles of the tines produces an alternating 
twisting moment about the axis of the torsion rod and thus causes 
the torsion vane to oscillate between the torsion pickoff coils. The 
vibrations of the vane cause electric signals to be sent out from the 
torsion pickoff coils. 

If the platform reverses its 'direction of turn (from clockwise to 
counterclockwise, or vice versa), all forces will be reversed in 
direction. This will cause the phase of the torsion vane signals to 
change by 1 80. The change in phase is detected by a phase detector, 
which compares the phase of the vane signals with that of the tine 
vibrations. 

The tuning fork is designed so as to have a high natural rate of 
vibration (of the order of 2000 oscillations per second), thereby 
helping to isolate it from extraneous vibrations. 

Since loss of energy in the torsion member is very small, the 
torsional vibrations must be damped in order to obtain a favorable 
time constant. 

To find the angular velocity of the torsion vane with respect to 
the platform, let us consider the motion of a particle m of one of the 
fork tines. Because of the extreme rapidity of tine vibration in 
comparison with the angular velocity of the platform, the path of 
m will be practically a very compressed sine curve, as shown in a 
very expanded and magnified form in Fig. 52. Let R denote the 
distance from the central axis of the fork to the center of each tine 
when at rest, and let A.fi denote the amplitude of the fork vibra- 
tions. Then the distance of m from the central axis at any time t is 
evidently given by the relation 

r=J8 + (AjR)siney 1 ^ (66-1) 

from which -=- = (AJS) c^ cos o)^, (66*2) 

Cut 

where 277/% denotes the vibration period of the fork. 

When the platform is turning, the particle m has a compound 
acceleration whose radial and transverse components are 

fflr 

a r=jp-~ m \ (66-3) 

do) rt dr 
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where a> denotes the angular velocity of the platform. The 
reaction force on m perpendicular to r is therefore ma e , the 
moment of this force about the central axis is 



T = rma e 




(66-5) 



Fig. 52 

Although the angular velocity of the platform is in general 
neither constant nor strictly periodic, it can be represented very 
approximately by a Fourier series, and for our purpose it is sufficient 
to use only the fundamental harmonic of the series. Hence we take 

(y = flioos6> *, (66-6) 

where O denotes the amplitude and oj^rr the fundamental 
frequency of the platform rotation. From this we get 

da) ~ . , 

-=T-= Q6 

dt 
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Substituting into (66-5) the values of r, drjdt, &> 5 d<y/dft, and simpli- 
fying slightly, we get 



T = -! 



+ 2(AJR) cos o^f cos (orf + 2(AJ2) i ^=- 1 sin % cos o^tf cos w 1 . 

(66-7) 



Now since o% < % and AjR <^ -R, the first and third terms in the right 
member of (66*7) are negligible in comparison with the second term 

and may therefore be neglected. The torque T thus becomes 

. (66-8) 



The product of cosines can be transformed to a sum by means of the 

formula 



CQg a CQg & = | cog ^ a + ^ + ^ cos ^^ 

Making this transformation, we get 

T = - mR(&R) Q,^ [cos (^ + a) ) t + cos (a)^ - 6> ) fj. (66-9) 

A similar equation holds for all particles of the tines. Hence the 
reaction moment for the whole fork is 



T = - MR(&R) Qw 1 [cos (a)! + w Q )t + cos (% - 6> ) t], (66-10) 

where Jf denotes the mass of the fork. 

We are now ready to write the differential equation of motion of 
the torsion vane. Let 

/ = moment of inertia of fork about central axis, 
<j> = deflection angle of vane at any time t, 
c = damping coefficient, 

A = torsion coefficient (torque required to twist the member 
through 1 radian). 

Then since the reaction torque, the damping torque and the 
moment of the torsional resistance all act to oppose the increase of 
the deflection angle, we have by equation (20-3) 



, 0: A, JfJ2(AJB)Q(u l 

or +-f~ = - -l - 1 



(66-11) 
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Since the vibration frequency of the vane would be the as 

of the tines if there were no damping, it is plain* that A/J = 6>f . Then 

Jf j?( AJ5) fib)* ~ 
ii we put c// = 2A/, = =Q, ft> 1 + (t>0 = -4 and ()%> (I)Q = 

for convenience, we can write (66-11) in the simpler form 
Here the characteristic equation is 



from which r= ki*J(to\ P) if 

Hence the complementary solution is 

<j> = e-** [Ci sin VK - P) + C% cos VK - fe 2 ) *]. (66- 1 3) 
The natural period of oscillation is therefore 

27T 

T "7(<=^)' 

Since the complementary solution and its derivative will soon fade 
out, we do not consider it further. 

To find the angular velocity of the vane due to the motion of the 
platform, we find a particular solution of (66-12) by putting 

= <7 3 sin At + C 4 cos At + C 5 sin St + G 6 cos -BL (66-14) 
Then: 

^ = <7 8 J. cos At-CtA sin J.^ + f7 5J B cos Bt - (7 6 sin Bt, 

$ = - (7 3 ^ 2 sin -4* - (7 4 J. 2 cos At - (7 5 J3 2 sin J5 - C B B* cos 5^. 

Substituting these into (66-12) and equating coefficients of like 
functions of t on the two sides of the equation, we get 



* If we put c = in (66- 11) and denote the right-hand member by/() the 
equation becomes -, 



The natural vibration period in this case is 27r/ A /(A/l). But the natural 
vibration period of the fork tines is STT/^. Since the two periods must be 
the same when c = 0, we must have V(A/7) = <y lf or A/1 = wf . 
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Solving these equations for (7 35 <7 4 , (7 5 (7 6? and substituting their 
values into (66-14), we find a particular solution to be 



which can be written 



sin (At + a} 

' (66 - 15) 



_ 

where c^tan-i-i-, / ?= tan -i-|^-. (66-16) 



These equations can be still further simplified by replacing A and B 
by their values and then neglecting all terms in which the negligible 
fraction <y /w x occurs. Thus 



- 



Likewise, 



I-) 



. 
24(6)! - w ) 2^(6)! - 6> ) k \ 2(a! - a> ) 

Hence ff 
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On substituting the above simplified expressions into (66-15) 
(66-16), we get 

^ = 2w^7(#^^ 

(66-17) 

When the sum of the sine terms in (66-17) is transformed into a 
product by means of the formula 



rt . - 

sin x -f sin y = 2 sin ~ cos -^ , 

2 2 



the equation becomes 



The angular velocity of the torsion vane with respect to the 
platform is therefore 

I^\ cosw^cos(w ^-/?)- O sinw^sin(ft>^-^) 
"i~"W L &>i J 



since ^) / w i * s negligible. 

Finally, when Q and k are replaced by their values, the angular 
velocity becomes 

)! * , * m /M ,x 

}^ (66 " 19) 

This equation shows that the amplitude of the angular velocity of 
the torsion vane with respect to the platform varies directly with 
the amplitude of the angular velocity of the platform with respect 
to fixed space. 

It is to be noted in (66-18) and (66-19) that <j> and <j> are modulated 
quantities* of essentially the same form, but that the amplitude of 
<f> is &>! times that of ^. A torsion pickoff sensitive to velocity may 

* A modulated quantity is a combination of two or more oscillating 
quantities which produces new frequency components not present in the 
original oscillating quantities. The modulating factor in (66-18) and (66*9) 
is sin (D-L t. 
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therefore be utilized. This type of pickoff is not subject to the 
resolution errors of a displacement-sensitive pickoff. 

Some of the advantages of the Gyrotron Vibratory Gyroscope 
are absence of bearings and other wearing parts, small size, rugged- 
ness, and long life expectancy. There is no maintenance problem 
and the device is sensitive only to rotation about the central axis. 

For a fuller description of the Gyrotron and further technical 
details of its construction and operation, the reader is referred to 
two papers in the Aeronautical Engineering Review for November 
1953. 



CHAPTER VII 
The Gyroscope as a Stabilizer 

A. THE GYROSCOPIC SPHERICAL PENDULUM 

67. Definitions 

A spherical pendulum is a pendulous body suspended from a 
fixed point and free to swing in any manner about that point, 
subject only to the condition that its center of gravity remain at a 
constant distance from the point of suspension. The locus of the 
e.g. is thus a spherical surface having the fixed point as center. The 
ordinary pendulum which swings in a vertical plane and the 
conical pendulum which swings in a horizontal circle are special 
cases of a spherical pendulum. 

A gyroscopic spherical pendulum is a spherical pendulum wherein 
the pendulous body is a gyroscope and its mountings, so mounted 
that its axis of spin coincides with the straight line joining the point 
of suspension and the e.g. of the gyro. The rotation of the earth has 
a slight effect on the motion of the gyroscopic pendulum, but that 
effect will be neglected in the following treatment. 

68. Differential equations of motion 

To derive the differential equations of motion of a gyroscopic 
spherical pendulum when slightly disturbed from its equilibrium 
position, we take a set of fixed rectangular axes X YZ and a set 
of moving rectangular axes X'Y'Z', the origin being at the 
point of suspension of the pendulum. The position of the moving 
trihedron at any instant can be specified by the two angles <f> and 6, 
as indicated in Fig. 53. If the moving trihedron originally coincided 
with the set of fixed axes, it can be brought to the position shown in 
the figure by first rotating it about OX through the angle <j) and then 
rotating it about Y r through the angle 6. The angles $ and 6 are 
assumed to be small,* and the gyro is assumed to be spinning with 

* The treatment of the general case, where <j> and 6 may be of any magni- 
tude, requires the use of elliptic functions and will therefore not be given here. 
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Fig. 53 

a constant high velocity &>. The angular velocities ^ and 6 are also 
assumed to be small and are indicated by vectors in the figure. 

The components of the angular velocity of the moving trihedron 
along its own axes are seen to be as follows: 

&V = - ^ cos 6, Wy, = 6, o)^ = - (j> sin 6. 
Hence the angular momenta along these axes are 

H x ' = A(f> cos 6 9 Hyt = A6, H^ = C(o) <j> sin 0) . 

Since sin# and $5 are both small, their product is negligible in 
comparison with o). Hence we may write 
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Then we have 



= A $ cos 6 , practically , 



Substituting into equations (16-4) the quantities found above, 
then neglecting all terms containing either $4 a or %5, and putting 
cos0=l,weget - Carf=J 



To find My, and jM^/, let W denote the total weight of gyroscope 
and mountings, and let I denote the distance from to the e.g. of 
the pendulum. Then since <p and & are assumed to be small angles, 
the moment arms of W about the x- and y-axes are Ifi and W, 
respectively. (The reader is reminded that in Fig* 53 the pendulum 
is in the octant diametrically opposite the first octant.) 

The components of W parallel and perpendicular to OX' are 
W sin 6 and W cos 6, respectively; and the components parallel and 
perpendicular to OT' are TFsin^ and IF cos 0, respectively. Only 
the perpendicular components W cos 6 and W cos (j> have moments 
about the two axes mentioned. 

Since Y r and Y are in the same vertical plane, the moment arm 
of W cos (f> about Y' is W. The moment arm of W cos d about OX' 
is l<p + W tan 6 sin 56, or l<f> practically. 

Now bearing in mind the convention as to the signs of moments, 
we have = m ^ practicallj? 



, = - (W) W cos <j> = - WW, practicaUy . 
Substituting these into equations (68-1), we get 
(a) 



(68-2) 
(6) Ad + Cuj + WI0 = 0.) 

These equations determine the subsequent motion of the pen- 
dulum when it is slightly disturbed from its equilibrium position. 
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69. Solution of the equations 

To solve equations (68-2), we eliminate one of the variables and 

thereby reduce the system to a single fourth-order equation in the 
other. To eliminate &, for example, we first differentiate (6) with 
respect to t and then eliminate 6 between the resulting equation and 

(a), thus obtaining 

WJ 

(c) A$+Cu$+j(A$ + Wty) = Q. 

Lf(s) 

We next differentiate (a) twice with respect to t and eliminate 
S between the resulting equation and (c), thereby obtaining 

Aty* + (CW + 2AWty$+W*P4> = Q. (69-1) 

The same procedure will give a similar equation in 6 and its 
derivatives. 
The auxiliary equation for (69-1) is 

A*r 4 + ((72^2 + 2A Wl) r* + WH* = 0. (69- 1 a) 

Hence 

2 _ _ (Qfyja + 2A Wl) 



2A* 

CQJ V(C 2 6> 2 + M. Wl) 

' 



Since the radical in the above expressions is positive, the two 
values of r 2 are both real; and they are both negative, because the 
radical is less than (7 2 &> 2 + 2 A Wl, Let these values of r 2 be 1c\ and 
-if. Then 

Wl) + Cw 



2 _ M _ - (CW + 2A Wl) - Co) V((7 2 ( t ) 2 + 4A Wl) 
r,--*,- 2l* 



Hence ^ = fci, r z = ik 2 , 

and 



C 2 2 + 2AWI + Ca> V(<7 2 w 2 + A W I) 
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Hence CW + 2AWI 

l* + K> = - _ - (69-4) 





(69-5) 



- 



- 
A 2 ' 



or Mi = -3-. (69-6) 

From (69-4), (69-5) and (69-6) we iind by elementary algebra 

(59-7) 



7 7 /M v 

Hence ^2 + ^ = -^ - -^ -- ', (69-8 a) 

*a-*i = ^- (69-86) 

The values of & 2 and i x can be reduced to simpler forms as follows : 



Expanding the radical into a binomial series and retaining only the 
first two terms of the expansion, we have 

2AWl ~ Gt AW 



A gyroscopic spherical pendulum may be of any length, but when 
used as a stabilizer it must have a long period of precession and this 
necessitates that the quantity AWl[C*b)* be negligible in com- 
parison with 1. This in turn necessitates that I be as small as 
practicable. On neglecting AWl\C*uP in the above value for & 2 , 
we get 

4 2 = . (69-9) 
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Proceeding in a similar manner with (69-8), we have 



Here we are not Justified in neglecting .4 WljG 2 ^ 2 in comparison 

with, zero. Hence we retain it and get 



A glance at these values of k^ and Jc 2 shows that when o) is large 
(as it always is in stabilizers) 2 is large and ^ is very small. Hence Tc^ 
is negligible in comparison with k 2 . 

Since the values of r are ik- L and iJc^ the solution of equation 

(69-1) is ^ =ClS i n jfc 1 + c 2 cosfc 1 ^ + C3sini 2 + C4Cosi 2 ^ (69*11) 

To find 9 without bringing in additional constants of integration, 
we first differentiate (68*2 (a)) with respect to t and then eliminate & 
between the resulting equation and (68-2 (6)), thereby obtaining 

0= --^[A*$+(C*(i>*+AWl)fa (60- 11 a) 

Then we find the values of $ and $ from (69-11) and substitute them 
into the above equation, thus obtaining 



+ A Wl)] cos 



- (0*a)*+AWl)] cosk 2 t}. (69-12) 

To simplify the coefficients of the trigonometric terms we replace 
&f and i a by their values from (69-2) and (69-8), respectively, and 
thus get 



and 1sA*t% - (G 2 a> 2 4- A Wl)] = - Co> Wl. 

From (69-3) and (69-7) we find in a similar manner that 

kJA*t% - (CW+AWlf] = Ca)WL 

On substituting into (69-12) these values of the coefficients, we get 
6 = c 2 sin Jc l i-c li cosfc^-^ sin k 2 t + c s cos k%t. (69-13) 
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Equations (69*11) and (69-13) each represent two 
and superposed simple harmonic motions of periods 

m %7T Cd> 

2^ = _ =: 2?r == , the precession period, 

and T 2 = = 2w ^ , the nutation period. 

TI is thus a long period and jP 2 a very short period. The motion for 
either (j> or 6 is therefore a pseudo-regular precession as discussed in 
Art. 38. 

If in (69-11) we decrease t in the i^-terms by a quarter period 

I - -r- j , and increase t in the i 2 -terms by a quarter period | - j- \ , 
the right-hand member becomes 



c x cos Tc^t + c 2 sin fc x ^ -f c s cos Jfc 2 ^ 

which is the value of 6 as given in (69-13), This result shows that in 
the long-period motion the phase of 6 is a quadrant behind ^ and 
in the short-period motion it is a quadrant ahead of ^. 

We now determine the constants of integration from assumed 
initial conditions. From (69-11) and (69-13) we get 



= c-Jc-L cos k^t c 2 ^ sin k l t + c 3 & 2 cos k 2 t c 4 & 2 sin k%t, 

(69*14) 



(69-15) 

Let us assume that while the pendulum was hanging in equili- 
brium vertically downward, with gyro spinning at constant speed 
a), it was struck a horizontal blow which imparted to it an angular 
velocity $ in the direction of the negative #-axis. Then the initial 
conditions are 

= 0, # = 0, $5 = 0, <9 = $ when = 0. 

On substituting these values into (69-11), (69-13), (69-14) and 
(69-15), we find that 

Cl==0> Ca= ~ 
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Hence equations (69-11) and (69*13) become 

A 



(cos kit cos k 2 t), 



A 



"Now neglecting ^ in comparison with 2 and replacing & 2 by its 
value cd)jA from (69*9), we have 

AS 

, (69-16) 



). (69-17) 

l_/U/ 

These equations show that the disturbed motion is stable and that 
the amplitude is very small when o) is large. 

It is of some interest to determine the constants of integration in 
(69-11) and (69*13) for a different set of initial conditions. Let us 
assume that the pendulum is swinging in the plane = and is 
struck a sharp blow at right angles to its plane of vibration at the 
instant it reaches the end of its swing. If the blow imparts to the 
pendulum an angular velocity in the direction of the negative 
t^-axis, the initial conditions are 

Substituting these quantities into equations (69-11), (69-13), 
(69-14), (69-15), we find 

_0o~Mo _ A _4 + Mo , _A 



Then equations (69-11) and (69-13) become 



1 ' 2 "'I "" 2 

Since i x and ^ are both small, their product is negligible in 
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comparison with $5 ; and since k t is negligible in l* 2s 

the above equations become 

2 - (69*18) 

r*AQ Z* / .1 rAQ Z* / 

V/WQ IV -IV T^ y \J\J'St /Vfttm 

These equations can be simplified still further in view of the fact 
that 5 /&2 ^ s negligible in comparison with . Hence we neglect 
this fraction in equations (69*18) and get 



= 00008^. J 



(69-19) 



These equations show that in this case the motion of the pendulum 
is a slow precession in a circle of radius W Q . 

70. The spherical pendulum without gyroscope 

If we put a) = in equations (68-2), we get the differential equa- 
tions of motion of an ordinary spherical pendulum. We thus have 



(70-1) 
The auxiliary equation for either of these equations is 



r fa 



;Wl 
r 



Hence the solutions are 
(a) (fr^c-L sin let + c 2 cos let, 



Then 

(c) ^^Cj 

(dJ) d c^lc cos & c 4 sin &. 
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For the Initial conditions 

^ = ? = 0, ^ = 0, $ = $o when = 5 

we find ^==0, c 2 = G, c 3 = ^ i, c 4 = 0. 

Hence the equations of motion are 






The pendulum thus vibrates with simple harmonic motion in the 
plane $ = 0. It will be noted that the amplitude is here much larger 
than in (69*17), because of the presence ofa) in the denominator of 
the amplitude of the latter. 
For the initial conditions 



we find from (a), (6), (c), (d), 

A. \ 

' C2=<>1 C3=0) c * =9 ' 



Hence equations (a) and (6) become 

(70-3) 




On multiplying these equations through by I and recalling that 
l(j>=y and 16 =x, we have 



X W Q GQS \-J-\t, 

which are the parametric equations of an ellipse with semiaxes 

I A \ 

I == I and Z# A . 
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The period of revolution is 

T = 2 



In the gyroscopic pendulum the long period is 2??^= , which is very 

yj I 

much longer because of the factor a). 

71. Numerical comparison of the amplitudes of the gyroscopic and 
spherical pendulums 

As a numerical case we consider a gyro or rotor 8 inches In 
diameter, having a rim 3 inches wide and 2 inches thick, and a web 
| inch thick. Assuming that the distance from the point of support 
to the e.g. of the pendulum is 3 inches, that steel weighs 490 pounds 
per cubic foot, and that the rotor is spinning at 10,000 r.p.m., we 
have the following data: 

Z = Jfoot, 
W =33.8 pounds, 

o)= 10,000 r.p.m. = 838 radians per second, 
(7 = 0.070, in pound-feet-seconds units, 

33 8 /1\ 2 

A = 0.040 + ' \~\ = 0. 106, in pound-feet-seconds units. 
32.2 \4/ 






To = = 0.011 second, 
554 

21 = J^-T = 43. 6 seconds. 
1 0.144 



It will be recalled that the equations of motion for both the 
gyroscopic pendulum and the spherical pendulum were derived on 
the assumption that <j> and 6 were to be small angles. Hence the 
assumed value of $ must be such that <p and 6 for both types of 
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pendulum will not exceed 3 3 say. A preHminary calculation shows 
that $ must not exceed | radian per second* Hence we take 
4j = | radian per second, and substitute this value in equations 
(09-16), (69-17) and (70-2), From (6946) and (69-17): 



(cos 0. 144 - cos 554) 
58,7 

= 0.00090 (cos 0.144$ - cos 5542), 
6 = 0.00090 (sin 0.144 + sin 554). 

Since the algebraic sum of the trigonometric terms in parentheses 

can never exceed 2, we have 



? 0.001 8 radian = 6', approximately. 
In the case of the ordinary spherical pendulum we have, from 
/0106 

= v ^rr ' 056 radian==3 - 2 = 192 '- 



The amplitude in this case is thus 32 times as great as in the 
gyroscopic pendulum under the same initial conditions. The great 
stability of the gyroscopic pendulum makes it of great value as a 
stabilizer. 
The period of this spherical pendulum is 



or about ten times as great as the short (nutation) period of the 
gyroscopic pendulum. 

72. Gyroscopic pendulum with point of suspension subjected to a 
periodic force 

When a gyroscopic pendulum is mounted on a ship, the rolling 
of the ship about its longitudinal axis subjects the point of suspen- 
sion to a to-and-fro motion athwartship. This motion is periodic 
and produces a forced vibration of the pendulum. 

To investigate the stability of the motion of the gyropendulum 
in this case, we reduce the problem to that of the nonrolling case 
by applying d'Alembert's principle (Art. 11) and then proceed as 
in Art. 68. In the case of the gyropendulum the external applied 
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force is unknown, but the inertia force is Wa/g acts through 
the center of gravity of the pendulum. It thus has a moment about 
any axis not parallel to it and not intersecting it. 

Referring to Fig. 53, let us assume that OX is parallel to the 
longitudinal axis of the ship. Then the periodic motion of the 
point of suspension of the pendulum is in the direction of T 
may be represented by a sine function of the time. Hence we 
assume that the motion is represented by the equation 

y=a saute. 
Then y 



The inertia force is therefore 



and acts to the right parallel to OF. Hence for the moment of this 
force we have 



Adding these moments to the gravity moments of Art. 68, we get 

W \ 

(a) A6 - Cud + WU = abH sin U, 

-9 \ (72-1) 



(b) 

as the differential equations of motion in this case. 
On eliminating 6 as in Art. 69, we get 

- sin bt. 



(72-2) 
To find a particular solution of (72-2) we put 



(72-3) 
where P and Q are constants to be determined. Then 

$ = - Pb 2 sin bt - Qb* cos bt, 
0iv = PW sin bt + Qb* cos bt. 
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Substituting into (72-2) these values of ^>, $ and ^ lv and then 
equating the coefficients of sin bt and cos bt on the two sides of the 
resulting equation, we get 

(c) 



(d) QAW - (OW + 2 A Wl) Qb* + Q WH 2 = 0. 

From (c) we get 

(e) Jr - 



- (CV + 24 Wl) 6 2 + fF 2 Z 2 ] ' 

and from (d) we have 

= 0. 



Here either Q or the expression in brackets must be zero. We first 
examine the bracketed expression. 
On putting r 2 = k 2 in (69- la), we get 

(CW + 2A Wl) F + WW = 0, 



which is the bracketed expression in (/) if b is replaced by k. Hence 
the bracketed expression cannot be zero unless b = Jc, which means 
that 6 X = A! and 6 2 = & 2 . But since the rolling period of a ship is not as 
slow as the period given by ^ nor as fast as the period given by & 2 , 
we must conclude that the expression in brackets is not zero. 

Hence we must have n 

y = u. 

Substituting into (72-3) the values of P and Q just found, we have 



as a particular solution of (72-1) for <j>. 

To find the corresponding particular solution for 6, we differen- 
tiate (72-1 (a)) with respect to t and then eliminate 8 between the 
resulting equation and (72-1 (6)). We thus find 



(72-5) 



We next find <f> and $ from (72-4) and substitute their values into 
(72-5), thus obtaining 

- - CcoWalb* _ 

** - (C V 5 4- 2 A Wl) 6 2 + F 2 Z 2 ] COS bL ( 6) 
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Now since 

AW - (CW + 2A Wl) b 2 + WH* = - [C 2 w 2 6 2 - ( Wl - Atff}, 
we can write (72-4) and (72-6) in the simpler forms 

, , m 

( } 



(72-8) 

These equations show that the rolling of the ship produces forced 
vibrations of the pendulum not only athwart the ship but also in 
the direction of the longitudinal axis. The presence of (*) in the 
numerator of (72-8) also shows that the amplitude of & is much 
greater than that of <j>. 

On replacing bt by bt |TT, we find that (72-7) takes the form of 
(72-8), which shows that the ^-vibration is a quarter period behind 
the ^-vibration. 

Because of the presence of w 2 in the denominators of (72*7) and 
(72-8) the amplitudes of the forced vibrations are very small. If we 
put o) = 0, the equations become 

WalW . ,\ 

:**,! (72 . 9) 



5 = 0, 

which show that in this case the pendulum vibrates only in the plane 
perpendicular to the axis of roll. 

When the gyro spins at very high speed, as it always does in 
stabilizers, the quantity (Wl-AV^ is negligible in comparison 
with 2 <w 2 6 2 . Hence on neglecting (Wl JL6 2 ) 2 in the denominators 
of (72-7) and (72-8), we get 

~7/TTr7 /fZ.2\ 

(72-10) 



__ 

(72-11) 



guv 

These equations show at a glance that the amplitudes of the 
forced vibrations of a gyroscopic pendulum are very small, and by 
comparison with (72-9) they are seen to be much smaller than the 
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forced vibrations of the same pendulum with the gyro not spinning. 
Hence the spin axis of the gyroscopic pendulum remains practically 
vertical in spite of the rolling of the ship. 

Note. The condition that a gyropendulum be unaffected by the 
motion of the ship or airplane on which it is mounted is that 
Go)/Wl = ^j(R/g), where B denotes the radius of the earth. For the 
derivation of this condition, see the Appendix. 

An important application of the gyroscopic pendulum is in the 
automatic firing of naval guns on a rolling and pitching ship. The 
gyroscopic pendulum maintains a horizon of reference. The aiming 
device of the gun is set at the required elevation above the horizon 
for hitting the distant target. When the roll of the ship brings the 
axis of the gun barrel to the preset angle of elevation, electric 
contacts close a circuit that fires the gun. In case the roll of the ship 
is not sufficient to raise the gun to the required elevation above the 
horizon, the preset angle of elevation is made great enough to allow 
for the insufficient roll. 



B. OTHER TYPES OF GYROSCOPIC PENDULUMS 

73. Gyroscopic pendulum suspended at center of gravity of gyro and 
with pendulous weight below 

Gyroscopic pendulums may take several forms. A common type 
is that illustrated in Fig. 54. It will be seen that the gyro is suspended 
in Cardan rings, the outer ring resting on pivots A and A' attached 
to supports not shown in the figure. The pivot axes A A' and BB r 
are perpendicular to each other and intersect at the center of 
gravity of the gyro. The plane of the gyro is kept approximately 
horizontal by means of a pendulous weight W. 

The mounting permits the pendulum to swing through an angle 
<j> about the axis A A' and through an angle about the axis BB' . 
Pig. 55 shows the position of the pendulum and pivot axes after the 
pendulum has been disturbed from its position of equilibrium. (In 
this figure the pivot axes of Fig. 54 have been rotated 90 to the 
right.) 

Comparison of Fig. 55 with Fig. 53 shows that the equations of 
motion in this case are exactly the same as those of the gyroscopic 
spherical pendulum found in Arts. 68 and 69. The derivation of the 
equations of motion will therefore not be repeated here. 
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In a numerical example, however, there is one Important 
difference between the gyroscopic pendulum of the present article 
and the gyroscopic spherical pendulum previously treated. That 
difference is in the transverse moment of inertia A. In the case of 
the gyroscopic spherical pendulum, A must be found by means of 
the parallel-axis theorem relating to moments of inertia, namely, 




Fig. 55 



where A@ denotes the transverse moment of inertia of a body of 
mass M about an axis through its center of gravity, and A L is the 
moment of inertia about any other axis L parallel to the gravity axis 
and distant d from it. The above formula was used in computing A 
in the example of Art. 71. In the present case, A is simply A o . 

Gyroscopic pendulums mounted as in Fig. 54 are used in 
stabilizing ships and airplanes. 

74. The inverted gyroscopic pendulum 

Gyroscopic pendulums with pendulous weight above the point of 
support have application in monorail cars. Fig. 56 shows such a 
pendulum when displaced from a vertical position. The gyro may 
be mounted as in Fig. 54 with its e.g. at the point of support, or it 
may constitute the pendulum itself as in the case of the gyroscopic 
spherical pendulum. The mathematics of the two cases is the same 
except for the transverse moment of inertia A. 
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From Fig. 54 we have 

<t) x , = - (j> cos 0, Mf = &, (1)*$ sin 0. 
Then fi^ = -4^cos0, Hf=A& 9 H s ^G(~<j>mnd) 

= Cfe>, practically. 



Hence 



=- - ^ ( jj cos 6 - $?$ sin 19) = - A $ cos 5, practically, 



Also, J4,= 

Substituting the above quantities into (16-4), we get 



A8+ Aft sin 6 cos 6 + Ca>< cos 6 = WW, 

- A6<j) cos + A^cos6 = Q. 
Neglecting the terms which contain ft and 6(j> and putting cos 6 1, 

^^ W ^-^=^1 

(6) ^+<7<w#= s TFM.J 

It will be seen that these equations differ from (68-2) only in the 
sign of I. 

To solve equations (74-1) it is not necessary to go through the 
process of eliminating 6. To get the fourth-order equation in <f> that 
results from such elimination, we simply change the sign of I in 
equation (69-1). We thus get 

- 2A Wl) $ + WH*<j> = 0. (74-2) 



An exactly similar equation holds for 6. 

The auxiliary equation for either equation is 



2^2 _ 2 A Wl) r z + WH* = 0, 
from which 

2 A Wl) 



_ 24 Wl) Co) V(C 2 co 2 - 4A Wl) 
2A 2 
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IQ 




Fig. 56 



Since both values of r 2 are real and negative when C*o) 2 > 4J. Wl, we 
write 





Inen 



_ 2A Wl + Co) V(C 2 OJ 2 - 4A Wl) 
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From these last two equations we get, by the same procedure as in 
Art * 69 ' 



2A 



*- 2 A 

By replacing the radicals by the first two terms of their binomial 
expansions, we get the simpler results: 



^ 

A ' 



The solution of (73-2) is therefore 

= c x sin ^ + c 2 cos ^ + c 3 sin k 2 t + c 4 cos A 2 L (74-3) 

An exactly similar equation holds for 0, except for the signs and 
positions of the c's (see Art. 69). 

We could determine the integration constants in (74-3) by the 
same procedure as in Art. 69, but it seems unnecessary to do so. 
The important result we have arrived at is that, since the values of 
<f> and 6 do not increase with time, the motion of the inverted 
gyroscopic pendulum is stable when the rotor spins fast enough to 
satisfy the condition C 2 o> 2 > 4 A WL 

75. Gyroscopic pendulum with gimbal axes rotating in azimuth 

When a gyroscopic pendulum is mounted on board ship or on an 
airplane and the ship or airplane turns in azimuth, the gyroscopic 
pendulum turns with the ship or airplane. Moreover, gyroscopic 
devices are sometimes mounted on platforms which can be rotated 
in azimuth in either direction with respect to the deck of the ship. 
Hence it is desirable to investigate the behavior of a gyroscopic 
pendulum when it is compelled to rotate in azimuth. 

Let us consider a gyroscopic pendulum of the type treated in 
Art. 73 and assume that the gimbal axes A A' and BB' (Fig. 55) are 
rotating in azimuth with constant angular velocity ijr in the same 
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direction as the gyro is spinning {see Fig. 57). To derive the differen- 
tial equations of motion in this case, we find the projections of ^ 
on the moving axes, add these to the angular velocities that 




Fig. 57 

when there is no rotation, and then use these total angular velocities 
in connection with equations (16-4). Thus, in Fig. 57 the rectangular 
components of ft in the 7^-plane are i/r cos along OQ and - ft sin ^ 
along Y'. Then ft cos <p has the components ft cos cos 6 along OZ f 
and -ftcos<fisin8 along OX'. Hence the projections of ft on the 
moving axes are ^ COS95sin on OX ' 9 

ftsinfi on OT f , 

on OZ'. 
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Since the angular velocities about OX\ OY f and OZ' when the 
gimbals are not rotating are $5 cos 6 9 6 and $ sin 6, respectively, 

the total angular velocities about these axes when the gimbals are 
rotating in azimuth are 

yf ~ "~ ^ COS d ty cos s * n $? 



<y s , = $ sin -f ^ cos cos 5. 
Hence JZ^ = -4 (0 cos 8 + i/r cos <j) sin 



fig, = <7(w 5 sin 8 + ^ cos $ cos $) . 
Then, remembering that t/r is a constant, 

JTT 

-TT^= -4(^ cos ^$sin0 ^ sin $5 sin <9 + ^0 cos $5 cos 
GW 

= A<j) cos 0, approximately, 

= J[ (^ i/r$ cos <fi)=zA@ y approximately, 



^ = O( ^ sin 6 <j)6 cos ^ sin ^ cos 6 i]f6 cos $5 sin 0) 
ttt 

= C^ sin #, approximately 

= 0, practically, since 6 is small and (j) changes slowly. 

On substituting the above quantities into (16*4) and neglecting all 
terms containing the squares and products of ^, 6 and ^, we get 

- A$ cos 9 + Cud - G^ sin < = M X ,= Wl<f>, 
A0 + Cdxj) cos 6 + OtL)^ cos (f) sin (9 = If^/ = - If Z0. 
Now putting sin <f> = <f>, sin 6 = #, cos ^ = 1 and cos = 1 , we have 

+ Wl) (f> = 0,^ 
+Ifi!)<9 = Oj 

It will be noted that if ^ = 0, these equations become identical with 
(68-2), as they should. 

It will be further observed that equations (75-1) are identical in 
form with (68-2), the only difference being that the constant 



STABILIZER 183 

coefficient Wl in (68*2) has been replaced by the constant 

C&i/f-i- Wl. Hence we may utilize the solutions of (68-2) to 

the solutions of (75*1). For example, we have, by equations (69-8a) 






< 75 ' 3 > 



And &! and k 2 in their final reduced forms are, from equations (69-10) 
and (69-9): 



The periods are therefore 



A 
2= n Oo)' 

These results show that rotation of the gimbal axes in the same 
direction as the spin of the gyro makes the long period shorter. 

If the gimbals are rotated at constant angular velocity ^ in the 
direction opposite the spin of the gyro, the equations of motion are 
obtained from the preceding by changing the sign of i/r. The 
differential equations of motion are then 

(75-4) 



The values of Jc v Jc%, and the periods in this case are 

., Wl-Guifr 



4/1 Co* 
fca= j7> 
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Here the long period is longer than is the case when the axes are not 
rotating. 

The reader will note that in the case of opposite rotation here 
considered, we must have Wl > C&ifr in order that k[ and T{ be 
positive. This fact puts a restriction on the magnitude of i/r if the 
gyroscopic pendulum is to serve as a stabilizer. Hence when turning 
in azimuth, the gimbal axes should be turned, if possible, in the same 
direction as the gyro spins. 

The solutions of equations (75-1) are 

= c x sin &! t + c 2 cos fc x t H- c s sin k 2 t + c 4 cos k%t, (75*5) 

, (75-6) 



the same as in Art. 69. To find the constants of integration we first 
find <j> and 6 from (75*5) and (75-6) and then assume the same 
initial conditions as in the first case in Art. 69, namely, 

= 0, = 0, = 0, # = # , when = 0. 

6 
We thus find = , , (cos ^ t - cos & 2 Z) , 

K>1 "T" "^2 

(sin 



1 ' 2 

On replacing ^1 + ^3 ^7 ^ s v ^ ue from (75-2), we get 

A6 



(75>7) 



) - (75<8) 

If the rotation in azimuth is opposite the direction of the rotor 
spin, the equations of motion are found from (75-7) and (75-8) by 
changing the sign of i/r and replacing Fs by k"s. Hence for this case 
we have . A 

(75 ' 9) 



(75 ' 10) 
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On comparing equations (75-7)-(75-10) with, the of no 

rotation (^ = 0), we see that rotation of the axes in the direc- 

tion as the spin of the gyro decreases the amplitude of a 'disturbed 
motion (makes the motion more stable), and that rotation in 
the opposite direction increases the amplitude of such motion 
(makes the motion less stable). These facts furnish another 
reason for making a turn in the same direction as the spin of the 
gyro. 



CHAPTER VIII 

The Gyroscope as a Stabiliser (Continued) 
The Ship StabiBxer 

A. SHIP ON A STRAIGHT COURSE 
76. Schlick or brake-type stabilizer 

The idea of using a gyroscope to reduce the rolling of a ship at sea 
is due to Otto Schlick of Hamburg, Germany, who first installed a 
large pendulous gyroscope on a ship and carried out experiments 
with it in 1903. The gyro was installed with its axis of spin normally 
vertical and with its supporting frame resting on horizontal 
trunnions athwart the ship (Fig. 58). A pendulous weight w was 
attached to the bottom of the frame in order to make the center of 
gravity of gyro and frame fall below the horizontal axis of support, 
thus making the gyro and frame gravitationaUy stable. 

From our previous study of the behavior of a gyroscope we can 
readily see that a rolling motion of the ship about a longitudinal 
axis* will cause a precessional to-and-fro motion of the gyro and its 
frame about the trunnion axis FF'. To limit the amplitude of such 
motion in the longitudinal plane, Schlick installed a braking 
mechanism B which applied a constant pressure on the axle. 
Although the brake type of stabilizer has been superseded by the 
active type invented by Elmer A. Sperry, the mathematical 
theory of the former type is instructive and will therefore be 
given here. 

* The rolling of a ship is a to-and-fro rotation about a longitudinal 
horizontal axis. If the water offered no resistance to rolling, the axis of roll 
would pass through the center of gravity of the ship. But when the water 
offers resistance to rolling, as in fact it always does, the axis of roll passes 
through what is called the /tranquil point', which is above the center of 
gravity. The rolling ship thus behaves somewhat as a pendulum swinging 
about a longitudinal axis through the tranquil point. See Peabody's Naval 
Architecture (1917), pp. 344-5. 
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When a skip floats in water, it displaces its own weight of water. 
The form of the displaced water changes as the ship rolls from side 
to side about a longitudinal axis. The center of gravity of the 




Fig. 58 



displaced water is called the center of buoyancy and changes its 
position as the ship changes its inclination from the vertical. 

Figure 59 represents a transverse section of a ship inclined from 
the vertical by an angle $ about a longitudinal axis. Here denotes 
the center of gravity of the ship and B denotes the center of 
buoyancy. When a ship is floating freely in equilibrium and inclined 
at a small angle from the vertical, a vertical line through B inter- 
sects the stable vertical line at a point M 9 called the metacenter. 
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The distance OM( = h) is called the metacentric It Is 

that the weight of the ship and the buoyant force of the 
water form a couple which tends to restore the ship to an 
position when M is above G. (If M were below G 9 the would 
capsize.) 

In the mathematical treatment of the ship stabilizer we must 
take account of friction in two places: (I) the skin and obstructional 
friction between water and ship's hull, and (2) the brake friction 
on the precession axle of the gyroscope. We must also include the 
pulsating action of the waves on the ship. 

The resisting moment of the skin and obstructional friction on the 
hull is proportional to the angular velocity of roll of the ship, and 
the moment of the brake friction on the precession axle is pro- 
portional to the angular velocity of precession. Hence with 
reference to Figs. 58 and 59, let: 

K<j> = moment of skin and obstructional friction on hull, 

W = moment of brake friction on precession axle, 
M sin bt = moment of wave action about axis of roll, 

1 = moment of inertia of ship about axis of roll, 

A = moment of inertia of gyro and frame about trunnion 
axis, 

W = weight of ship, 

w= pendulous weight attached to gyro frame, 

h = metracentric height of ship, 

I = moment arm of w. 
Then 

W h sin ^ = W h<j> = righting moment about axis of roll, 

wl sin 6 = wW = moment of w about trunnion axis. 

Now since the axis of roll is in the direction of the #-axis and the 
trunnion axis is in the direction of the y-axis (see Figs. 55 and 58), 
we have 



Hence 

H y .=A6, 
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since <j> and 6 are assumed to be small. Then 

JTJ 

rp: = - /($ cos - <$ sin 0) = - J^" cos ^ y 




Fig. 59 

Now bearing in mind the directions of the several moments and 
substituting the above quantities into equations (16-4), we have 



AS + Caxfi cos = - wW - k6, 
or, since cos#=l approximately, 
(a) 
(6) 



(76-1) 



These are the differential equations of motion of the ship about its 
axis of roll and of the gyroscope about the trunnion axis. 
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To eliminate $ we first differentiate (76-1 (a)) and get 

(c) 1$ + K$ - Cod + Wl4 = Mb cos bt. 
Then from (76-1 (6)) we get: 

4=~ 



Substituting these into (c), we have 



+ (Kwl + Whk) d + WJiwW = - CuMb cos 61. (76-2) 
Here the auxiliary equation for the left member equated to zero is 
Air* + (Ik + A K) r 3 + (Iwl +Kk + CW + WAh) r 2 

+ (Kwl + Whk) r + WJiwl = 0. (76-J) 

To eliminate ^ and get a fourth-order equation in 5, we proceed in 
a similar manner by first differentiating (76-1 (6)) and then substi- 
tuting into the resulting equation the values of 6, 6 and 8 obtained 
from (76-1 (a)). The result is 



+ ( Whk + Kwl) <f> + Whwl<j> == M(wl - AV) sin W + Mbk cos 61. 

(76-4) 

The auxiliary equation for (76-4) is seen to be the same as (76-3) 
above, as it should be. 

In order to study the motion completely in a numerical case, it 
would be necessary to find all the roots of the auxiliary equation 
corresponding to (76-3). 

The characteristic equation for the system (76-1) can also be 
found by a shorter and more general method. Let us put 

^ = ae^, 6=be rt . 
Then $5 = are**, ^ 
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Substituting these into the system 



and dividing the results throughout by e ri , we get 
(I r 2 + Kr+ Wh) a - Gw b = 0, 
Gwa + ( Ar* + kr + wl) b = 0. 

These homogeneous equations in a and b will have a common 
solution if, and only if, the determinant of their coefficients is equal 
to zero, or /fj ^ 



Cur 



= 0. 



Expanding this determinant and arranging the terms in descending 
powers ofr, we get (76-3). 

In order that the motion be stable, the values of and 6 must not 
increase indefinitely with time. This means that the roots of (76-3) 
must either be pure imaginaries or that the real parts of the complex 
roots must be negative. A necessary condition for stable motion is 
that all the terms in (76-3) must be positive.* Then since Whwl 
must be positive, A and Z must have the same sign. This in turn 
means that the weight w must be placed below the trunnion axis 
(see footnote, Art. 76). 

Continuing with the solution of (76*2), we find a particular 

solution by putting 

G^Esmbt + Scosbt, (76-5) 

where E and 8 are constants to be determined. From (76-5) we get 
6 = Kb cos bt Sb sin bt, 
= - ffi* sin bt - Sb 2 cos bt, 
=:-Eb* cos bt 4- Sb* sin U, 
6^ jR6 4 sin bt + Sb* cos bt. 

* See Routh's Stability of Motion , p. 14, or Advanced Rigid Dynamics 
(6th edition, 1905), pp. 219-23. 
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Substituting into (76-2) these values of 0, 8, etc., 
equating the coefficients of sin bt and cos bt on the two of the 

equation, we get 

[Alb* - (Iwl + Kk + C*a)* + WAh) b* + Whwl] M 

+ [(Ik + AK) 6 3 - b(Kwl + kWh)] S = 0, 
[(Kwl + kWh)b- (Ik + AK) 6 3 ] E 

+ [Alb* - (Iwl + Kk -f C*o> 2 -f TM A) 6 + Whwl] S = - Cto&lf . ^ 

(76-6) 
For brevity, let us put 

E = ^4/6 4 - (Iwl +Kk + C 2 w 2 + TF 'A h) b 2 + FAwi, 



Then equations (76-6) become 



Solving these for R and S, we find 



Then (76-5) becomes 



which can be written in the form 

sin(W - a) ' (76 ' 7) 



E 
where a = tan" 1 ^ . 

Ju 

To find a particular solution of (76-4) we put 

<f>=Lsmbt+Neosbt, (76-8) 

where L and N are constants to be determined. Then 
<f>=Lb cos bt Nb sin bt, 
<j>=- Lb* sin bt - Nb* cos 6*, 
$ = - .L& 3 cos bt + Nb* sin U, 



13 
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Substituting into (78-4) these values of ^ 5 $ 9 etc., and equating 
coefficients of sin& and cos U on the two sides of the equation, we 

get 

[ A Ib* -(AWh + Kk + (Pop 4- Iwl) b z + Whwl] L 

Ah*}. 

. (76-9) 

i-(AWh+Ek + CW + Iwt)1P+Whwl\N~Mbk.) 

Replacing the bracketed quantities by E and F as before, we have 



Solving these for L and N, we get 

M(wl - Ab*) E - (Mbk) F 

^^ 



Then (76-8) becomes 

M 

{[(wl - AV) E - bkF] sin U + [(u - Ab*) F + bJcE] cos &}, 



(76-10) 



which can be written in the form 
M 



, , . 

where p = tan 

H 



Equations (76-7) and (76- 10) give the forced oscillations due to 
the pulsating action of the waves on the ship. 

77. Roll of the ship with gyro clamped 

For purposes of comparison we now investigate the rolling of the 
ship when the gyro is clamped. On putting a) = in (76-1), we get, 
since in this case $ = $:=$=0, 

(77-1) 



SHIP STABILIZES, I5 

Here the auxiliary equation is 



_ __ -K 

.. r - 



Hence 



Here the natural period of vibration is 



To find a particular solution of (77-1) we put 

<j) = (7 3 sin 6 + ^4 cos U. (77-4) 

Then 5 = (7 3 6 cos 6^ 46 sin &, 

^* = - C^ft 2 sin U - G 4 6 2 cos 6^. 

Substituting these into (77'1) and then equating coefficients of 
sin bt and cos bt on the two sides of the equation, we get 



Solving these for (7 3 and 64, we find 

M(Wh-Ib*) 



Hence (77-4) becomes 

M 



which reduces to 

9 *-*' (77 ' 5) 



Kb 

where r 



Equation (77-5) gives the forced oscillation of the ship due to the 
action of the waves. 
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78* Numerical example 

The following data, according to Perry/ 11 are approximately 
correct for a ship of 6000 tons displacement, having a metacentric 
height of 1.6 feet, a natural roll period of 14 seconds, and equipped 
with a gyro 12 feet in diameter, weighing 10 tons, and spinning at 
100 radians per second, all data being in foot-pound-second units: 

I^IO 8 , 4 = 7x10* Ow = 2.5xl0 6 , JT = 4xi0 6 , 



Substituting these data into (76*3) and then dividing through by 
j W6get * ** 0. (78-1) 



Our first problem is to determine the nature of the roots of (78-1) 
and then find them. We could proceed to find the roots directly by 
Graeffe's root-squaring method, without previously determining 
their nature, but in dynamical problems it is better to determine the 
nature of the roots before trying to find them. 

If the characteristic or auxiliary equation of a fourth-order 
differential equation be written in the form 

x*+px* + qx*i-rx + s = Q, (78-2) 

the necessary and sufficient conditions that the motion determined 
by the given differential equation shall be stable (that is, that the 
real parts of the complex roots of (78-2) shall be negative or zero) 
is that all terms in (78-2) must be positive and that 

pqr r 2 ~^?% > 0.*)* 
These conditions are fulfilled in equation (78- 1), for 

(7.18) (19.41) (1.83) - (1.83) 2 - (7.18) 2 (2) = 148.6. 
Furthermore, if a quartic equation be written in the form 

Q, (78-3) 



_ ry 

let H = ac o , 

/=ae-4& 

J = ace + 2bcd - ad* - eb* - c 3 , 



* 'The Use of Gyrostats', Nature, Lond., 12 March. 1908. 
f Routh, Stability of Motion, p. 14, -or Advanced Rigid Dynamics (6th 
edition, 1905), pp. 219-23. 
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Then the conditions that all the roots of (78-3) be complex are 

A>0 3 
and H and a 2 ! 12U 2 not both negative.* 

Applying these conditions to equation (78-1), where 0=1, 
6 = 1.795, c = 3.24, <2=0.46 5 e = 2, we find 

H = 0.02, I = 32.67, J = - 28.83, 
A = 12,430, # 2 I-12J? 2 = 32.67. 



Hence all the conditions for complex roots are satisfied, and we 
now know that all the roots of (78-1) are complex and that the real 
parts of the roots are negative. 

In order to find the roots of (78-1) we shall not apply Graeffe's 
method but shall find them with less labor by utilizing all our 
knowledge concerning them. Since the roots are all complex and 
their real parts are negative, we represent them by 

#1 ifii and a 2 %> 
or -c^ + ifii, -%-i^, -a 2 -f-ijff a , -a a -i>ff 2 . 

Then from the well-known relations between the roots and coeffi- 
cients in a rational integral equation,f we have 

(a) a 1 + a 2 = 3.59, 

(6) oj+/?f + ai+/0i + 4a 1 a 2 = 19.41, 



(78-4) 
(c) 

(d) 

We now utilize some additional knowledge concerning the roots 
of equation (78-1). We know from our previous study of the 

* Burnside and Panton, Theory of Equations, vol. 1 (1899), p. 145. 
If a l9 a 2 , a 3 , a 4 are the roots of the quartic equation 



these relations are 
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gyroscopic pendulum that /? 15 which corresponds to Jc^ of Art. 69, Is 
very small; and since the slow vibrations are not quickly damped 
out, we may infer that % is also small. We also know that /? 2 , which 
corresponds to k t of Art. 69, is fairly large; and we may also suppose 
that a 2 is also large, since the rapid oscillations are soon damped out. 
We can therefore get approximate values for the smaller roots of 
equations (78*1) by neglecting the first two terms, and we can get 
approximate values for the larger roots by neglecting the last two 
terms. 
Hence to find the approximate values of the smaller roots we have 



r= 



Neglecting the last two terms of (78- 1 ) to find approximate values 
of the larger roots, we have 

r 2 + 7.18r+ 19.41 = 0. 



Approximate values for a v a 2 , l9 /? 2 are therefore 

% = 0.047, a 2 = 3.59, ^ = 0.317, /? 2 = 2.55. 

We now proceed to find more accurate values by the method of 
iteration. Starting with a x = 0.047, we substitute this into (78-4 (a)) 

andg6t a 2 = 3.59-0.047 = 3.543. 

Then substituting the value of 4ce 1 a 2 = 4(0.047) (3.543) = 0.67 into 
78*4 (5), we have 



(of +/?!) + (o| +/?!) = 19.41 - 0.67 = 18.74. 
For brevity, put a 2 +]3* = x 9 a| +y?l = y. Then 

x+y= 18.74 

and xy = 2, 

from (78*4 (d)). From these we get x = 0.11, y = 18.18. That is, 
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Substituting these and a 2 into (78-4 (c)), we have 

3-543 (0.11)4-a 1 (18.I8) = 0.915 3 
whence o x = 0.0289. 

We now repeat the iteration by substituting this value of % into 
(78-4 (a)), thus getting a 2 = 3.56. Then 



= 4(0.0289) (3.56) = 0.41. 
Substituting this into (78-4(6)), we have 

(of +y?f) + (o| -f /?|) = 19.41 - 0.41 = 19.00. 
Then 



from which # = 0.106, y= 18.87. That is, 

of +yfff = 0.106, a|-fy?l= 18.87. 
Now substituting these into (78-4 (c)), we have 

o^lS.87) + 3-56(0.106) = 0.915, 
whence a x = 0.0285. 

We now iterate once more by substituting this value of a^ into 
(78-4 (a)), giving a 2 = 3.56. Then 4a 1 a 2 = 0.406, and (78-4(6)) 
becomes 



) = 19-41 . Oi406 = 19-004- 

Hence we have x + y = 1 9. 004, 

xy = 2 9 

from which # = 0.106, y 18.87. Then substituting these and a 2 
into (78-4 (c)), we get <%! = 0.0285 as before. Since these values can 
not be improved by further iteration, we take them as correct. 
From a| -f /?| = 0. 106 and a t = 0.0285, we get 

)fff = 0.106- (0-0285) 2 = 0.1052, 
& = 0.324. 
Also, al+/?| = 18.87 and a 2 = 3,56. Hence 

y?|=18.87-(3-56) 2 = 6 1964, 

,#2=2.49. 
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The equation for the natural vibration of the ship with gyro 
spinning is therefore 

.= e ~Mm (A cos 0.324$ + B sin 0.324$) 

+ e-* m (C cos 2.49$ + D sin 2.490, (78-5) 

where A 9 J3, C, D are constants of integration to be determined from 
initial conditions. 

The periods of oscillation are 

19 ' 4 seconds > T * = = 2 - 62 seconds. 



For purposes of comparison to determine the effect of the gyro 
on the natural oscillation period of the ship, we compute the 
oscillation period when the gyro is clamped. On substituting into 
equation (77-2) the numerical values of/, K and Wh, we find 



cos 0.447 + C 2 sin 0.447$). (78-6) 

Here the period is 

2n 
T = = 14.06 seconds. 

The gyro stabilizer thus lengthened the natural oscillation period by 
37.7 % and also caused the amplitude of the oscillations to decrease 
more rapidly. 

We have yet to examine the forced oscillations. From equation 
(76-7) we found the amplitude of 6 due to the pulsating force of the 
waves to be 



In order to indicate the effectiveness of the stabilizer in preventing 
the rolling of a ship in a rough sea, we shall consider the case of a 
ship rolling among high waves having a period near the natural 
(free) oscillation period of the ship. In this case the effect of the 
waves will be cumulative, thereby causing the ship to roll with wide 
amplitude. We take b = 0.455, thus giving a period of 

= 13. 8 seconds 



0.455 
for the forced oscillations of the ship. 
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On substituting the numerical values of the quantities 

in E and F as given in Art. 76, we find 



Then from equation (76-7) we have 





x 10x 0.465) 



113.7511 
"141.4x10** 

From equation (76-10) we have 

M 



7.223 Jf 0.051 iJf 



141.4 x 10 6 10 6 * 
Here we notice in passing that 

(Amp.) g 118.75 
(Amp.), 7.223 1O -'*' 

or (Amp. )0= 15.75 (Amp. L. 

The amplitude of the forced oscillations when the gyro is 
clamped is given by equation (77*5). Hence we have 

M 



* VflfPA-Jft 8 ) 2 - 

M _Q.513Jf 
" 1.950 xlO 6 "" 10 6 ' 

which is ten times as great as the amplitude when the stabilizer was 
acting. 

The foregoing study shows two marked effects of a gyroscopic 
stabilizer on the rolling of ships: 

(1) The stabilizer reduces the angular velocity of rolling and 
thereby increases the period of rolling, thus making the period of 
rolling greater (in most cases) than the period of the waves.* This 

* One of the longest periods of waves on record was 16-5 seconds, and this 
occurred with waves 46 feet high. See Gaillard, D JD ., Wave Action in Relation 
to Engineering Structures, Washington, 1904. 
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increases in the period of rolling prevents in large measure the 
cumulative effects of the oncoming waves. 

(2) The stabilizer makes a drastic decrease in the amplitude of 
rolling, especially when the period of the waves is nearly the same 
as the natural oscillation period of the ship. 

Fig. 60 shows graphically the marked effect of the stabilizer on 
the amplitude of rolling. 





-. Stabilizer locked 


* Stabilizer acting * 


15 




1 




10 


A 


A 


n 


A 


5 


11 


fl 


A 


v A 














U ^\ r r> S\ f\ ^ 




5 


II 


/ u 


IIU 




\J* L> V ^^^ ^^\/ *^ \I\J \/ ^*s 




U 


I 


v U 




10 


~ 








15 


- 


^ 





Fig. 60 

79* Speny or active-type stabilizer 

Water waves do not beat against a floating ship as they do against 
a solid sea wall, because of the fact that a ship is free to rise and fall 
and thereby permit the bulk of a wave to pass under it. Hence any 
single wave imparts only a small amount of roll to the ship. But 
when each of a continuous train of waves imparts a small increment 
of roll, a large roll is soon built up if the period of the waves happens 
to be nearly the same as the natural period of roll of the ship. 

The idea of the Sperry stabilizer is to quench each increment of 
roU as soon as it occurs, thus preventing the building up of a large 
roll. This is accomplished by causing a forced precession of the 
gyro about the trunnion axis in the proper direction to oppose 
the roll. 

The Sperry stabilizer consists of a large gyro with spin axis 
normally vertical, a small control gyro with spin axis horizontal and 
athwart the ship, a precession motor, and a magnetic brake, plus 
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the motors required for spinning the gyros. The of the 

control gyro is free to precess about a vertical axis. Horizontal 
pins PI and p% project from the casing of the control gyro, as 
indicated in Fig. 61. Centralizing springs are attached to pin jv 
Pin p-L is free to swing back and forth between electric contacts 
c and 




Fig. 61 

The operation of the Sperry stabilizer is as follows: When the 
ship rolls to starboard, for example, by as much as a single degree, 
the control gyro processes about a vertical axis as indicated in the 
figure. This precession causes pin p^ to make electrical contact at 
c lr thereby starting the precession motor and causing it to exert a 
torque about the trunnion axis in a clockwise direction as viewed 
from the starboard side. When the ship rolls back in the opposite 
direction, the casing of the control gyro processes about the 
vertical axis in the opposite direction, causing pin p l to make 
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electrical contact at c 2 soon after the ship has rolled past its vertical 
position. This starts the precession motor in the opposite direction, 
causing it to process the main gyro in the reversed direction so as to 
oppose the roll again. 

The Sperry stabilizer also retards the rolling in other ways. If 
the rolling of the ship is fast enough to overdrive the motor, it 
thereby converts the motor into a generator for the time being. 
Since a part of the kinetic energy of rolling is expended in driving 
the generator, the rolling is reduced to some extent. 

Furthermore, as the current taken by the motor diminishes, the 
magnetic brake comes into action and reduces the roll. 

In order to obtain the differential equations of motion of the 
ship under the action of the Sperry stabilizer, let us suppose that the 
main or stabilizing gyro is forced to precess about the trunnion 
axis in a clockwise direction as viewed from the starboard side 
(see Fig. 61). Then since the precession vector is directed toward the 
port side, the rule for the direction of the gyroscopic reaction 
moment shows that this moment tends to roll the ship toward the 
port side. Hence if the ship begins to roll to starboard, this rolling 
can be counteracted by a forced precession of the main gyro in the 
direction of negative 6. It is to be noted that this is also the direction 
of the G.B.M. induced by the rolling. 

Since the main gyro and the control gyro are both nonpendulous 
in the Sperry stabilizer (pendulosity is not needed), and since no 
brake friction moment is constantly applied about the trunnion 
axis, we can obtain the differential equations for the Sperry 
stabilizer by putting w = and Jc = in equations (76- 1 ) and writing 
in a term for the forced precession moment about the trunnion 
axis. Hence the differential equations of motion are 

1$ +K<j>- Ca>0 + Wh(f> = M sin bt,} 

(79*1) 



where N denotes the variable moment about the trunnion axis. 

Because of the slight gaps with respect to time in the train of 
operations of the Sperry stabilizer, the rolling of a ship under the 
action of the stabilizer probably cannot be represented exactly by 
a continuous relation between <j) and t. However, since the action 
of the stabilizer is controlled by the rolling of the ship, we may 
reasonably assume that the applied precessional moment about the 
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trunnion axis is proportional to the angular velocity of roli of the 
ship. The applied precessional moment will then supplement the 
gyroscopic moment opposing the ship's rolling. The value of JV 
in this case would be Q<}>, where Q denotes a positive constant. 

This assumed relationship leads to the result that by making Cw 
sufficiently large the amplitude of rolling can be kept as small as 
desired while the ship is following a straight course or a course 
which curves in the same direction as the spin of the gyro; but the 
result also shows that the ship will be unstable if its course curves 
in the direction opposite to the direction of spin of the gyro. In 
view of the fact that ships are stable under the action of the Sperry 
stabilizer when the course curves in any direction, the assumption 
that N = Q<j> is not valid. 

Since a roll in the direction of positive can be quenched by a 
forced precession of the gyro spin axis in the direction of negative 
6, we may reasonably assume that 6= m<j>, where m denotes a 
positive constant. This assumption leads to the result that the ship 
will be stable when its course curves in any direction. In this case 
the precessional moment about the trunnion axis is assumed to be 
a function of t and will be denoted by E(t). It can be found, if 
desired, after ^ has been found. 

With the assumptions stated above, equations (79-1) are replaced 
by the following: 



(a) I<f> + K<j> - Cud + Wh<fi = M sin bt, 
(6) 



(79-2) 

(<0 
On substituting (c) into (a), we get 

I<j) + (J? + O^m) $$ + WFA^ = M sin W. (79*3) 

Here the characteristic equation is 



-(K + Cum) 

whence r= - ' 



If (K + G(j)mf ^ 41 W h, both values of r are real and negative, and 
there is no natural oscillation. The motion is overdamped. 
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If (K + Cwmf < 4IWh, r becomes 

_ - (JT + Cam) . ^{4JWh - (JT + Cam)*} 
r _ 2j * 27 ' * 

The free oscillations are then represented by the equation 

, r ( + Cam)n t . ^ ^. 

p = exp ~ -j ......... ............. I (^sinA^-f c 2 cosA^), 

where x J^Wk- ( K + C^}_ 

2jt 

Note that the presence of Cto> in the exponent ofe makes the damping 
much more rapid. 

To find a particular solution of (79-3), we put 

$5 = PI sin U + P 2 cos 6*. 
Then $5 = P l b cos U-P^b sin Ztf , 



Substituting these into (79-3) and equating the coefficients of 
sin U and cos "bt on the two sides of the equation, we get 



Solving for P and P 2 , we find 



(Wh- 



A particular solution is therefore 
M 



x[(Wh-~ Ib*) sin bt-b(K + Own) cos 6*], 
which can be written as 

. M 

*=1%WK^)*+IW ( 

. 
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Here it is evident that the amplitude of $ is less it be if 

a) = 5 and that it can be decreased at pleasure by Cm 

sufficiently. It is also to be noted that when to = 3 equation (79-4) 
reduces to (77-5). 

We can now find R(t), if desired, by means of (6) and (c) of (79-2). 
From (79*4) we have 



^ 

Then from (c), 



/} mbM 



Substituting into (79-2 (&)) these values of 8 and ^, we get 



x Jyimb sin (bt a) + Co) cos (bt a)] 



where 



- 1 r. 
Amb 



The precessional moment applied by the stabilizer is therefore 
periodic with the same period as ^ but differing in phase by the 
angle /?. Note that ft is likely to be near 90, because its tangent 
is directly proportional to CM. 

Because of the fact that the Sperry stabilizer quenches the 
increments of roll as fast as they occur, a smaller gyro suffices for 
stabilizing a ship than is the case with the brake-type stabilizer. 
Throughout a cruise of 29,000 miles by a yacht of 3000 tons 
displacement, a Sperry stabilizer kept the amplitude of roll down 
to 1 (2 from one extreme to the other) by means of a rotor 8 feet 
in diameter, weighing 22 tons (long tons) and spinning at 1300 r.p.ni. 

An active-type stabilizer can not only quench the rolling of a 
moving ship, but it can also produce rolling of a ship at rest. By a 
change in the electrical connections at the switchboard, the stabilizer 
can be made to apply to the ship a periodic torque that will produce 
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rolling with increasing amplitude. If, for example, a ship equipped 
with such a stabilizer should run aground, the stabilizer might 
be able to loosen It from the ground and enable it to keep going. 
Fig. 62 shows how effectively a Sperry stabilizer controlled the 

rolling of a small ship on one occasion. 



Stab, on Stabilizer off t Stabilizer on 

'26! 




Maximum roll : 2 (total) 

7,46 p.m. Rotor speed 1050 r.p.m.. Time 10 sec. {* *f 



RECORD TAKEN IN A SHIP OF 2000 TONS DISPLACEMENT 

Unsiabilised Stabilized 

Maximum roll, side to side, 28 Maximum roll, side to side, 2 

Average roll, 13 Average roll, 1 
Fig. 62 



B. SHIP ON A CURVED PATH 

80* Brake-type stabilizer 

Let us assume that a ship equipped with a gyroscopic stabilizer 
is making a turn In a circular path. Let the ship have a constant 
angular velocity ijr and a constant linear velocity V. Then from the 
relations ^2 

and a = -, 



where H denotes the radius of the circular path and a n the centri- 
petal acceleration of the ship, we find (on eliminating B) 



Then the magnitude of the centrifugal force acting at the center of 
gravity of the ship is w 



and the moment of this force about the longitudinal axis of roll is 



where p denotes the distance from the center of gravity to the axis 
of roll. Since the axis of roll is above the e.g., this moment will tend 
to heel the ship inward toward the center of curvature. 



SHIP STABILIZER 209 

Since the ship is moving on a curve, the axes of the gyro and 
frame are turning in azimuth with angular velocity i/r. We shall 
assume that the gyro is spinning in a counterclockwise direction as 
viewed from above, and that the ship is moving in the same 
direction (swerving to the left) on the curve. Hence to derive the 
differential equations of motion we proceed as in Art. 75. 

Referring again to Figs. 57 and 58, let the longitudinal axis of the 
ship lie in the direction of OX. Then the trunnion axis will lie in the 
direction of OF. Then as in Art. 75 the projections of ft on the axes 
of the moving trihedron are 

i/r cos <j> sin 6 on OX', 

ijfsm<j) on Y', 
t]f cos <l> cos d on OZ'. 

In the disturbed position of the gyro its angular velocities about 

the moving axes are 

<j> cos 6 about OX', 

6 about 07', 

about OZ f . 



Hence the total angular velocities about the moving axes are 
o) x = $GOs6 T/fG 
a)y = 6 T/r sirup, 



The corresponding angular momenta are therefore 
H x = 1($ cos 8 + ^Jr cos <fi sin 8), 



H z = G((j) <f>sixid + i/r cos ^ cos 6). 

Hence _^= /^ approx., 

d/t 

dH v .* 
-j2=A6, approx., 

dH, A 



after neglecting negligible products and replacing sines by angles 
and cosines by unity. 

14 SG 
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On substituting the above quantities into (16-4), we get 



(a) If - Gu& + (Ov& + Wh) $ + K$= JL + Jf sin bt 

9 



(80-1) 

Differentiating (6), we have 



(c) .4$+ Cfo$ -f (Cto ^ -f wZ) ^ -f M = 0. 

Solving (a) for $ and then differentiating the result twice, we have 



^ 



Substituting these into (c) and rearranging slightly, we get 



Wfi* + Ca)i]f( Wh + wl) + Whwl] <f> 

vWVilr 
= -- - -(Ca)fr + wl) + (C6)fi+wl) M sinbt 

y 

-f kbM cos U - AMb 2 sin bt. (80*2) 

Here the characteristic equation is 



= 0. (80-3) 

Since all coefficients and the constant term are positive, the first 
condition for stability is fulfilled. 
To find a particular solution of (80-2), we put 
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Then <j> = Qb cos U - Kb sin U 9 

$ = ~ Qb* sin bt - Rb* cos 6*, 
$ = - QW cos & + EW sin &*, 
0i v = Qb* sin & + -K6 4 cos bt. 

Substituting these into (80-2) and equating coefficients of like 
functions on the two sides of the equation, we get 

-pWVft(wl + Cui/r) 
~~ g [CWfi* + Ca)fi( Wh + wl) + Whwl] * 
and 



where 



)^( Wh + wl) + Whwl, 
F = (AK + Ik) & 3 - [Co)ft(K + k) + Kwl + JcWh] b, 



Solving equations (80*4) for Q and R, we find 

M(ES-kbF) M(JcbE+F8) 



A particular solution of (80-2) is therefore 

M 



x {(ES - Tcb F) sin 6* + ( JcbE + FS) cos bt] 



ft 2 + G^( Wh + wl) + Whwl] 

Z.2/j2\ 

), (80-5) 



. .JcbE+FS 
-i- 



, 
where 

The negative sign of the constant term shows that the ship is 
inclined inward toward the center of curvature of the path, and the 

14-2 
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result farther shows that the ship oscillates about the inclined 
position. 

If the path of the ship curves in the direction opposite the 
direction of spin of the gyro, the sign of ijr must be changed. Then 
(80*3) becomes 

All* + (AK + Ik) r 3 + [CPoP - 



2 - Ca)ijr( Wh + wl) + Whwl] = 0. (80-6) 
In this case the first condition for stability requires that 
0*>$(A + 1) <(CW + AWh+Iwl+KJc), 

C<*)TJf(K +k)< (Kwl + k Wh) 
and Ca)i/r( Wh + wl) < (O 2 ^ 2 + Whwl). 

Changing the sign of ^ also causes a change of sign in some of the 
terms in E, F, 8 and P 3 and therefore affects the angle of heel of the 
ship and the amplitude of the roll. The behavior of a ship is thus 
less certain when its path is curving in the direction opposite the 
direction of spin of the gyro. 

81. Active-type stabilizer 

To find the differential equations of motion under the action of 
the Sperry stabilizer, we put w = and 4 = in equations (80-1) 
and then make the assumptions of Art. 79. The equations are then 

(a) 

(b) 

(c) 6= -m<}>. 

On substituting (c) into (a) we get 

)$ = ^ + M&ibt. 

(81-2) 
eristic equation is 

r* + (K+ Com r+Wh+ Ca>ir = 
whose roots are 



Here the characteristic equation is 

Ir* + (K+ Com) r+Wh+ Ca>i/r = 0, 



- -4J(Wh + C 

r _ . _ 
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If (K + Cwm)* > 4/( Wh + Cwifr), both roots mil be real and negative. 
The motion will be overdamped and there will be no natural 
oscillation. 

If (K + Cum)* < 4/( Wh + Cufi), the values of r will be 

own) 2 } 



21 21 

In this case the free motion of the ship is a damped roll of period 

47TI 

""" V{4/( Wh + Go)i/r) -(K + Owm) 2 } ' 

which could be very long. 

If the ship is turning in the direction opposite the direction of spin 
of the gyro, i/r will be changed in sign. The characteristic equation 
in this case becomes 



If the constant term ( Wh Ctoijr) is negative, that is if Co)i/r > Wh, 
both values of r will be real and one of them will be positive. Hence 
the motion will be unstable in this case. 

To find a particular solution of (81-2), we put 

(p == P l + Qi sin bt + Q 2 cos &* 
Then = 



<j)=~-Qi b* sin bt - Q 2 6 2 cos 6*. 

Substituting these into (81-2) and then equating coefficients of like 
functions on the two sides of the equation, we get 

-pWVifr 

1 



b(K + Com) Qi + ( W h + Cuft - J6 a ) Q 2 = 0. 
Solving for Q^ and Q 2 , we find 



coi/r - 



'(Wh + Ca)i/r - /ft 2 ) 2 
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A particular solution of (81-2) Is therefore 



"" g( Wh + Oo>^) 

x [( Wh + Cto^ 1& 2 ) sin bt-b(K + Cum) cos 6q, 

which, can be written in the form 

M _ . b 
b*)* + b*(K+ 7) ' 

(81-3) 

, , 
where y = tan" 1 

If the ship is turning in a direction opposite the direction of spin 
of the gyro, the sign of i/r must be changed. Equation (81-3) then 
becomes 

M . 

)j 



^ 

(81-4) 



_ 
where 



It will be seen that in this case the ship is heeled inward toward the 
center of curvature of the path and at a greater angle of inclination 
than when it is turning in the same direction as the spin of the gyro. 
The amplitude and phase angle of the roll have also changed. 

If we put J/T = in (81-3) and (81*4), we get (79-4). The amplitude 
and phase of roll on a curved path, aside from the heel, are thus 
not the same as on a straight course. 

The unsymmetrical effect of stabilizers on curves could be 
corrected by using two identical gyros rigidly connected by 
toothed gears somewhat in the manner shown in Fig. 65. If both 
gyros rotated in opposite directions at the same speed and precessed 
in opposite directions at the same speed, the one-sided effect of 
either gyro when traversing a curve would be annulled by the other 
gyro. 



CHAPTER IX 

The Gyroscope as a Stabilizer (Continued) 
Monorail Cars 

A. MONORAIL CABS ON A STRAIGHT TRACK 

Railway cars designed to run on a single rail, with the center of 
gravity above the rail and held in dynamical equilibrium by 
gyroscopic action, were invented during the first decade of the 
present century. The first car of this kind was invented by Louis 
Brennan of England about 1905, and further developed and 
improved by him in subsequent years. On this car the gyro was 
mounted with its spin axis horizontal and at right angles to the 
track (rail). 

In 1909 August Scherl of Germany brought out a monorail car on 
which the gyro was mounted with its spin axis vertical, as in the 
case of the ship stabilizer. The Scherl type of car was studied and 
possibly improved a little later by P. P. Schilowsky of Eussia.* 
More recently (1952) a monorail car of improved design has been 
demonstrated in Germany. 

82. Monorail car with gyro axis vertical (Scherl and Schilowsky 



The type with vertical spin axis is represented schematically in 
Fig. 63. This type of car has a weight mounted above the trunnion 
axis so as to make the center of gravity of gyro and frame lie above 
the trunnion axis, thus making the combined gyro and frame 
gravitationally unstable about the trunnion axis. The car and 
mounted gyroscope thus constitute an inverted gyroscopic pendu- 

* ' The present writer had a considerable advantage over both Brennan 
and Scherl, in that he entered the field after them and was therefore in a 
position to observe and avoid all the deficiencies of their systems.' 
P. P. Schilowsky, The Gyroscope: Its Practical Construction and Application 
(1924), p. 137. 
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lum of the type treated in Art. 74. To find the differential equations 
of motion, let 

W = weight of car and all equipment, 
w = weight of gyro and frame, 
h = height of e.g. of car and equipment above rail, 
I = height of e.g. of gyro and frame above trunnion axis, 
1= moment of inertia of car and equipment about rail, 
A = moment of inertia of gyro and frame about trunnion axis. 
Then 

Wh sin <j> = Wh(j) = moment of W about rail, 
wl sin d = wld = moment of w about trunnion axis, 
K$ = moment of friction about rail, 
W = damping moment about trunnion axis. 

Hence with axes as in Fig. 56 we have 

a) x , = ^ cos 9 = <j> approximately, a)^ =6, (L>^ = ^ sin 6 ; 

H x , = 1$ cos 6 = I<j> , practically ; 

H yf = Ad, H^ = C(o) <j> sin 6) = Co), practically. 

-- 



Substituting the above quantities into (16-4) and keeping in mind 
the directions of the moments, we have 



Since ^, 6 and sin d are all small, we neglect the products containing 
d<j> and put cos 6 = 1 . Then we have 



(a) 
(b) 
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These equations can also be derived by the G.R.M. method as 
follows:* 

A forced precession of angular velocity <j> about the rail induces a 
G.R.M. of magnitude Co)(j) about the trunnion axis, and this G.R.M. is 
negative (its vector points in the negative direction of the t/-axis). 




Fig. 63 

Then the forced precession about the transverse axis in the 
negative direction induces a G.B.M. of Oo>( - 6} about the rail axis in 
the backward (negative in this case) direction. That is, the G.B.M. 
in the direction of the rail is - Ca>( - 6) = Cud. Then by the funda- 
mental equation of rotation about an axis (Art. 20) we have: 



or 



which are the same as (824). 

* Note that the forward direction of motion is in the negative direction of 
the ic-axis and that, as increases, the velocity vector # points in the 
negative OS-direction. 
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To find the characteristic (auxiliary) equation for the system 
(82-1), we assume ^ = ^ ^^ 

Then ^ = ar e^, 9 = bre*, 



Substituting these into (82- 1 ) and then dividing the results through- 
out by G H , we get 



G(t)ra + (Ar 2 + Jcr wl)b = 0. 

These homogeneous equations in a and b will have a common 
solution if, and only if, the determinant of their coefficients is equal 
to zero. Hence we have 



-Cur 

= 0. 
Co)r 



Expanding and arranging in descending powers of r, we get 
A Ir* + (K A + M) r 3 4- (<7 2 co 2 - A Wh - Iwl + KTc) r 2 

l = Q. (82-2) 



A first condition for the stability of the monorail car is that all the 
coefficients and the constant term in (82-2) shall be positive. This 
means, first of all, that I and h must have like signs, which means 
that the center of gravity of the gyro and frame must be above the 
trunnion axis. 

Since the moment of rail friction, K<j>> is necessarily a retarding 
moment, the coefficient K is necessarily positive and the sign of the 
product wlK cannot be changed. In order that the coefficient of r 
be positive, the term Whk must be changed in sign and must be 
numerically greater than wlK. This means that the moment kd, 
which we assumed to be a retarding moment, must actually be an 
accelerating moment. Hence if the monorail car of the Scherl type 
is to be stable, there must be some provision for accelerating the 
precession about the trunnion axis. Scherl surmounted this 
difficulty by installing a servomotor to apply automatically an 
accelerating moment about that axis, as did Sperry in the case of his 
ship stabilizer. We therefore change the differential equations (82-1) 
so as to provide for the accelerating moment. 
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Let us assume that the servomotor applies an accelerating 
moment of magnitude R<j) about the trunnion axis, where R 
denotes a positive constant. Then since there is no braking moment 
led (we neglect bearing friction), equations (82-1) become 

'LEA } (82 ' S) 

The characteristic equation for (82-3) is found as above by putting 

(j> = ae^y 6 = 66^. 
The result is 

A Ir* + AKr* + (CW -AWh- Iwl) r 2 + (Co)R - Kwl) r + Whwl = 0, 

(82-4) 
These coefficients will all be positive if 



and Co)R>Kwl. 

These conditions can evidently be satisfied by increasing Ca> 
sufficiently. 

Because of the fact that the gyro in a monorail car must spin at 
very high speed (higher than in the case of the ship stabilizer), we 
may treat the quantity A Wh 4- Iwl as negligible in comparison with 
<7 2 o> 2 . On neglecting this quantity, equation (82-4) becomes 

AIr*+AKi* + C*w*r* + (C(>R-Kwl)r + Whwl^O. (82-5) 

We know from our previous studies of the gyroscope that two 
roots of (82-5) are small and that the other two roots are con- 
siderably larger. We get approximate values of the larger roots by 
neglecting the last two terms of (82-5) and solving the quadratic 
equation 



from which we get 

K 

v ~ ~ _ -4- ti 

r ~ - 



i 

2AI 



if A Z K Z < 4AIC 2 (t) 2 . We may further neglect the second term under 
the radical in comparison with the first term. Then we have 

K . Cu 
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Here the period of oscillation is 



which may be very short. 

To get the two smaller roots of (82-5), we neglect the first two 

terms and solve the quadratic 

<7 8 a a r 2 + (CcuR - Kwl) r + Whwl = 0, 

which gives 



- Kwl) . ^^Whid- (Ca)R - 

^ i 



if 

If we neglect the quantity (Ca)R-Kwl)* in comparison with 
4:0^ Whwl, the value off becomes 

_-(Ga)R-~Kwl) -VO^M 
r ~" 2(7 2 a> 2 l Co) ' 

In this solution the damping coefficient is small and the period 



is long. 

It is to be noted that if Co) R > Kwl, the real parts of all the roots 
of (82-5) are negative and the motion of the monorail car is therefore 
stable. 

83. Monorail car with gyro axis horizontal (Bremian type) 

The Brennan monorail car is represented schematically in Fig. 64. 
The gyro is mounted with its spin axis horizontal and perpendicular 
to the rail, as previously stated. The frame is made unstable about 
the precession axis by application of an unstabilizing moment 
proportional to d. This unstabilizing moment is applied auto- 
matically by a mechanical device and corresponds to wW in Art. 82, 
the coefficient of 6 having the same sign as WJi. In the following 
treatment of the Brennan car we shall use the same notation as in 
Art. 82, even to the point of using wl for the coefficient of the 
unstabilizing moment acting on the frame. We shall derive the 
differential equations of motion by the G.B.M. method (see footnote 
on p. 217). 
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Fig. 64 

Referring to Fig. 64, we see that if the ear tilts to the right with 
angular velocity <j>, this is a forced precession about the rail axis. 
Hence this forced precession induces a G.B.M. of Caxj> about the 
vertical axis OZ and in a counterclockwise (negative) direction as 
viewed from above, or COMJ>. 

The forced precession about OZ in the negative direction induces 
a G.B.M. of Cu(-6) about a line parallel to the rail and directed 
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backward, that is, in the negative direction. Hence the G.E.M. in 
this case is <7et>( $), or Cw6. Then from the fundamental 

equation of rotation about an axis (Art. 20) we have 



and 

Ad = wl0 + k& Co)(j> (here we take Jed to be an accelerating 

moment), 

or (a) If +K(j>- C(*>6 = TPfyM 

(b) Ad 



On putting fi ae 1 *, 6 be 9 * and proceeding as in Art. 82, we find 
the characteristic equation for the system (83*1) to be 



lz=Q. (83-2) 
All coefficients in (83-2) will be positive if 

AK>Ik, Whk>wlK, 
and (7V 2 > Kk + A Wh + Iwl. 

Since I> A, the above conditions lead to the further conditions 

. wl A 
and < . 



Because of the fact that the gyro in a monorail car must spin with 
very high speed (higher than in the case of the ship stabilizer), we 
may assume that the sum of terms Kk+AWh+ Iwl is negligible in 
comparison with (7 2 6> 2 . When those terms are neglected, equation 
(83-2) becomes 

Air* + (AK - Ik) r 3 + <7 2 o>%- 2 + ( Whk - wlK) r + Whwl = 0. 

(83-3) 

We know from our previous studies of the gyroscope that two 
roots of (83-3) are small and that the other two are considerably 
larger. We get approximate values of the larger roots by neglecting 
the last two terms of (83-3) and solving the quadratic equation 

Air* + (AK - /i) r + <7 2 w 2 = 0, 
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from which we get 

- (AK-Ik) 



~ 2AI 

Now treating the quantity (AK-Ilc) z as negligible in. comparison 
with 44/CV, we get 



To get approximate values of the smaller roots of (83*3) we neglect 
the first two terms and solve the quadratic 



+ ( WKk - wlK) r + Whwl = 0. 
From this we get 



We now neglect the quantity (Whk-wlK)* in comparison with 
4:C*a) 2 Whwl and thus obtain 



The values of r given by (83-4) and (83-5) show that if the real 
parts of the complex roots are to be negative, we must have 

AK>Ik from (83-4), 
and Whk>wlK from (83-5), 

as previously stated on the basis of Eouth's conditions. 
The periods of the two types of oscillations are 

27rCa) and r- 
2 ~ 



Hence the long period (2\) can be lengthened by increasing a) or by 
decreasing I; and the short period can be shortened by increasing 
o) or by decreasing A. 
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B. MONORAIL CARS ON CURVES 

84. Vertical axis type 

Since the gyro on a monorail car of vertical axis (of spin) type is 
mounted in the same manner as the gyro of a ship stabilizer, the 
differential equations of motion of such a car on a curve will be very 
similar to those for a stabilized ship on a curve. We find them by 
substituting the appropriate quantities into equations (16-4). The 
left-hand members will be the same as for the ship stabilizer in 
Art. 80. To find the moments M x and My, let E<}> denote the 
accelerating moment about the trunnion axis, where R denotes 
a positive constant. Then the moment about the precession axis will 

be 



Since the centrifugal force tends to tilt the car outward away from 
the center of curvature, the moment about the re-axis will be 



The differential equations of motion for the Scherl monorail car 
running on a curve are therefore 



g 

=zwld-M<f>, 
or 

(a) 

" ' (84-1) 

(6) Ad+Ca)<t> + (Ca>i/r-wl)6+It<f>=:Q. 
To solve these equations for $5, we first differentiate (6) and get 

Solving (a) for 6 and then differentiating the result, we have 



- - Wh) fa 

CCt) 

8=* r [/0 lv + K$ + (Co)T]f - Wk) $\ . 
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Substituting these into (c), we get 
AI$* + AK$ + [<7 2 w 2 + Ca)i/r(A +I)-(A Wh + ltd)] $ 

I - Kwl) < + [C 2 ^ 2 - Ca)i/r( Wh + wl) 

wl). (34-2) 



u 

Here the characteristic equation is 

Air* + AKr* + [C%; 2 + C<*>$( A + I)~(AWh + Iwl)] r* 

$ + Ca)B - Kwl) r 

- Co)fi( Wh + wl) + Whwl] - 0. (84*3) 
For stability, all coefficients must be positive and we must therefore 

)>(AWh + Iwl), 



CWft* + Whwl > Co)i/r( Wh + wl). 

These conditions can be met by making Co) sufficiently large. 
To find a particular solution of (84-2) and thus the amplitude of 

oscillation, we put , ~ , , 

r <p = Q } a constant. 

Then ^ = ^ = ^=^ = 0. 

Substituting these into (84*2), we get 



! ^ 2 - C<t)ijr( Wh + wl) + Whwl] " 

We have already found that the denominator must be positive if 
the car is to be stable. Hence (j> will be negative only when 



This means that the car will lean inward toward the center of 
curvature only when Cd)i/r < wl. 

A necessary condition for the dynamical equilibrium of the car is 
that the resultant of the centrifugal force and the weight of the car 
must pass through the rail. The requirement that Coiijr < wl puts a 
severe limitation on the velocity of spin of the gyro and may be 
sufficient to explain the failure of the Scherl car to run satisfactorily 
15 SG 
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on curves. The defect was corrected by using two gyros of the same 
size, spinning and processing at the same speeds but in opposite 
directions, and having their frames rigidly connected by toothed 
gears (see Fig. 65). 

If the car is running on a track which curves in the direction 
opposite the direction of spin of the gyro, the sign of i/r must be 
changed in the preceding equations. When the sign of fi is changed, 
the characteristic equation (84-3) becomes 

A IT * + AKr* + {CW - [Co)fi(A +I)+AWh + Iwl]} r 2 
+ [C&R - (Ca)Ki/r + Kwl)] r 
+ [O 2 ^ 2 + Ca)$( Wh + wl) + Whwl] = 0. 
Here the conditions for stability are 

0*a>* > [Co)fi(A +I) + AWh+ Iwl], 



which can be met by making Co) and R large enough. 
The angle of inclination in this case becomes 



)fi( Wh + wl) + Whwl] 

and is definitely positive. The car therefore always leans toward 
the centre of curvature of the track, making it always possible for 
the resultant of the centrifugal force and the weight of the car to 
pass through the rail and thus fulfill the fundamental condition 
for dynamical equilibrium. 

The Schilowsky car was forced to precess by means of a pen- 
dulum. 

85. Horizontal axis type 

When the spin axis of a gyro is horizontal, a rotation in azimuth 
augments or diminishes its precession about a vertical axis, 
depending on the direction of the rotation. Hence when a monorail 
car of the Brennan type runs on a curved track, the rate of preces- 
sion 6 is increased or decreased by the angular velocity with which 
the car changes its direction. Furthermore, the precession angle 6 
is also increased or decreased by the angle fit. If the track is 
curving to the right (see Fig. 64), the angular velocity of precession 
will therefore be 6 + fi and the precession angle after time t will be 
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Fig. 65 
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Q + fa. We can therefore write the differential equations of motion 
for the Brennan car on a curve by replacing 6 by 8 + ^ and replacing 
by 6 + fa in equations (83-1), and then adding to the right-hand 
side of (83- 1 (a)) the moment of the centrifugal force. The differential 
equations of motion are therefore 



(d) ~ T , ^ ^ , j- , .. .- T y. F)*l\ 

(b) Ad + k(6 4- ^) -f C<s)<j> - wl(6 + fa) = 0. 

The term on the right in (85*1 (a)) is negative because the moment 
of the centrifugal force tends to overturn the car outward (to the 
left) and thus to decrease 0. 

To solve equations (85*1) we first differentiate (6) with respect to 
time, thereby obtaining 

(c) AS+ Jed + Ca)$ - wl(6 + i/r) = 0. 
Then from (a) we have 



and 



On substituting these into (c), we get 

2 w 2 +JKk-AWh-Iwl) $ 



- ( Whk + wlK) <$> + Wkw\4 = . ( 85 . 2 ) 

a 

To find a particular solution of this equation, we put 



Then ^ = ^ = ^=^=0. 

On substituting these into (85-2), we get 

Vilf 
Q = = -Jtl. (85-3) 

a 

The car thus inclines inward toward the center of curvature. 
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To find the corresponding particular solution for 0, we substitute 
the above particular value of $ into (85-1 (a)), giving 



integration of which gives 

0=-jk 

The precession of the frame about the vertical axis is therefore not 
oscillatory but increases with time. Hence the car would soon 
overturn if left to itself. 

If the car is running on a curve that turns to the left, we find its 
behavior by changing the sign of f. Then (85-3) becomes 



which shows that the car is tilted toward the center of curvature in 
this case also* 

In order to make his monorail car traverse curves either to the 
right or to the left with equal facility and remain dynamically 
stable while doing so, Brennan equipped Ms car with two exactly 
similar gyros which rotated at the same speed and in opposite 
directions. The gyros were mounted side by side with their spin axes 
in the same straight line when the car was in equilibrium, as indi- 
cated in Fig. 65. The frames were coupled by means of spur-toothed 
gears so that they precessed at the same rate but in opposite 
directions. Whatever effect the azimuthal motion produced on one 
gyro was therefore annulled by the other. The car always leaned 
toward the center of curvature, however, because the centrifugal 
force acted on the whole car and equipment. 



CHAPTER X 
Astronomical Applications 

86* Precession and nutation of the earth's axis 

The earth, may be regarded as an enormous gyroscope whose 
center of mass is constrained to move in an elliptic orbit around the 
sun. The form of the earth is an oblatejpheroid whose equatorial 
diameter is about 27 miles greater than its polar diameter. The 
polar axis (axis of spin) is not perpendicular to the plane of the 
earth's orbit (plane of the ecliptic) but is inclined from the per- 
pendicular by about 23^. The direction of the polar axis in space 
changes but very little (about 50" A a year) as the earth proceeds in 
its orbit around the sun. 

Because of the equatorial bulge and the inclination of the polar 
axis just mentioned, the resultant attraction of the sun or moon 
does not pass through the center of mass of the earth. The resultant 
attraction is therefore equivalent to a force acting through the 
center of mass and a couple which tends to straighten up the polar 
axis and make it perpendicular to the plane of the ecliptic (see 
Fig. 66). Since, however, the spinning earth is free to turn in any 
manner about its center of mass, this couple can do nothing but 
cause the spin axis to precess (with a slight nutation) around the 
perpendicular to the plane of the ecliptic, the type of motion, 
studied in Chapter IV. 

By reason of the fact that the moon is comparatively near to the 
earth (238,800 miles away) and the plane of its orbit is inclined only 
about 5 from the plane of the ecliptic, the precession due to the 
moon's attraction on the equatorial bulge is more than twice as 
great as that due to the sun. 

To find a mathematical expression for the couple causing the 
precession and thence the differential equations for the precession 
and nutation of the polar axis, let us assume that the sun is the 
attracting body. Assume that the earth is at the winter solstice, so 
that the precessional couple has its maximum value. Take a set of 
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rectangular coordinate axes with, the origin at the center of mass 
of the earth and with the i-axis along the line of intersection of the 
plane of the ecliptic by the plane of the earth's equator. Then a 
plane through the polar axis of the earth and the center of the sun 
will be perpendicular to the &-axis 3 and the line in which this plane 
intersects the plane of the equator may therefore be taken as the 
y-axis (see Fig, 67). 

Let (x, y, z) be the position of any particle dm in the earth, and 
let p denote the position vector of the particle. Let R be the 
position rector of the sun and let r be the vector drawn from 
(x, y, z) to the sun, so that R = p + r (86-1) 

Let M denote the mass of the sun and y the constant of gravita- 
tion. Then by the law of universal gravitation the particle dm is 
attracted along the vector r with the force 

jMdm 
ajf -~~^ > 

or, as a vector, 
. 

CJ,JP 



where c denotes a unit vector along r. 

Now from (86*1) 

r = jtv p. 

And from Kg. 67 p = xi + yl + zk, 



Hence r = jJ? cos 6 - kB sin 6 - ix - jy kz 

y) k(J?sin# + z). 



Then dF = { _ ix+j(E cos 6 - y) - k(jR sin 6 + z)}. 



The moment of this force about the origin is, by (12*2), 

i j k 

x y z 

-x Mcos6-y -(Rsind + z) 



, , yMdm 

dM = p x dF = - g 



= \R(y sin ^ + z cos 5) + jJJa; sin 6 + kEx cos 6, 
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Fig. 66 




Fig. 67 
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jRMdm . 

L 
r* 



vMMdm 
* L - g 



(86-2) 



Since r is a function of x, y> z through (86-1), equations (86*2) 
cannot be integrated until r is expressed in terms of these variables. 
By the law of cosines we have 

- 2Rp cos (R 9 p) 



- 2(jS cos 6 - kB sin 6) * (ix -f jy + fcg) 
) ) by (5-6) 



--2fsin5) j, 
since (p/R) 2 is negligible in comparison with 1.* Hence 



~ 



and 



1 1 f 2 1 

^=-^1 l-^(ycos^~^sin^) 



Now expanding the bracketed quantity into a binomial series 
and retaining only the first two terms of the series, we have 



* Even in the case of the moon the value of (p/i?) 2 is (rfijgp) 8 = 0.000277. 
In the case of the sun the quantity is far smaller than this. 
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On substituting this value of I/r 3 into equations (86-2) and inte- 
grating, we have 

M x =-^l I (y sin 0-1-3 cos 0) l+-g(2/cos0~-zsin$) dm 

yMi f r 3 r r 

= - 4^2 sin y^ + oos0 zcZm + -=j sin cos \(y*-z*}dm 

+ cos20 U^ |L 

The first two integrals are zero because the origin is at the center of 
mass of the earth. The last integral, a product of inertia, is zero 
because the coordinate axes are principal axes of inertia. The value 



of M x is thus 



f 
J 

sin 6 cos d f [> 2 + jr 2 - (x* + z 2 )] dm 

sin cos 1 f (a; 2 + y 2 ) dm - J (x 2 + z 2 ) dm 1 . 



The first integral on the right is the moment of inertia of the earth 
about the 2-axis (polar axis), and the second integral is the moment 
of inertia about the y-axis (an equatorial axis). Denoting these 
moments of inertia by C and A, respectively, we get 

J4=-^^sin0oos0(0-4). (86-3) 

The values of M y and M z as found from (86*2) are zero, because they 
are given in terms of the integrals \xdm y \xydm&nd p2<#m,eachof 

which is zero for the reasons stated above. The moment of the 
attraction is thus entirely about the cc-axis (axis of nodes) for the 
position of the earth here assumed & result that was evident from 
the outset, because of the position of the coordinate axes. Equation 
(86*3) gives the maximum moment that can be exerted by the sun. 
A similar expression holds for the maximum moment that can be 
exerted by the moon. 

Since M x as given by (86-3) is a vector pointing in the negative 
direction along the #-axis, both because of its sign and because of 
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the right-handed-screw rule for the vector product p x R, it may 
be written in the form 

M x =-to^(C-A)maO<nae. (86-4) 

JK 

The presence of i, sin 8 and cos 6 in (86-4) suggests that M x may be 

expressed in a more general form as a vector dot product times a 
vector cross-product. Hence we multiply both numerator and 
denominator of (86*4) by i? 2 and obtain 



J? 5 
SyJf 

' 



(C-A)(R cos A) ( - LR sin A) (see Fig. 67) 
(C--4)(R-k)(Rxk). (86-5) 



In order to find the moments about the axes of the moving 
trihedron and thence the differential equations of motion of the 
earth's axis, we must express the vector M in terms of unit vectors 
along the moving axes. This is accomplished by a process analogous 
to a transformation of coordinates. Referring to Fig. 68, let 
X l TL Z denote a set of fixed axes, with OX and OT l lying in the 
plane of the ecliptic and OZ l being perpendicular to that plane. 
Let XTZ denote the moving trihedron, OZ being the axis of 
spin (polar axis) of the earth, OX the line of intersection of the 
ecliptic plane by the plane of the earth's equator, and T lying in 
the equatorial plane and at right angles to OX and OZ. If we 
assume that the moving axes originally coincided with the fixed 
set, they can be brought to the position shown in Fig. 68 by first 
rotating through the angle i}r about OZ : , thus bringing OY^ to the 
position OM ; then rotating through the angle 6 about OX , thereby 
bringing the y- and z-axes to the positions shown. 

Let ii and j\ denote unit vectors along OX l and OT^ respectively, 
and let i, j, k denote unit vectors along the moving axes. Let j' 
denote a unit vector along OM. Then from the vector right triangles 
having unit vectors as hypotenuses, we have 



But j' 

Hence i 1 = icos^~(jcos5 ksin 6) smt/r 

(86-6) 
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And 



= i sin ijf + ( j cos 6 - k sin 6 ) cos ijr 

= i sin ^r -f j cos ^ cos 5 k cos ^r sin 5. (86-7) 

The position vector R changes in direction continuously as the 
earth travels around the sun. Let n denote the angular velocity of 
the earth in its orbit (assumed circular) around the sun. This is also 




the angular velocity of R, which after a time t will have swept over 
an angle nt. Hence if R coincided with OX l at time t = 0, its position 
at time t will be given by the vector equation 

R = i 1 S cos nt -f J-L E sin nt. 
Replacing i x and j x by their values as given in (86-6) and (86-7), weget 



= LR (cos nt cos i/r 4- sin nt sin ijr) 

+ IE cos 6 (sin nt cos ^ cos nt sin ijf) 

- kJB sin 6 (sin nt cos \Jr cos ft sin i/r) 
= iJ? cos (?tf - ^) + j-B cos (9 sin (nt i/r) - kR sin 5 sin (n^ i/r). 
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Hence by (5-6) 



and by (6-8) 

I j 



Rxk = 



JZcos(?i ^) JB cos 5 sin (nti/f) J?si 



00 I 

= IE cos 8 sin (w ifr) jJ2 cos (% ft) . 
Substituting into (86-5) these values of R . k and R x k, we have 



x [1R cos 6 sin (n^ i/r) jR cos (nt ^)] . 
Hence 14= - i -(0-4)sin6'cos5sm 2 (^-^) J (86-8) 






M x and M y can evidently be expressed in terms of double angles in 
the form 



M x = -((7-4)sin5cos5{l-cos2(^-^)}, (86-8 a) 



(86-9a) 

The angle nti/ria the longitude of the sun. It is also the angular 
distance of the earth from the line of nodes. If we put nti/r = a, 
equations (86-8) and (86-9) become 

in^cos^sin 2 ^ (86-86) 

-~- (86-96) 



These equations show at a glance that the precession moment is 
zero for a = and a = 180; that is, when the earth is at either of the 
equinoctial points. M x is always negative (always tending to make 
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the earth's axis perpendicular to the plane of the ecliptic) and has 
its greatest valne when a = 90 and a = 270; that is, at the summer 
and winter solstices. M y is zero at the equinoxes and solstices. It 
has its greatest values at a = 45, 135, 225 and 315. 
The mean value of M x from a = to a = 2?r is 



M x = - (C - A) sin cos (8640) 



(since the mean value of sin 2 a for this interval is |), whereas the 
mean value of M y over the same interval is zero. 

To find the differential equations of motion of the polar axis we 
use equations (25-5), because M Z =Q. From Fig. 68 we see at once 



o) x ^0 and a) y = TJr$m0. (86-11) 

Then from (25-5), (86-80) and (86-9a) we have 

^ 

~~ 



where O denotes the angular velocity of the daily rotation of the 
earth. Since the angular velocity of the moving trihedron is so 
exceedingly slow (the axis of nodes skews around at only 50".4 a 
year), the first two terms in the left-hand members of the above 
equations are negligible in comparison with the third terms. 
Hence on neglecting these terms and replacing OJ X and a> v by the 
values given in (86- 1 1 ), we have 



^- -Ztti^TF ooB*{l-ooB2(nMfl}, (86-12) 



d# SyM (C-A\ . . 0/ . /N /Q*I% 

T- "~^Ar>^ 7T~) sin a sin 2 (n^ 1 ^). (86'13) 

a# 2i2JB 8 .\ 6 / 

These are the differential equations for the precession and nutation 
of the earth's axis. 
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Equation (86' 12) shows that dt/rjdt is always negative and that J/T 
is therefore always decreasing. This means that the line of nodes is 
moving backward and thus the equinoxes have a retrograde 
motion. The rate of precession varies slightly and in a periodic 
manner owing to the term cos 2(nt ^). Its average value is found 
by omitting the periodic term cos2(n ^). Hence the average 



The nutation is seen from (86-13) to be periodic, with a period of 
TT radians per year. The spin axis of the earth thus traces on the 
celestial sphere an undulating curve whose mean arc distance is 
about 23 J from the pole of the ecliptic. The motion is a pseudo- 
regular precession and is very similar to that discussed in Arts. 37 
and 38. 

To find the numerical value of the annual precession due to the 
sun's attraction, we substitute the appropriate numerical data into 
(86*14). But since we do not know the numerical value of M, we 
must eliminate it. 

In the revolution of the earth around the sun, the sun's gravita- 
tional attraction on the earth is balanced by the centrifugal force 
acting on the earth. Hence we have the equation 

(86-15) 

where m denotes the mass of the earth and n denotes the angular 
velocity of revolution. From the familiar equation 6 = nt for 
circular motion, we get 

27r=nT or ^=-7=-, 

where T denotes the period of revolution (time for describing the 
angle 2n). Substituting this value for n into (86-5) and dividing 
through by m, we get 



Now substituting this value ofyM/R 3 into (86-14), we have 

3 
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We take the year as the unit of time. Then with reference to the 
fixed stars the angular velocity of the earth's rotation is 

} = 2n x 266 J radians per year.* 

Taking ^ = ^-^ (Tisserand), # = 23 27', ?=!, and substi- 
tuting these quantities into (86-16), we get 

3(4?r 2 ) 0.9174 

,, i - . ' , ^rz = - 0.00007725 radian per year. 

at / s 4-TT x 366^- x 305.6 

Dividing this result by 0.000004848, the radian value of 1*, we get 
0.00007725 15 ,. Qperyear- (86 . 1?) 



0.000004848" 

Up to this point we have assumed that the sun is the attracting 
body. The same analytical results hold in the case of the moon. In 
this case equation (86- 14) becomes 

tr\ A\ 

cosd, (86-18) 



where m 2 denotes the mass of the moon and r 2 is "the distance from 
the earth to the moon. 

To eliminate ym 2 /rf , we utilize the fact that the earth's attraction 
on the moon is balanced by the centrifugal force acting on the moon 
in its orbit. But here the attracting force is not simply ymm 2 /rf but 
is {y (m + m 2 ) m 2 }/rf , as is shown in works where the problem of two 
bodies is considered.f (The second expression was not used in 
dealing with the attraction of the sun because the mass of the earth 
is negligible in comparison with that of the sun.) Hence we have 

y(m-fm 2 )m 2 . 

^ - 2/ a = m 2 r 2 K|, (86-19) 



where n% denotes the angular velocity of revolution of the moon 
around the earth. IT 2 denotes the period of the moon's revolution 
around the earth, the equation for circular motion gives 



* The earth actually rotates about its axis 366i times a year, 365J times 
with respect to the sun and once in going around the sun. 

t See, for example, Appell's Traiti de Micanique Rationnelle, vol. i ( 19 1 9), 
pp. 413 and 416. 

16 so 
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from which n% %jrjT^ Hence (86-19) becomes 



(86-20) 



Since the mass of the earth is 81-5 times the mass of the moon, we 

have 

m=8L5m 2 . 

Substituting this value of m in (86-20), we have 
y(82.5) m 2 



ym 2 _ I f2nV 
Then substituting this value of ym 2 /rf into (86-18), we get 



The period of the moon's revolution around the earth is 27.32 
days. Hence 



Using this value for T 2 in (86-21), taking 

C-A_ I 

C 
and 6 = 23 27', we have 



1 . /365\ 2 0-9174 
x 



82.5 \27-32/ 305-6 

= - 0.0001674 radian per year 
= 34". 5 per year 

Hence the total precession due to both sun and moon is 
-15".9-34".5 = ~50".4 per year. 

* Here we use 365J instead of 366 J because we are concerned with the 
time required for the earth to revolve around the sun. The required time is 
365J days. 
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For a precession of 1 the required time is 3600/50.4 = 71.4 years. 
The time required for a complete revolution of the vernal equinox 
around the sun is 

360x60x60 rt ^ n 

= 25,700 years. 

87. Other astronomical applications 

(a) Retrogression of the Moon's Nodes. The moon's nodes are the 
points where its orbit intersects the plane of the ecliptic. They 
correspond to the earth's equinoctial points (equinoxes). They have 
a retrograde motion and make one revolution in the plane of the 
ecliptic in 18.6 years. 

The retrogression of the moon's nodes can be explained by 
gyroscopic principles. The interested reader will find the matter 
treated in Klein and Sommerfeld, Heft 3, in GreenhilTs Report of 
Gyroscopic Theory, p. 176, and in A. Gray's Gyrostatics and Rota- 
tional Motion, pp. 219-25. 

(b) Variation of Latitude. The earth's axis of rotation does not 
coincide exactly with its geometric axis of figure. The north 
geographical pole, for example, wanders in a circle of about 50 feet 
in diameter. This wandering of the pole causes a slight variation in 
the latitude of any place. 

If the earth is assumed to be an absolutely rigid body, the period 
of the pole variation is found by gyroscopic principles to be 

C 
- - = 305.6 days. 

But the earth is not an absolutely rigid body. It is slightly elastic, 
and its elasticity causes the period of the pole's wandering to be 
about 428 days. The interested reader will find the matter fully 
treated in Klein and Sommerfeld, in GreenhilTs Beport, and in 
Gray's Gyrostatics and Rotational Motion, Chapter 11. 
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APPENDIX 

Schuler Tuning of Gyroscopic Compasses 
and Gyroscopic Pendulums 

1. Introduction 

When gyroscopic compasses and gyroscopic pendulums are 
installed on ships, airplanes and other vehicles, the proper func- 
tioning of these gyroscopic instruments is liable to be disturbed by 
the motions of the vehicles on which they are installed. We saw in 
Art. 72, for example, that the rolling of a ship will produce a forced 
vibration in a gyroscopic pendulum suspended on the ship. Such 
facts and considerations led Max Schuler, of Germany, to a study 
of the conditions which must be satisfied if the performance of 
gyroscopic instruments is not to be disturbed by the motions of the 
vehicles on which the instruments are installed. He found that if the 
instruments are not to be affected by the motions of the vehicles, 
the vibration periods of the instruments must be equal to 2?r ^/(iJ/gr), 
where R denotes the radius of the earth and g is the acceleration of 
gravity. This vibration period is about 84 minutes. Schuler tuning 
consists in so designing gyroscopic instruments that their vibration 
periods shall be about 84 minutes. Such design is easily attainable 
for gyroscopic compasses, but is scarcely within the limits of 
practicability for gyroscopic pendulums. 

To get some idea of the principle underlying Schuler tuning, let 
us consider a gyroscopic pendulum suspended in a vehicle that is 
free to move in any direction over the earth's surface, and let E 
denote the earth radius drawn to the moving vehicle, as indicated 
in Fig. 69. We shall disregard the earth's rotation and consider the 
motion of a vehicle on a fixed earth.* 

* The earth rotates through 15 every hour, and this alone would tend to 
divert the spin axis of the gyro from the vertical by 15 per hour. But, 
because the earth's rotation is perfectly uniform, a suitable mechanism can 
compel the spin axis to rotate at the same rate and thus remain vertical when 
the point of support is at rest with respect to the earth. 
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Fig. 69 

If the point of suspension moves from A to B with velocity v, the 
earth radius drawn to the vehicle rotates through the angle 6 
with angular velocity dO/dt, or v say, so that v*=R(d0ldt)~Mv. 
Because of the tendency of a spinning gyro to maintain its spin 
axis in the same direction while the vehicle moves in any manner, 
the spin axis will not be vertical when the vehicle reaches -B, but 
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will have departed from the vertical by the angle 8. Hence if the 
spin axis is to be kept vertical while the vehicle is moving, it must 
be rotated simultaneously with the same angular velocity and 
angular acceleration as the earth radius drawn to the vehicle. The 
motion of the pendulum must therefore be a translation combined 
with rotation about the point of suspension. 

2. Schuler Tuning of a Gyroscopic Pendulum. 

As an aid in explaining the derivation of the condition for Schuler 
tuning in a gyropendulum, we have drawn in Fig. 70 a set of 
rectangular coordinate axes with origin at the point of suspension, 
but no use will be made of the axes except for reference in the 
explanation. Assume that the vehicle moves along the path A 
(Fig. 70) and let 

v = velocity of vehicle along AB (this path may be considered 
an arc of a great circle of the earth), 

R = earth radius drawn to the moving vehicle, 

v = angular velocity of earth radius drawn to the moving 
vehicle, 

a = acceleration of vehicle along AB, 
m = mass of pendulum, 
W = weight of pendulum = mg, 

I = length of pendulum. 

Then the required moment for rotating the spin axis in the #2-plane 
will be supplied by gravity if the pendulum is tilted away from the 
vertical by the small angle ft in the j/2~plane, as indicated in Fig. 70. 
The moment produced by this tilt is Wlsinfl about the re-axis, and 
it causes the pendulum to precess in the #2-plane with angular 
velocity v. Hence by equation (22-3) we have the relation 

W I sin ft = CW sin (90 - ft) 



f 1 T X O C^ 

from which tan/r=-^- =1^- 
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Since ft is small because of the large value of B, we may replace 
tan/9 by ft and get C(j)V 

P=WIR (1) 

as the equation of equilibrium of the pendulum for any speed v of 
the vehicle. The rate of change of y? from the equilibrium position is 



di 



G(i> dv Co) 
WlR~dt = WlR t 



(2) 



and thus varies with the acceleration of the vehicle. 






Now, when the vehicle moves in the #2-plane with acceleration a, 
it gives the point of suspension the same acceleration. But the 
pendulum, because of its inertia, resists such acceleration with the 
inertia reaction force (Wig) a acting through the center of gravity 
of the pendulum. This inertial reaction force exerts a moment 
(W/g)al about the point of suspension and in the #z-plane. This 
moment, because the vector representing it points in the direction 
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OF', causes the spin vector o) to precess in the t/z-plane with an 
angular velocity (d/}/dt) 2 toward OY f . Then by equation (22-4) we 
have the relation 

(3) 



If the pendulum is to remain vertical while the vehicle is moving 
in any manner, we must have 



Then on making these substitutions in (2) and (3), we get 

Co) _Wla 
Wm a ~~gCa) 






But, since C = mk*=(W/g) k 2 , where k is the radius of gyration of 
the pendulum about the spin axis, we get the simple relation 



as the condition for the Schuler tuning of a gyropendulum. 

We found in Art. 69 that the precession period of a gyropendulum 
is given by the relation T = 2n(Co}/ Wl) . Hence on replacing (Gw] Wl) 
from (4), we get 



= 27T I*. 

*J g 



(6) 



Taking .5 = 20,890,000 feet and ^ = 32.16 ft./sec. 2 , we find from (6) 

T = 5064 seconds 
= 84.4 minutes. 

On substituting in equation (1) the value of Cu/Wl from (4), we 
get 

e 75 



as the equilibrium angle for a Schuler tuned gyropendulum. 
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Numerical Example. To get an idea of the magnitude of ft as given 
by (7), assume that v = 60 miles per hour = 88 feet per second, 
5 = 20,890,000 feet, and g = 32.2 ft./sec. 2 . Then from (7) we get 

yff = 0.003392 radian 
= 11' 40*. 

If v = 300 miles per hour, we find /?== 58' 20", or 1 approximately. 
The required tilt is thus very small. 

In Art. 69 we found that when a gyropendulum is disturbed 
while its gyro is spinning rapidly, the pendulum bob will slowly 
process in a circular path, the supporting rod describing a conical 
surface. The results in the above example show that in a Schuler 
tuned pendulum the radius of the circular path is exceedingly small, 
being /? = 0.00339? for a vehicle speed of 60 miles per hour. 

In order to examine the practicability of Schuler tuning for gyro- 
pendulums, let us consider a numerical case. Assume that the rotor 
of a gyropendulum is 6 inches in diameter, that the rim is 2 inches 
wide and 1 inch thick, and that the web of the rotor is J inch thick. 
Then if 8 denotes the density of steel, we find 

moment of inertia of rotor = - 



* * 
mass ol rotor = - , 

144 



"96x864' 



nS _ 23 
~96x 864""" 144^576* 
If 

0) = 12,000 r.p.m. = 4007T radians per second, 

JR = 20,890,000 feet, 
g = 32. 16 ft./sec. 2 , 

we find from (5) I = 0.000308 x 2n feet 

= 0.023 inch. 

This short distance between the point of suspension and the center 
of gravity of the pendulum makes the construction of such a 
pendulum a rather difficult matter. 
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3. Schuler Tuning of the Gyroscopic Compass 

In Art. 56 we found that when a ship is sailing in a northerly 
direction with speed v and on a course making an angle with the 
north meridian, the gyrocompass points slightly to the west of 
north at an angle 7 given by the equation 

VGOB0 V 

( } 



where v n ( = v cos 6) denotes the northern component of the velocity. 
This equation gives the equilibrium position of the gyro axis for 
any speed and direction of the ship. The rate of change of 7 from 
the equilibrium position is 

*._ *. (2) 



JJOcosA dt .KQcosA' 

where a n denotes the northern acceleration of the point of suspen- 
sion of the compass and 7 increases westward. 

When the point of suspension is given a northward acceleration, 
the inertia of the pendulous weight w resists such acceleration by 
the inertia force (w/g) a n , and this exerts a moment r(w/g) a n about 
the point of suspension and in the plane of the meridian, where r 
is the moment arm of w. Since the vector representing this moment 
points toward the west, the gyro axis precesses toward the west at 
angular velocity (dy/dt) 2 in the plane elevated at the angle ft above 
the horizon (Fig. 46). Hence by equation (24-4) we have 



Then, if there is to be no disturbance in the indication of the 
compass, we must have 



Equating the values of these derivatives from (2) and (3), we get 



* __ W n 



.BQcosA gCa)' 

from which = i( ) cos A . (4) 

wr \g) 

This equation expresses the condition for Schuler tuning of a gyro- 
compass. 
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In Art. 52 we found that the long vibration period of a gyro- 
compass is given by the formula 

T^i 



' wQcosA" 
On substituting the value of wrQ. cos A from (4), we get 

*-%/ < 5 > 

which is the same period as found for the gyropendulum and there- 
fore gives a period of 84.4 minutes. 

Fortunately, it is possible and practicable to construct gyroscopic 
compasses that will satisfy condition (4), and modern gyro- 
compasses are so constructed. Gyrocompasses therefore indicate 
true north regardless of the motions of the ships on which they are 
installed. 

When equation (47-1) is combined with (4) above, the result is 

2 sin 2A 



or, since ft is very small, 

,_**** m 

which gives the elevation of the north end of the gyro axis when in 
its equilibrium position. For A = 45 we get 

ft= 0.00173 radian 
= 6' = 0. 1 , approximately. 
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